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December 26, 1963 



To the Members of the 1620 Users Group: 

It is with sincere regret that I must announce my resignation, effective during 
this meeting, as President of the Western Region. The executive council of 
the group, acting in accord with the by-laws, has appointed Paul Bickford to 
serve the remainder of the current term. 

Paul's appointment was made upon my recommendation, and I feel that his 
experience with the group will make him eminently qualified to continue its 
growth and activity. 

In turning the administration of the Western Region over to Paul, I can only 
say that I shall sincerely miss working with all of you, I express my thanks 
and appreciation to those of you, too numerous to mention individually, who 
have worked to make the meetings which I have conducted successful, and I 
extend my best wishes to the officers and members of the 1620 Users Group 
for the continued growth and success which you have enjoyed over the last 
four years. 



Sincerely, 
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1620 USERS GROUP MEETING 
TEMPE • ARIZONA 
DECEMBER 11*12*13* 1963 
ROSTER OF ATTENDEES 



USERS I NST ALLAT I ON 

GROUP NO* REPRESENTATIVE 



1118 NANCY PAQUIN 

U*S* PUBLIC HEALTH SERVICE 
ROCKV I LLE ♦ MARYLAND 

1118 G* ROBERT ORNDORFF 

U*S* PUBLIC HEALTH SERVICE 
ROCKV I LLE • MARYLAND 

1216 CARLIS TAYLOR 

A F R R I 

BETHESDA • MARYLAND 

1302 T* R* HOFFMAN 

UNION COLLEGE 
SCHENECTADY* N* Y* 

1334 DR* REIMUT WETTE 

M*D* ANDERSON HOSPITAL AND TUMOR INST* 
HOUSTON* TEXAS 

1346 EL IAS C* TON I AS 

ERDMAN AND ANTHONY - CONS* ENGRS* 
ROCHESTER ♦ NEW YORK 

3082 PAUL BICKFORD 

OKLAHOMA UNIV* MED* RESEARCH COMP* CTR* 
OKLAHOMA CITY* OKLA* 

3082 CAR A MITCHELL 

OKLAHOMA UNIV* MED* RESEARCH COMP* CTR* 
OKLAHOMA CITY* OKLA* 

3166 DR* HERMAN B« WEISSMAN 

UNIVERSITY OF ILLINOIS 
CHICAGO* ILLINOIS 

3175 HELEN LIGON 

BAYLOR UNIVERSITY 
WACO* TEXAS 

3182 ROBERT G* LANGE 

AUTOMATIC ELECTRIC LAB OR AT OR ICS* INC* 
NORTHLAKE* ILLINOIS 

3240 MELFORD E* MONSEES 

U*S* ARMY ENGR* DISTRICT 
KANSAS CITY* MISSOURI 



3261 



GREGORY J* SHANAHAN 

CECO STEEL PRODUCTS CORP* 

CICERO* ILLINOIS 



USERS I NSTALLAT ! ON 

GROUP NO* REPRESENT AT I VE 



3273 ROBERT C* BAB I ONE 

ACIC 

ST* LOUIS* MISSOURI 

3273 CHARLES WEISS 

USAF AERONAUTICAL CHART AND INFO CENTER 
ST* LOUIS* MISSOURI 

5001 R* C* WENRICK 

ACF INDUSTRIES* INC* 
ALBUQUERQUE • NEW MEXICO 

5001 G* J* REED 

ACF INDUSTRIES* INC* 
ALBUQUERQUE* NEW MEXICO 

5005 R* C* WEAVER 

BEAR CREEK MINING CO* 
SALT LAKE CITY* UTAH 

5014 WALTER DAVIS 

GENERAL DYNAMICS/ASTRONAUT ICS 
SAN DIEGO* CALIFORNIA 

5016 RICHARD tf* PUGSLEY 

COMPUTERMAT INC 
LOS ANGELES* CALIFORNIA 

5019 EDGAR M* BLIZZARD 
JET PROPULSION LAB 
PASADENA* CALIFORNIA 

5020 MARILYN DOIG 

COLORADO STATE UNIVERSITY 
FORT COLLINS* COLORADO 

5021 J* W« LAFON 

US ARMY ENGINEER 01 STRICT ALBUQUERQUE 
ALBUQUERQUE* NEW MEXICO 

5027 MARVIN RUBINSTEIN 

ELECTRO OPTICAL SYSTEMS* INC* 
PASADENA* CALIFORNIA 

5032 BOB MANNING 

GOODYEAR AEROSPACE 
LITCHFICLD PARK* ARIZONA 

5032 N* A* KUFFEL 

GOODYEAR AEROSPACE 
LITCHFIELD PARK* ARIZONA 

5032 J* MOSS 

GOODYEAR AEROSPACE 
LITCHFIELD PARK* ARIZONA 
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USERS 
GROUP NO* 



INSTALLATION 
REPRESENTATIVE 



5032 DAVID M« •HERREN 

GOODYEAR AEROSPACE CO 
LITCHFIELD PARK ARIZ 

5041 DAVID KEY 

MOTOROLA 

PHOENIX* ARIZONA 

5045 W* WILCOXSQN 

U*S* NAVAL CIVIL ENGR* LAB* 
PORT HUENEME • CAL IFORN 1 A 

5053 ROBERT L* SHUTT 

SACRAMENTO PEAK OBSERVATORY 
SUNSPOT* NEW MEXICO 

5053 FRANK BIRD 

SACRAMENTO PEAK OBSERVATORY 
SUNSPOT* NEW MEXICO 

5057 DAVID G# KITZINGER 

SANDS A CORP 

ALBUQUERQUE* NEW MEXICO 

5057 ELIZABETH L» FROST 

SAND I A CORP* 
ALBUQUERQUE* NEW MEXICO 

5056 BOB BABCOCK 

SUNDSTRAND AVIATION DENVER 
DENVER* COLORADO 

5064 ALLEN L* GRAVITT 
SIGNAL OIL AND GAS CO* 
LOS ANGELES* CALIFORNIA 

5065 ROBERT W* WILLSON 
SALT RIVER PROJECT 
PHOENIX* ARIZONA 

5065 ERNEST NICHOLS 

SALT RIVER PROJECT 
PHOENIX* ARIZONA 

5065 MAX A* HAYES 

SALT RIVER PROJECT 
PHOENIX* ARIZONA 

5078 DR« MORRIS J* GARBER 

UNIV* OF CALIFORNIA AT RIVERSIDE 
RIVERSIDE* CALIFORNIA 

5078 THOMAS M* LITTLE 

UNIVERSITY OF CALIFORNIA AT RIVERSIDE 
RIVERSIDE* CALIFORNIA 
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USERS INSTALLATION 
GROUP NO* REPRESENTATIVE 



5079 SAMUEL K» PR INGLE 

MAGNOLIA PIPE LINE CO* 
DALLAS* TEXAS 

5084 E* B. LOOP 

UNION OIL CO OF CALIF 
RODEO* CALIFORNIA 

5085 ROBERT D* MOFF1TT 

U*S* ARMY ENGR* DIVISION* N*P* 
PORTLAND* OREGON 

5089 DONALD J* MARTIN 

U.S. PUBLIC HEALTH SERVICE 
LAS VEGAS* NEVADA 

5089 DAVID L* BAER 

U*S* PUBLIC HEALTH SERVICE 
LAS VEGAS* NEVADA 

5095 HARRY D* REN IOC 
WEYERHAEUSER CO 
TACOMA* WASHINGTON 

5096 BOYD C* MORRIS 

U.S. BUREAU OF RECLAMATION 
SACRAMENTO* CALIFORNIA 

5096 R* BRUCE BROWNRIGG 

U*S* BUREAU OF RECLAMATION 
SACRAMENTO* CALIFORNIA 

5099 BEVERLY DOIG 

W S M R 

WHITE SANDS* NEW MEXICO 

5104 JOHN R. GUI NN 

TEXAS COLLEGE OF ARTS + INDUSTRIES 
KINGSVILLE* TEXAS 

5120 WILLIAM L* REUTER 

SO. DAKOTA SCHOOL OF MINES AND TECH* 
RAPID CITY* SOUTH DAKOTA 

5125 CHIN MO LEE 

UTAH POWER AND LIGHT CO. 
SALT LAKE CITY* UTAH 

5133 JOSE RAMIREZ 

MASON AND HANGER - SILAS MASON CO.. INC 
AMARILLO* TEXAS 

5139 R I CHARD A. HARRIS 

NORTH TEXAS STATE UNIVERSITY 
DENTON* TEXAS 
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USERS 
GROUP MO* 



INSTALLATION 
REPRESENTATIVE 



5144 JAMES J* STANLEY 

U S WEATHER BUREAU RFC 
SACRAMENTO t CALIFORNIA 

5146 BILL WOLLENHAUPT 

GOLDSTONE TRACKING STATION 
BARSTOW* CALIFORNIA 

5146 JOE SIBLEY 

GOLDSTONE TRACKING STATION 
BARSTOW* CALIFORNIA 

5150 SAM THOMPSON 

HALLIBURTON CO 
DUNCAN* OKLAHOMA 

5150 GEORGE A* LARCA0E 

HALLIBURTON CO 
DUNCAN* OKLAHOMA 

5158 A* GARLAND 

NECHES BUTANE PRODUCTS CO* 
PORT NECHES* TEXAS 

5165 CHARLES R* HEB8LE 

CORPS OF ENGINEERS 
WALLA WALLA* WASHINGTON 

5165 CECIL L* ASHLEY 

U S ARMY CORPS OF ENGINEER 
WALLA WALLA WASH 

5179 L* E* HARVEY 

FOOTHILL COLLEGE 

LOS ALTOS HILLS* CALIFORNIA 

5161 ROBERT R* WHITE 

LOS ANGELES DEPT* OF WATER AND POWER 
LOS ANGELES* CALIFORNIA 

5163 JOSEPH P* SNOW 

UNIVERSITY OF WYOMING 
LARAMIE* WYOMING 

5190 WILLIAM G* LANE 

CHI CO STATE COLLEGE 
CHI CO* CALIFORNIA 

5195 S* V* BURKS* JR* 

PITTSBURGH PLATE GLASS CO* - CHEM* DIV* 
CORPUS CHRIST I. TEXAS 

5199 CHARLES S* WALKER 

SCHOOL OF ENGINEERING 
TEMPE* ARIZONA 



USERS 
GROUP NO* 



INSTALLATION 
REPRESENT ATI VE 



5210 ANNA L* NIEHUES 

PALO ALTO UNIFIED SCHOOL DISTRICT 
PALO ALTO* CALIFORNIA 

5210 ARLINE K* KAPPHAHN 

PALO ALTO UNIFIED SCHOOL DISTRICT 
PALO ALTO* CALIFORNIA 

5211 KENNETH KRIEGE 

CALIFORNIA STATE POLYTECHNIC COLLEGE 
POMONA* CALIFORNIA 

5215 ROSEMARY PETERSEN 

UCLA - WESTERN DATA PROCESSING CENTER 
LOS ANGELES* CALIFORNIA 

7007 A* C« R» NEWBERRY 

UNIVERSITY OF ALBERTA CALGARY 
ALBERTA* CANADA 

SETH P* EVANS 
PHOENIX COLLEGE 
PHOENIX* ARIZONA 

WARREN BUXTON 
PHOENIX COLLEGE 
PHOENIX* ARIZONA 

JOHN P* MCCALLI STER 
U S WEATHER BUREAU 
FT» WORTH* TEXAS 

IBM W* H* DUKELOW 

IBM 

KANSAS CITY* MISSOURI 

IBM CHARLES £• BERRY 

IBM - WESTERN REGION OFFICE 
LOS ANGELES* CALIFORNIA 

IBM ROBERT A* EBERT 

IBM • WESTERN DATA PROCESSING CENTER 
LOS ANGELES* CALIFORNIA 

IBM JAMES E* MORGAN 

IBM - WESTERN REGION OFFICE 
LOS ANGELES* CALIFORNIA 

IBM BRUCE J* SOCKS 

IBM CORP* 
CHICAGO* ILLINOIS 

IBM ANGELO ARENA 

IBM 

WHITE PLAINS* N»Y* 



USERS 
GROUP NO* 



INSTALLATION 
REPRESENTATIVE 



IBM GERALD R* HOGSETT 

IBM - WESTERN REGION OFFICE 
LOS ANGELES* CALIFORNIA 
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WESTERN REGION 1620 USERS GROUP 
MINUTES OF THE SEVENTH MEETING 



The group assembled on the campus of Arizona State 
University in Tempe, Arizona, December 11, 12, and 13, 1963. 

All papers listed m the agenda were presented as 
scheduled except paper A-l, Generalized SPS Routines for 
Handling Simple Problems, by Thomas L. Yates, which was not 
given, paper A-2, Expotential and Sinusoidal Curve Pitting, 
by E. P. Hilar, which was given at Technical Session n G", and 
paper G-l, Ray Trace Program for a General Lens System, by 
D. H. O'Herren, which was given at Technical Session "A", 
Paper E-l, Network Analysis, was presented by Mr, Gerald Hogsett, 
IBM, Los Angeles, Paper C-2 was read by Bob Ebert, 

Copies of the abstracts and papers presented are 
enclosed with the exception of Paper G-l which will be included 
with minutes of a later meeting. 

The general meeting was opened on December 11 by 
President Bob Ebert and after opening remarks was turned over 
to Jim Morgan of IBM who introduced the other IBM representatives 
to the group. The sound -off session was moved up to this time 
to fill in the time for paper A-l, A summary of questions raised 
and answers to them follows the minutes of the business meeting. 

Business Meeting 

The business meeting was opened by Bob Ebert on 
December 12 and under old business, the request was made from 
the floor that the Roster of Members be arranged alphabetically 
in order of the name of the installation. The request was 
referred to Angelo Arena of IBM. 

The first item of new business was the resignation 
of Bob Ebert as president of the Western Region 1620 Users 
Group. This was necessary under the by-laws as he is no 
longer connected with a 1520 installation, A letter by Bob to 
the Users Group is included with the minutes of the meeting. 

The executive council, acting under the provisions 
of the by-laws, appointed Paul Bickford, Western Region Secre- 
tary, to complete the unexpired term of office of President, 
and he appointed Robert R, White to fill the vacant position 
of Secretary. 

With Paul now presiding as President the new business 
continued with the selection of Denver, Colorado, as the site, 
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and the third week in June, 1964, as the probable date of the 
next meeting. It was announced that the next fall meeting will 
be a joint meeting with the Midwestern Region, and will probably 
be held on the campus of Oklahoma University at Norman, Oklahoma, 
on November 9 or 22. The program committee will distribute the 
information regarding dates and agenda as early as possible, and 
they request abstracts of papers be sent in as soon as possible 
to aid them. A list of pre -registrants will be available early 
on the first day of the meeting. 

Paul then announced results of the Executive Council 
meeting at Pittsburgh and they are as follows : 

1. Jim Davidson is soliciting suggestions for 
methods for removing programs from the library. 
A criterion for value of the program and a means 
for the removal of inadequate programs is needed. 

2. The cost of running the users group has risen. 
Dues paid by each region to the National 1620 
Users Group were $.25 per registrant at each 
regional meeting. This has now been raised to 
$.50. Most of this increased cost has been in 
publication of the Newsletter. This must now 
be distributed to over 1,000 members while the 
average number of registrants at each meeting 
has not risen. It was stated that the regis- 
tration fees for the meetings may have to be 
increased . 

This concluded the new business and the business meet- 
ing was adjourned. 

After an excellent luncheon, Mr. Melford E # Monsees, 
ADP Co-ordinator, U.S. Army Engineer District, Kansas City, 
Missouri, spoke on "The Impact of Automation on Professional 
Engineering." A copy of the talk is enclosed. 

Two special interest groups met for discussion of 
mutual problems. The Civil Engineering group was presided over 
by Elias C. Tonias of Erdman and Anthony, Consulting Engineers, 
Rochester, New York, and minutes of the discussion are enclosed. 
The Power Utilities Engineer group was very informal but of 
benefit to all participants. 

Sound-Off Session and Comments 

It was announced that Version 2 of FORTRAN/FORMAT is 
available and has these advantages over Version 1 - 

1. Multiple Format Specifications (also in the 
Pre -Compiler) 
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2» Source decks are available (on magnetic tape) 

3, There are more subroutine options. 

Ninety hours or more are required for field installation 
of 1311 disk drives. 

Some operating difficulties with the disk were discussed 
and these points were brought out. 

1. In order to use a FORTRAN IID or SPS IID object 
deck the Monitor I must be on a working drive. 
This is because no loader or subroutines are 
punched . 

2. There is no current provision in Monitor I to 
load object decks to disk unless they are 
Monitor I compiled. It is therefore necessary 
to recompile all programs. 

There is a great deal of interest in FORTRAN/FORMAT. 
Requests were made for a disk version and a version to batch 
compile. It can now be changed to allow free form input similar 
to the first 1620 FORTRAN compiler. Chuck Berry of IBM Wilshire 
office can supply instructions for this change. 

Requests were made for the following from IBM: 

1. Relocatable I/O and arithmetic subroutines for 
FORTRAN so only those used would be called. 

2. A reduction in noise level on the Model 2 1622. 

3. A FORTRAN pre-compiler which would give indi- 
cation of the amount of storage the compiled 
program takes. 

4. A Tabulate command in FORTRAN. 

5* A loader for Systems Output Format decks punched 
under Monitor I control. 

6. FORTRAN processors to take advantage of extended 
machine capabilities. 

7. A Report Program Generator for the 1620-1443- 
1311 combination. 

8. COBOL for 1620 

9. ALGOL for 1620 
10. FORTRAN IV. 
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Some of these requests were answered in the comments 
session by Jim Morgan and the other IBM representatives. The 
questions answered are: 

1. Processors specifically designed for larger 
machine configurations are being investigated. 

2. There are no present plans for IBM written 
COBOL, ALGOL, or Report Program Generators, 
however, there is an ALGOL processor available 
from Southern Illinois University. 

3. The request for FORTRAN IV was noted and Jim 
asked for more opinions on the need for it. 

The new hardware announced for the 1620 since the last 
Western Region meeting is the CalComp Plotter in two models, 
binary capabilities for the Model 2 - 1620, and index registers 
for the Model 2 - 1620. 

It was noted that the FORTRAN II sine subroutine does 
not handle small angles correctly. IBM programming systems is 
working on the problem now. 

Angelo Arena discussed the program library and said 
that the KwTC index is still in the process of being developed 
and may be changed even further, so that it will be easier to 
find programs by number. It will be published every six months with 
supplements issued every month. It will be 3-hole punched. 
Angelo requested that when programs are ordered, unless the user 
is sure the program fits his need, he order only the documenta- 
tions. IBM can supply this quickly while program decks and tapes 
take longer. More descriptive titles for programs would also 
help this problem. 

There is now available from IBM a manual which lists 
the RPQ 1 s available (A26-5799) . Some of these mentioned were: 

1. The ability to punch one character on paper tape. 

2. A real time clock. The clock IBM is now install- 
ing on all equipment is an elapsed time clock 
only and is not addressable. 

3. An addressable IR-2. 

Version 2 of FORTRAN II will be available in the near 

future . 

The 1443 print commands are handled in the same manner 
as 7000 series machines with Column 1 being used for carriage 
control . 

In order to get IBM publications without fail, have 
the IBM Systems Engineer put the installation name on the SRL 

4 
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list. Mailing is then automatic. Be sure that the installation 
receives the 1620 Bibliography, A26-5692. 

The meeting was concluded with a demonstration of a 
Model 2 - 1620 with a 1311 disk file, on the evening of 
December 12, at the IBM Phoenix office and tutorial sessions 
on December 13 for MONITOR 1, FORTRAN II, and SPS 1620/1710. 



1620 USERS GROUP WESTERN MEETING 
PHOENIX, ARIZONA, WINTER 1963 
CIVIL ENGINEERING SESSION 
WEDNESDAY. DECEMBER 11. 1963. 3:30 PM 



The highlights of this session were discussions 
on the disc file, the use of Fortran vs. SPS in Civil 
Engineering programs, the introduction of COGO and the use 
of the CALCOMP plotter. 

It was the general consensus of the participants 
that SPS provided a better and faster object program than 
Fortran; Fortran however may be considered as the language 
for installations where speed in programming is preferred 
over speed in program runs. 

The use of COGO in Civil Engineering installations 
should be reserved to engineers who are neophytes in the 
field of electronic programming and should not by any means 
replace existing programs or prevent the development of 
special programs. 

Also of great interest was the use of the CALCOMP 
plotter now adopted by IBM. Successful programs have been 
written for contour plotting either with or without photo- 
grammetric equipment and in plotting cross sections and 
profiles. It has also successfully been used in structural 
design in such cases as in the plotting of influence lines. 
Another field of successful use is that of hydraulic engineering 
and hydrology. The use of other plotting devices such as the 
Digital Scale and the Wilde T7 has tremendously helped in 
reducing conventional type of work and in tightening the inter- 
relationship of photogrammetry, plotting and computing. 



EL IAS C. TONIAS 
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1620 USERS GROUP 
WESTERN REGION 

AGENDA FOR I963 WINTER MEETING 
ARIZONA STATE UNIVERSITY 
TEMPE, ARIZONA 

DECEMBER 11, 12, 13, 1963 



WEDNESDAY — DECEMBER 11, I963 
8:00 - 10:00 Late Registration 

9:00 Welcome and Opening remarks - R. Ebert, Regional President 

and IBM Representatives 

10:15 Coffee Break and informal "get acquainted" session 

11:00 Technical Session "A" 

A-l Generalized SPS Routines for Handling Simple 
Problems . 

Thomas L. Yates, Director, Statistics Computing 
Lab,, Oregon State University, Corvallis, Ore. 

A-2 Exponential and Sinusoidal Curve Pitting. 

E. P. Hilar, Goodyear Aerospace, Litchfield Park, 
Ariz. 

12:00 Break for Lunch 

1:30 Technical Session "B" 

B-l Automated Design Engineering. 

W. W. Rogers, IBM, Los Angeles, Calif. 

B-2 A Payroll and Labor Distribution Program Package. 
Elias C. Tonias, Erdman and Anthony, Consulting 
Engrs., Rochester, N.Y. Richard C. Devereaux, 
IBM, Rochester, N.Y. 

B-3 The SPIRE System - Salaried Personnel Information 
REtrieval . 

(Jary J. Reed, Proj. Engr., ACF Industries, 
Albuquerque, New Mexico 

3:00 Coffee Break 

3:30 Technical Session "C" 

C-l A 1620 Program for Minimization of Boolean Functions 
Expressed as Sums of Minterms. 
Thomas R. Hoffman, Prof, of Elect. Engr., Union 
College, Schenectady, N.Y. 
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Agenda - Tempe - Page 2 



WEDNESDAY — DECEMBER 11, I963 (Cont'd.) 

C-2 Critical Speed, Stress, and Bearing Reaction Cal- 
culations for a General Shaft, Using Numerical 
Integration. 

Ralph B. Bates, Mgr. of Engr. Computing, Industrial 
Div. of American Standard, Detroit, Mich. 

C-3 Three Dimensional Surface Pit. 

David G. Kitzinger, ACP Industries, Albuquerque, 
New Mexico 

5:00 Adjournment of Day f s Sessions. 

8:00 p.m. New Users Meeting, to- 3?e followed by Sound-Off 
Session at Approximately 8:30. 



THURSDAY -- DECEMBER 12, 1963 

9:00 Technical Session "D" 

D-l Maximum Likelihood Resolution of Two Mixed Normal 
Distributions. 

Reimut Wette, D.Sc., Asst. Biometrician, The 
University of Texas, M.D. Anderson Hospital and 
Tumor Inst., Houston, Texas 

D-2 Comparison of Two Methods of Finding Significant 
Contributors in Multiple Regression. 
M. J, Garber, Director, Computing Ctr., University 
of California, Riverside, Calif. 

10:15 Coffee Break 

10:45 Technical Session "E" 

E-l Network Analysis. 

H. N. Tyson, Jr., IBM, Los Angeles, Calif. 

11:15 BUSINESS MEETING 

12:00 LUNCHEON "The Impact of Automation on the Professional 

Engineer" 

Melford E. Monsees, ADP Co-ordinator, U.S. 
Army, Corps of Engineers, Kansas City, Mo. 

1:30 Technical Session "P n 

P-l An Integrated Earth Work System. 

Cecil L, Ashley, ADP Co-ordinator, U.S. Army 
Engineering District, Walla Walla, Wash. 

F-2 Hydro-System, A Daily Operation Analysis. 

Charles R. Hebble, Civil Engr., U.S. Army Corps 
of Engineers, Walla Walla, Wash. 
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THURSDAY — DECEMBER 12, 1963 (Cont'd. ) 

3:00 Coffee Break 

3:30 Technical Session "G" 

G-l Ray Trace Program for a General Lens System, 

D. H. O'Herren, Goodyear Aerospace, Litchfield 
Park, Arizona 

4:00 IBM Reports and Discussion of Sound-Off Session, 
5:00 Adjournment of Day's Sessions. 

8:00 p.m. Demonstration of the 1620, Model II, with 1311 
Disc Pack at the Phoenix IBM Branch Office. 
(Transportation will be arranged) . 



FRIDAY DECEMBER 13, 1963 

9:00 Workshop Session "A" 
A-l MONITOR I 

1:00 Workshop Session "B" 
B-l FORTRAN II 
B-2 SPS 1620/1710 

4:00 Conclusion of Workshop Sessions, and Final Adjournment 
of Meeting. 
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1620 USERS GROUP 
WESTERN REGION 

ABSTRACTS FOR I963 WINTER MEETING 
ARIZONA STATE UNIVERSITY 
TEMPE, ARIZONA 

DECEMBER 11, 12, 13, 1963 



NO. TITLE, AUTHOR, ABSTRACT 

A-l GENERALIZED SPS ROUTINES FOR HANDLING SOME SIMPLE PROBLEMS . 

Thomas L. Yates, Director, Statistics Computing Lab., ~~ 
Oregon State University, Corvallis, Ore. 

Two-way frequency distributions; output editing and 
formatting: card reproducing. (No further abstract 
available . ) 

A-2 EXPONENTIAL AND SINUSOIDAL CURVE FITTING . 

E. P. Hilar, Goodyear Aerospace Corp., Litchfield Park, 
Arizona. 

A specified number of exponentials are fitted to 
given equally spaced data points. The frequencies of 
the exponentials can be a mixture of real frequencies 
and complex conjugate pairs. The least mean square 
error criterion is used to find both the frequencies 
and the coefficients of the exponentials. The program 
is written in FORTRAN language. The limitations of the 
program will be discussed. 

B-l AUTOMATED DESIGN ENGINEERING . 

W. W. Rogers, IBM Corp., Los Angeles, California 

A presentation of Automated Design Engineering and 
how to achieve it. The role of a computer, ability to 
capture design logic, use of decision tables, comparison 
to manual design methods, and considerations involved 
with establishing Automated Design Engineering will be 
discussed. 

B-2 A PAYROLL AND LABOR DISTRIBUTION PROGRAM PACKAGE . 

Elias C. Tonias, Head of Data Proc. Dept., Erdman & 
Anthony, Consulting Engineers, Rochester, N.Y., and 
Richard C. Devereaux, IBM Corp., Rochester, N. Y. 

The objective of this paper is to demonstrate how 
some free lo20 time may be utilized in a relatively 
small scientific or engineering installation through 
the use of a package of commercial programs. The 
successful operation of such a program package since 
the first of this year (1963) has helped to justify 
the installation of a 1620 in another firm. This 
package, designed for use with the basic 20K 1620 

(Cont'd.) 
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Computer, with paper tape 1/0 and without any peripheral ( 
equipment, produces the payroll reports, and complete 
labor cost distribution and breakdown reports. Written 
in FORTRAN for easy maintenance, the ideas from this 
program package might prove a worthwhile tool in justi- 
fying the installation, or in increasing the production 
ratio of a 1620 in a small scientific or engineering 
account. With few or even no alterations, parts of the 
package may be used to handle other various time and 
cost distributions. 

B-3 THE "SPIRE" SYSTEM - SALARIED PERSONNEL INFORMATION RETRIEVAL . 

Gary J. Reed, Project~~Engineer, ACF Industries, AlbuqUerque, 
New Mexico 

The Albuquerque Division of ACF Industries has devel- 
oped a powerful management tool for use in personnel 
administration. This tool is the SPIRE system, utilizing 
an IBM 1620. 

The system involves over 30 programs designed to 
utilize information from a master file of information 
maintained on magnetic tape. This master file contains 
over 300,000 different items of information. By Judicious 
choice of information blocking, search time has been kept 
to a minimum for all applications. 

The SPIRE system current applications include: in- 
plant recruitment, providing resumes of personnel qualified 
to fill vacant positions; preparation of quarterly reports 
containing information on employees eligible for merit 
considerations; creation of summary reports on salary 
increases for any time period; man power inventories; 
projections for salary budgeting; preparation of data for 
salary surveys; current salary status reports; and many 
other similar reports. 

The information contained in the SPIRE system has 
proved to be complete for all applications thus far. The 
approach and system layout have provided an economic way 
of producing information that no hand techniques could 
supply at any cost. In the applications where the results 
could be obtained by clerical methods, enormous cost 
savings have resulted. 

The SPIRE system has been operational since February 
1963. The use of existing programs and the creation of 
new applications consistently increase reflecting over- 
whelming management acceptance. 

C-1 A 1620 PROGRAM FOR MINIMIZATION OF BOOLEAN FUNCTIONS 3 

mftgasfip as sums og mintms . 

Thomas R. Hoffman, Prof, of Elect. Engr., Union College, 
Schenectady, N.Y. 

A FORTRAN program implements a major portion of the 
Boolean minimization problem by reducing a sum of minterms 
to a logically equivalent set of prime implicants. 

Minterms are entered as 3-digit octal -coded fixed A 
point numbers. With the help of a special table, digit ^ 
comparisons at the octal level reveal terms combinable 
according to the Boolean identity XA + X# = X. Systematic 

(Cont'd.) 
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determination of all possible combinations of this type, 
coupled with a bookkeeping system that keeps track of 
the different subsets produced, leads to the desired 
result. Prime implicants are printed out in octal code, 
to complete the processing. 

The program can handle minterms having as many as 
nine variables, although the 1620 memory capacity (20K) 
may be exceeded in problems involving both large numbers 
of variables and long lists of minterms. 

C-2 CRITICAL SPEED, STRESS, AND BEARING REACTION CALCULATIONS 

FOR A gSNgftAL SHAFT, PSDta NUMERICAL MSgtUITOC — 

Ralph B. Bates, Mgr. of Engr. Computing, Industrial Div. 
of American Standard, Detroit, Michigan 

Critical speed, stresses, and bearing reactions can 
be calculated on a digital computer for a general two- 
bearing shaft, any number of cross sections, variable 
loading, and any length. 

Besides eliminating the tedious labor of the calcu- 
lations, the computer provides flexibility. A number of 
calculations may be rapidly made to optimize design, or 
to check out application variations on a standard design. 

The method is illustrated in detail by an example 
calculation. Briefly it consists of dividing the shaft 
into increments of length, determining the load and shaft 
moments of inertia in each ire rement and the computer 
calculates critical speed, stresses, and bearing reactions. 

C-3 M-151, THREE-DIMENSIONAL SURFACE FIT . 

David G. Kitzinger, ACF Industries, Albuquerque, New Mexico 

This code uses multiple interpolation techniques in 
combination with extensive transformation of variables to 
effect accurate fits of most smooth functions of three 
variables. With limited experience in selecting the form 
of curve fits in two dimensions, the programmer can fit 
complex three-dimensional surfaces. Output is designed 
to facilitate successive better approximations to the 
function in terms of additional transformation of variables. 
Second and third order fitting of functions is aided by a 
statistical error analysis. Two magnetic tapes and a 60K 
memory are required, although modifications can be easily 
made for smaller machines. 

D-l MAXIMUM LIKELIHOOD RESOLUTION OF TWO MIXED NORMAL 

DISTRIBUTIONS" : 

Reimut Wette, D.Sc, Asst. Biometrician, Univ. of Texas, 
M.D. Anderson Hospital and Tumor Institute, Houston, Texas 

Iterative estimation of the five parameters from a 
sample taken from two normal distributions by the method 
of maximum liklihood seems preferable over the method of 
moments, because generalization to more than two parent 
distributions appears easier. Two approaches to estimate 
the maximum likelihood information matrix were abandoned 
in favor of the third, which cut down on computing time 

(Cont«d.) 
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(1620 FORTRAN II) for the complete estimation procedure 
considerably. Computing time is, aside from variations 
due to behavior of the procedure depending on goodness of 
initial estimates and structure of the sample, directly- 
proportional to the size of the sample. Therefore, a 
grouping-of-data program was developed specifically for 
this problem, and large samples, preceding the estima- 
tion procedure and reducing computing time to reasonable 
limits. 

COMPARISON OP TWO METHODS OP FINDING SIGNIFICANT CON- 
TRIBUTORS IN MULTIPLE REGRESSION. 

M. J. Garber, Director, Computing Ctr., University of 
California, Riverside, Calif. 
(No abstract available.) 

NETWORK ANALYSIS . 

H. N. Tyson, Jr., IBM Corp., Los Angeles, Calif. 

Discussion will center about the general technique 
of network analysis, and the capabilities of a system of 
the 1620 programs utilizing this technique for the analysis 
of electronics circuits. The programs operate under MONITOR 
I on a 40K 1620. The capabilities allow for AC, DC, and 
transient analysis. 

AN INTEGRATED EARTH-WORK SYSTEM . 

Cecil L. Ashley, ADP Co-ordinator, U.S. Army Engineering 
District, Walla Walla, Washington 
(No Abstract Available.) 

HYDRO SYSTEM DAILY OPERATION ANALYSIS 

C. E. Hildebrand, L. A. Duns tan, C. R. Hebble, R. D. Moffitt, 
U.S. Army Engr. District, Walla Walla, Washington 

The program is a mathematical model of a system of 
hydro-electric projects. It is a general program appli- 
cable to any river system and scheme of development. It 
is capable of accurately simulating the hour-by-hour 
operation of a hydro system for as long a real time period 
as desired. It will determine the effects produced by 
existing hydro stations in regards to reservoir levels, 
river stages, and alternative distributions of system 
load among a group of hydraulically and electrically 
integrated hydro stations. The program makes it possible 
to determine the operating characteristics of planned 
future projects in regards to backwater encroachment on 
upstream reservoirs, pondage requirements, the effect of 
peaking discharges on downstream river stages and reser- 
voirs, and effects of added power installations. Two 
versions of the program exist: One for an IBM 1620 with 
40K memory, and one for an IBM 1620 System with 60K in 
the IBM 1620 and 4K in the IBM 1401. 
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G-l RAY TRACE PROGRAM FOR A GENERAL LENS SYSTEM . 

D. H. ! Herren, Goodyear Aerospace Corp., Litchfield 
Park, Arizona 

This program, written in FORTRAN with Format, and 
requiring 40K core storage, is designed to trace rays 
through complex lens systems, outputting intersection 
points at selected surfaces in the system. Cylindrical, 
conic, and toric surfaces can be handled and can be arbi- 
trarily oriented with respect to the optical axis. Lens 
surfaces can be virtually any contour describably by 
second degree (or less) equations or at least divisable 
into sections which can be so described. A sample ray 
tracing problem is given to illustrate the use of the 
program. 



The compiler workshop sessions to be held on Friday will 
assume that those attending will have a working knowledge of the 
externals of the programming systems (i.e. how to write programs 
in FORTRAN or SPS, etc.). The purpose of these sessions will 
be to present some of the "Hows and "Whys" of the internal 
workings of the compilers. The session on MONITOR I will include 
more information on using the monitor, as well as delving into 
some of the internals. 
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SOCIAL AND ECONOMIC ASPECTS OF AUTOMATION 



There are many differing views regarding the impact that the elec- 
tronic computer will have on the development of our social and economic 
lives, but most of us can agree that the potentialities of automatic 
data processing (ADP) are limited only by the boundaries of our individ- 
ual imaginations. The outlook, therefore, is open to wide-ranging 
speculations. 

The truly great impact of digital computation will be a dramatic 
speedup in the rate of technological progress and concurrently in the 
evolution of our social and economic lives. 

Electronic data processing is one of the most powerful catalysts of 
technological development yet discovered. This is so because ADP has 
the ability of extending the capabilities of man's intellect. 

The human mind is the most powerful, most versatile, most useful 
natural gift bestowed upon man by his Creator. Any instrument that can 
substantially increase its capabilities is certain to have a profound 
effect on our future development. 

Man succeeded in building our technology to an extraordinarily sophist 
cated level during World War II without an aid to his .intellectual capa- 
bility even remotely approximating the power of ADP. Now, with the aid 
of computers that increase the productivity of his intellect in many 
areas by factors running into the tens of thousands, we are certain to see 
significant advances in the tempo of technological, economic and social 
progress. 



Another broad area of impact will be in fields of business, indus- 
try and communications. We are all familiar with the progress ADP has 
already made in automating the business office. And computers are now 
being used increasingly in manufacturing— not only in process control 
applications, but also as an aid to the efficient management of the over- 
all manufacturing operation. 

We are entering a battle for our very survival in the market places 
of the world. Much emphasis is being placed on economic development. 
Technological progress, rising productivity and ascending standard of 
living are the true sources of economic strength. They are vital to 
national survival in today's competitive world. 

As we move forward, we will encounter the problems always inherent 
with social and economic change. These problems should in no way warrant 
artificial restrictions on technological development for this is vital 
to the success of any business in a free economy. 

It may well turn out that the efforts of the new technology will 
be far more lastingly felt in its impact on many of the traditional prin- 
ciples and practices of management. Many traditional personnel practices 
are obviously going to be automated or abolished and various leader 
groups will change in power and prestige. 

While there are no precise means as yet of measuring the speed of 
technological change, it is reasonable to assume that by the mid-1960' s, 
as those born during the Second World War establish families and the 
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challenge of foreign competition becomes more intense, the present rate 
of change will increase. At the present time, productivity increases in 
the nonfarm sector amount to roughly 2.5 percent a year. Even a moderate 
speedup in this rate would mean that by 1990, a relatively short span of 
less than 30 years, industry could double its production with the same 
labor as it employs today. 1 

THE ENGINEERING EMPLOYMENT SITUATION 2 

Engineering, the second largest professional occupation, is exceeded 
in size only by teaching; for men, it is the largest profession. The 
approximately 875,000 engineers in the United States in mid-1960 have 
made major contributions to the design, construction, and efficient 
utilization of the machines, equipment, roads, and buildings used by 
the Nation's 180 million people. Engineers provide technical, and fre- 
quently, managerial leadership in industry and Government. They develop 
new products and processes, design many types of machines and structures, 
and contribute in countless other ways to the technological progress of 
the country and to the national defense. 

The outlook is for continued rapid expansion of the engineering 
profession. Engineering has been one of the fastest growing professional 

* Joseph A. Raffaele, "Automation and the Coming Difussion of Power 
in Industry," Personnel , (May- June 1962), 30. 

2 

U. S. Bureau of Labor Statistics, Employment Outlook for Engineers , 
(Washington: U. S. Government Printing Office, 1961), pp. 101, 103-104. 
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occupations in the United States in the past 50 years, and there is every 
indication that the demand for engineers will continue to grow. As in 
recent years, there will probably be a particular need for engineers with 
advanced degreas to teach and to do research. 

Some of the major factors expected to raise the demand for engineer- 
ing personnel are: Continued high levels of Government spending for 
defense, accentuated by the increasingly large amount of engineering 
time necessary for the development of modern weapons; growth of popula- 
tion and expansion of industry; increasing complexity of industrial 
technology, as such the trend toward automation of industrial manufacturing 
processes; and further growth in expenditures for research and development. 
In particular, the large sums spent for research and development in recent 
years by both industry and Government — total research and development 
expenditures in the United States amounted to more than $13 billion in 
1960-61 — have broadened existing areas of employment for engineers and 
opened up new ones, such as those concerned with computers, missiles, and 
nuclear energy. As scientific frontiers are extended, more areas of work 
for engineers will be provided. 
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THE EMERGENCE OF AUTOMATION IN ENGINEERING 



That the systems concept is having a reaction in the engineering 
field is generally recognized. The question is how much. The change in 
a lot of features surrounding the engineering profession is beginning to 
assume large proportions and that rate of change is increasing rather 
than leveling off. 

For instance, we are told that the amount of engineering information 
and scientific information which is directly relevant to engineering 
problems has doubled within the last fifteen years. Even now the volume 
of knowledge related to engineering is so great that no one man can 
possibly know it all, even though he does nothing but study from boyhood 
to senility. For us engineers that means two things. We must specialize 
more than was formerly necessary or desirable, and second, we must be 
diligent students throughout our active lives. If we fail to study and 
keep abreast of the developments which directly apply to our chosen 
fields of engineering, we shall quickly become back numbers and soon there- 
after become useless to an advancing civilization. 

As a civil engineer, I am more familiar with the developments in that 
field than in some of the others. I recall the observations I have made 
and the discussions I have had with various men over the past 25 years, 
relative to the extreme reluctance of engineers to adapt labor-saving 
devices in their own work . Our fellow engineers were accomplishing much 
in the industrial world in devising and perfecting labor-saving equipment, 
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but we civil engineers were extremely slow in demanding any sort of 
machinery or equipment that would make our work easier o 

For example, up to a few years ago, we were using the same methods 
of surveying that had been devised in ancient times. The transit was a 
little better than the old surveyor's compass, but the long and tedious 
process of making ground surveys was basically unchanged for 3,000 years. 
Now electronic measuring devices for survey parties are becoming standard 
equipment. 

The electronic computer is an item which cuts across all the fields 
of engineering and many of the areas of science as well. They have been 
in use now for only about 15 years, but in that time they have increased 
the computational ability of mankind one million times. The industry has 
grown in 15 years from zero to the production of a billion dollars worth 
of equipment during the year 1960. Since that date there has been an 
upsurge in the number of manufacturers of electronic computers. Competi- 
tion has become keener than it was a few years ago and we can confidently 
look forward to a greater variety of computers, composed of more depend- 
able pieces of equipment, at a cheaper price than they command today. 3 

The electronic computer systems have eliminated much of the routine 
drudgery that has long been the bane of an engineering office. Also, they 
have eliminated the need for those men who were fitted for nothing more 
than routine work. There will be less and less need for the man who can 
only run a calculator or a slide rule, and who has been noted in an office 

3Murray A. Wilson, "Change or Progress," American Engineer , (June, 
1962), 29. 
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primarily because he could remember all the formulae that applied to his 
line of work, able to recall and use them accurately as the occasion 
demanded. This work, the machine can do better and much more rapidly. 

The other feature is that the computer can solve the basic, theoret- 
ical equation and eliminate the necessity for the approximations and short- 
cuts that have been in common use in so many phases of design, simply 
because the basic formulae were so long and complicated that the process of 
solving them longhand took so long that no one could afford to use them. 

On the other hand, the new facilities will put a premium on imagination 
and ingenuity . These have always been desirable qualities in an engineer, 
but in the past a great many engineers have been kept gainfully employed 
on jobs that require little of either. This situation is changing and 
the prizes in the future are going to those men who are well endowed with 
these two important attributes. 

EFFECT OF AUTOMATION ON ENGINEERING EDUCATION 

Actual and tangible changes in the practice of the art and science 
of engineering are being reflected in our schools and are under constant 
discussion in publications, so it seems probable that we have created a 
new concept of the profession of engineering, or if you please, a new 
image of the engineer. The curricula in our colleges are in a state of 
fluidity with considerable differences in their means of meeting the 
challenge of the changing conditions. On one thing all seem to agree — 
that is, that the engineer of tomorrow and the day after, will need to be 
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more thoroughly grounded in the basic physical sciences than was thought 
necessary for his predecessor. Others have felt that an equally important 
need of the coming engineers is a knowledge of the past such as is gained 
through a study of the "humanities." In an attempt to make room for these 
additional courses in a curriculum already overcrowded, something had to 
give, so the subprofessional subjects such as shop practice, material 
testing, laboratory work of various kinds, surveying and similar courses 
are being eliminated. The logical justification of this elimination is 
that these functions are actually to be performed by subprofessional men 
anyway, or as we propose to call them, the engineering technicians. 

Modern problems of engineering are no respecters of traditional 
boundaries between the specialities. Accordingly, changes are being made 
by schools of engineering. It is understood that Dr. Keith Glennan has 
made broad moves which will go far to establish interdisciplinary approaches 
in engineering education at Case Institute of Technology. At the under- 
graduate level he has consolidated the departments of chemical, civil, 
electrical and mechanical engineering into a single administrative unit — 
the Engineering Division. The engineering faculty are re-grouping in a 
natural way according to their common professional interests such as 
systems, design, energy conversion, materials, information processing 
and other emerging fields. 

Degrees in electrical, mechanical, civil and chemical engineering 
will continue to be granted, but a new degree at Case — probably named 
Bachelor of Science in Engineering — will be offered. It will give the 
student the opportunity to plan his elective program — with faculty 



advice — to suit his career interests. Such a program can be designed 
to lead more effectively into advanced work. 

There is little doubt that the four-year graduate will continue to 
play an important role in industry, but there is an ever increasing need 
for the engineer with the depth of knowledge and experience produced by 
work at the advanced graduate level. 

Also, Dean B. R. Teare, Jr. of the College of Engineering and Science, 
Carnegie Institute of Technology, said recently in a letter to me — 
"The electronic computer has certainly made an impact on the individual 
courses in our engineering programs but it has not yet been the reason 
for extensive curriculum changes. Some engineering departments, pressed 
with the lack of time in a four-year program, have had to decide between 
continuing courses in engineering graphics and courses in computer logic."^ 

The electronic computer is also having its impact on the curriculum 
at M.I.T. In a letter I received last January from Dean Gordon Brown's 
office some of the directions in which the computer was leading were out- 
lined. ^ For example the inclusion in the curriculum at M.I.T. of the 
following subjects: 

Digital Computer Programming Systems 

Mathematical Methods in Civil Engineering 

Digital System Application. 



^B. R. Teare, Jr., Carnegie Institute of Technology, in letter to 
author dated January 3, 1963. 

^Gordon Baty, Adm. Asst. to Dean of school of Engineering, M.I.T., 
in letter to author dated January 18, 1963. 

10 



3<* 



Other related subjects have been included iri the curriculum, but it 
pointed out that the automatic computations is not presented as an end in 
itself, but as another tool of analysis in the engineer's kit. The subjects 
merely focuses upon the techniques available to the engineer for exploiting 
the power of electronic computations. 

Dean Brown's office also advised me that the impact of the computer 
upon research activities has been enormous. A copy of the semiannual 
Report, available from the M I«T» Computation Center, can give you some 
idea of its magnitude. Yet, however important these computer-related 
activities have become to the School, there is a danger involved in 
attributing any of them simply to the availability of computer technology. 
For this is only one of the influences which have converged to create many 
of our most exciting research projects and subject offerings. Others include 
new methods in statistics and operation analysis, Systems analysis and 
synthesis, theory of learning, and the information technologies 

TANGIBLE BENEFITS FROM AN ELECTRONIC COMPUTER SYSTEM 

First hand knowledge of the impact of automation on professional engi- 
neering has been obtained as a member of the staff of the U. S. Army 
Engineer District, Kansas City, Missouri. This Corps of Engineers office 
has civil works engineering and construction in States of Kansas, Missouri, 
Nebraska, Iowa, and Colorado, and military engineering and construction in 
the States of Kansas and Missouri. The total work of this district averages 
about 70 million dollars per year and includes the design and construction 
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of large multipurpose dams, levees, floodwalls, channel Improvement works, 
pumping plants, and necessary utilities, highways, railroads, bridges, etc., 
required to be relocated in connection with flood control projects. Also, 
included in the assigned work of the district is the construction of 
military facilities and structures for the Army and Air Force. 

Computer facilities for engineering applications have been available 
since January 1958. The initial computer was the Burroughs E-102, which, 
although of limited capability, was used until replaced in December 1960 
with an IBM 1620 paper-tape system. Utilization of the paper-tape system 
for engineering applications increased rapidly and in January 1962, the 
system was augmented to provide for high speed punched card input-output. 
During this period, the utilization of electronic computer systems has 
provided tangible benefits through reduced cost of construction, savings 
in engineering and clerical manpower, and by providing a superior end- 
product or flood control structure. To date the major effort in implementing 
ADPS procedures has been directed toward high- benefit engineering applica- 
tions. Only preliminary phases of planned implementation of ADPS in areas 
of personnel administration, property accounting, real estate activities 
and fiscal responsibilities has been possible. As of December 1963, over 
75 computer programs were being used in the fields of mechanical engineering, 
structural engineering, hydraulic engineering, hydrology, reservoir regula- 
tion and earthwork and soil mechanics. 

The currently installed system in the Kansas City District Office has 
provided the following: 
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a. Improvement of the engineering and design product. 

b. More timely and accurate data which is fully responsive to the 
needs of management and engineering, including data previously not economi- 
cally obtainable. 

c. Savings in costs by maintaining continued evaluation and balance 
of equipment and personnel. 

d. Savings in costs and engineering manpower by application of ADP 
principles. 

e. Simplified and reduced manual data handling and eliminated dupli- 
cation of files, reports, and entrance of source data. 

Further benefits are being obtained through realization of the follow- 
ing objectives: 

a. Expansion of hydropower and pumped storage-power study programs to 
be used in connection with reservoir regulation, hydrology and hydraulic 
engineering programs used for design, construction, and operation of multi- 
purpose projects. 

b. Expansion and refinement of structural design analysis programs 
used for preliminary and final design of various civil engineering projects. 

c. Expansion and refinement of earthwork and soils mechanics programs 
used for quantity computations and stability analysis of large and small 
embankments. 

d. Implementation and expansion of Critical Path Scheduling techniques 
to be used for coordination of construction activities as well as coordina- 
tion and scheduling of engineering and design programs. 
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e. Refinement of personnel administration and reporting procedures. 

f. Implementation of one additional phase of the engineering budget 
management data to eliminate manual posting and to provide more accurate 
readily accessible data for management and estimating purposes. 

g. Exploitation of the principle of "management by exception" through 
the potentialities of data processing equipment and techniques by continued 
education and training of personnel in the use of machine oriented report- 
ing procedures and elimination of duplication of detail. 

h. Full utilization of presently installed data processing equipment 
by insuring that all data processing activities are essentially high- 
benefit programs. 

CONCLUDING OBSERVATIONS 

Last year, 1962, was a year which may very well be recognized as the 
beginning of the first plateau of maturity for the industries that automatic 
control has helped to create. 

Signs of maturity are also evident among the scientists and engineers 
who created these new industries and who must continue to act as whole 
partners with management and finance in continuing to create and exploit the 
new scientific breakthroughs which will firmly establish automatic control 
in its ultimate position as the greatest servant of mankind. This maturity 
takes a number of forms: 

a. The number of engineers and scientists enrolled in post-graduate 
or extension courses in management is testimony that the importance of 
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market and management factors in technical decisions is not widely recognized 
by the technical experts. The idea of technical performance, cost, schedule, 
physical characteristics, and reliability has now become widely accepted as 
a factor in technical decisions. Finally, an increasing number of scientists 
and engineers has come to realize that their technical brilliance is wasted 
if their ideas cannot be sold, and that it is a hollow satisfaction to be 
able to prove that something new and wonderful can be done unless means are 
found to ensure that it will be done. 

b. Electronics engineers, dynamicists, and even civil engineers have 
discovered the benefits of computers; and the electrical, mechanical, and 
field service engineers are recognizing the importance of automation in 
the translation of their diagrams and equations into operating realities. 

c. The overriding importance of reliability in concept, design, manu- 
facture, operation, and maintenance of automatic control systems has been 
thoroughly recognized. 

d. Our new ability to use high-speed, high-capacity digital computers 
as controllers for automatic control systems, plus the advent of micro- 
minaturization, has provided the flexibility for universal application to 
systems of almost any complexity, involving almost any combination of 
scientific disciplines. This new flexibility provides the base from which 
automatic control can be adapted to applications ranging from space vehicles 
to automatic factories — from complex air-traffic-control systems to the 
most microscopic of biological measurements and processes — from the unmanned 
vehicles of oceanological exploitation to the complex man-machine systems 

of industry and sociology. 

15 



e. In the areas of technology, today, because of the broadening appli- 
cation of this industry into all fields, the engineer and the automatic 
control technologist must be able to understand and communicate with technolo 
gists from almost every field of endeavor. This universality of automatic 
control science application has brought about a lowering of the barriers 
of disciplinary specialty to permit an intermingling of the most widely 
diversified technologies. 

Finally, there is no doubt that the electronic computer has had its 
effect and will continue to have an impact, not only on professional 
engineering but on the curriculum and educational program in our schools of 
engineering throughout the United States. 
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INTRODUCTION 



This paper presents a program written around the method of fitting exponental 
curves to data as given by F. A, fillers in his book, "Practical Analysis"^. 
The program fits exponentals with either real or complex conjugate pair fre- 
quencies. The frequencies and the coefficients of the exponentals axe fitted 
using the least mean square error criteria. The number of exponentals to 
be fitted is specified as part of the input data. The program is written 
in Fortran language 



THEORY OF METHOD 



Given N observations equally spaced in the independent variable x by an 
amount h and originating at >: - 0, it is desired to fit a sum of n exponentals 
to them. The observations will be represented by the dependent variable y. 
The desired fit will be written 

y(x) - Co * £ C.eV (1) 
1 * 

where the coefficients, C, are real, ar-d the frequencies, a, are real or 
occur in complex conjugate pairs which give rise to terms of the form 

C e** cos (2) 

in Equation (1). 

The restriction on N, the number of data points, is 

N > n + 1 (3) 

The derivation of the method of fitting the exponentals to the data will 
begin by assuming that the points actually fit the exponental representation 
exactly. An expression involving only the frequencies, a, and the successive 
differences of the data points will be derived. The fiction of an exact fit 
will then be removed by introducing an error term into this relationship e 
The above assumption may be written 

n 

y(m) - Co * I C,e a j (ln - 1)h (t) 
tli 

where y(m) is the m data point. 

Defining the difference between successive data points as 

> y(j*D - y(J) (S) 
Using Equation (It), Equation (5) may be written as 
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Defining 

D(j) 5 e a j h 



f(3) « C.e a j (m - 1)h 



Equation (6) may be written 

n 

- £ (U(j)-l) f(j) 
Using Equations (£) through (8) 

D(m+k) * y(m+k+l) -y(m+k) 



" t f(3) (U(J)-l) U(j) k 



Now n+1 equations of the form (10) will be written out 
D(m) - Ef(j)(U(j)-l) 



D(m+n) - £ f ^J )(ul j ) " 1)u (j) n 



where the sum runs from 1 to n on j, 

Under the assumption that the curves are an exact fit, all of these 
equations (11) hold true. Now if we consider them as n+1 equations 
in the n unknown f(j) then the determinate of the equations must be 
zero. This restriction on the determinate yields the following 
equation 
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n 

\— D (m+k) S(n-k) - (12) 

h 



The functions S(j) are the symmetric functions and are defined in the two 

o 

equivalent forme shown below*". 
S(0) - 1 



s(i) - - J u(j) 



l 



n . n-i » 

(2) ■ S u(3) X u(k) 

1 \2 / 



n-1 

(13) 



or, given an n^ 1 order polynomial whose zeros are the U(j)s, then the 
polynomial can be written as 



n 

i 

o 



S(k)z n " k - (lh) 



where S(0) - 1 

In the problem of fitting the exponentals to the data we do not know the 
value of the symmetric functions before hand. If we did the required 
frequencies could be calculated from them. Mit by removing the fiction of 
an exact fit through introducing an error term into Equation (12) and 
writing it as 



n 

I 





D(m+k) S(n-k) - £(m) 



(1$) 
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the error can be minimi zed by the proper choice of the symmetric functions. 
The means of choosing them will be the least mean square error criteria. 



Squaring Equation (1$) and summing over the permissible values of m yields 

N-(n+l) n N-(n+l) 

) m f) k D(m*k) S(n-k)f - \ n tW (16) 

Since it is the symmetric functions that are being fitted, the partial 
derivatives with respect to the symmetric functions of the left hand side 
of Equation (16) are equated to zero yielding n equations of the form 

n 

S S D(hi+jO D(nn-k)| S(n-k) - (17) 

where JL runs from 1 to n and all sums on m run from 1 to N-(n+l) unless 
otherwise noted. 

The T)(j)s are calculated from the data and the n simultaneous Equations (17) 
are used to solve for the best fit symmetric functions. Equation (lii) is 
then solved yielding for its roots the U(j)s. For a real root the frequency 
is calculated from 

■I n u(j) 

V ~n (18) 
For a complex conjugate pair of roots of the form 

atib (19) 
the corresponding frequencies are calculated from 

a - Tjjj- £n(a 2 +b 2 ) (20) 

P ■ £ arctan (b/a) (21) 
where the corresponding terms in Equation (1) are now of the form (2). 
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Once the frequencies have been found the fitting of -the coefficients of 
Equation (1) to the data is a straight forward problem using the least 
mean square error criteria on the following equation: 

Co ♦ Y Cf'f*- 1 * -y(m) - e(") 

where the coefficients, C, are now to be fitted. The development from 
this point on is Fortran oriented and will be presented in the discussion 
of the program itself. 



METHOD AND PROCEDURE 

The description of the program will emphasize the means used to speed up 
the calculations of the various portions of the problem. Once these are 
understood the listing becomes self-explanatory. 

The first major task is to generate the coefficients of the symmetric 
functions in Equations (17). To accommodate the DO loops the following 
definition will be made 

DE m ( j) = D(m+j-l) - y(m+;5+l) - y(m+j) (23) 

In the program the symbol is written 

DE(j) and the value of m is kept track of by a DO loop. 



Now defining 
m 

J, 

The n equations (17) may be written 



m 

Bj k 5 DET(j)* DE m (k) (2h) 



n+1 i s 



*S(n+l-k) - (2$) 



where j runs from 1 to n 

The desired coefficients are the sums 



E^ k - ^ DE Bl (j)*PE m (k) (26) 

A rapid method of calculating these coefficients will now be described. 

This method makes use of the following properties which are easily derived 

from Equations (23) and (2l). (j) 
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de" 1 ^)- raft. 1+1) 



(27) 



■S " >>i.w (28) 



" < fj (") 

^ " £ <.3 (30) 

The quantities DE, E, and ^ E are stored as they are calculated. The 
procedure is as follows: m is set to one, and all the required values 
of the DEs and the Es are calculated and stored. The values of the Es 
are entered into the coefficient suras of the form (26) as the first term. 
Then m is set to 2 and the DEs are shifted one space to the left using 
the property shown in Equation (27). The missing term is calculated from 
Equation (23). The stored Es are shifted up and left one space using the 
property shown in Equation (28). The missing terms are calculated using 
Equation (2h). The values of the Es are entered into the coefficient sums 
as the second term. The program proceeds as for m ■ 2 until all the terms 
of the coefficient sums have been entered. A separate but concurrent aspect 
is that of storage. The properties of Equations (29) and (30) show that 
only terms on and to the right of the diagonal need be calculated. The 
remaining space is used for storage. The terms E^ ^ are stored in E(j,k+1) 
and the terms ^ E*J k are stored in E(k,j). At the end of the procedure 
described above Equation (30) is used to complete the array of coefficients. 

The n simultaneous linear Equations (25) are solved by a Gaussian reduction 
routine. The solution yields values for the n symmetric functions (13). 

The n^ h order polynomial (Ik) is now solved using a stripped down version 

of the Barstow method „ The roots of Equation (Hi) may be either real 

or occur in complex conjugate pairs. In the program the roots are stored 

in the following manner: Two fields R( j) and M(j) are defined where j runs 

from 1 to n. If the j^ 1 root of Equation (lk) is real then it is stored in 

th 

R(j) and M(j) is set to zero. If the j and j+1 roots are a complex con- 
jugate pair then the real part of the roots is stored in R(j) and M(j) set 
to one and the imaginary part of the roots is stored in R(j+1) and M(j+1) 

is set to two. The field M(j) serves to identify the contents of the field R(j). 
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The program does not use the variable S(j) directly but stores the solutions 
of the n linear equations in the field ::( j) where 

S(J) - x(n+l-j) < j < n (31) 

The entry into the polynomial solver places the symmetric functions in the 
field A(j) where 

A(j) - S(j-l) - x(n+2-j) 1 < J < n, A(l) - S(0) - 1 (32) 
The transfer from x to A is done directly. 

Once the roots have been found the frequencies may be calculated using the 
following form of Equations (18), (20) and (21) 

ctj - ^ J?n(R(j)) wnen M( j) - 

or 

. - * r- arctan | 

If a real root or the real part of a complex root is negative then the 
remainder of the program cannot give meaningful results and the program halts 

The remainder of the program is a straight-forward determination of the 
coefficients of Equation (22), The discussion will be; concerned with a 
rapid way of generating the needed numbers. The most general form for 
n * 3 will be considered here. The generalisation of the results is quite 
straight-forward. Equation (22) is rewritten 

Co ♦ C x e^* 1 - 1 ^ ♦ C 2 e a 2 (ni - 1)h Cos(p 3 (ia-l)h) 

♦ C 3 e a 2 (m " 1)h sin(p 3 (m-l)h) -y(m) - €(m) (3h) 

-8- 



when M(j) ■ 1 
and M(j+1) - 2 



Pefiniiip; DO loop notation 

m 



V(l) • 


1 






in 

V(2) - 








m 

V(3) 2 


e a 2 (m-l)h 


Cos(p^(m- 


•l)h) 


m 

v(L) - 


e a 2 (m-l)h 


Sin(^(m- 


•l)h) 


m 

V(5) 2 


-y(m) 






x(j+l) 


- C. 

3 







(35) 



where m is a superscript or V. 

Using Equations{35) Equation Oh) is written as 

m m m m n ^ * 

x(l)*V(l)+x(2)»V(2)+x(3)W(3)*x(l)^(l4)+V(5) - €(m) W 

Fitting the Cs or the xs by the least mean square error criteria as in the 
first part of this report yields n+1 linear simultaneous equations of the 
form 

n+1 mm mm 

y |^ V(j)*V(k) «at(k)j* ^ V(j)*V(n*2) - (37) 

where j runs from 1 to n+1 and all the sums on m now run from 1 to N unless 
otherwise noted. 

Defining 

m m m ■ f>R\ 

E(j,k) - V(j) * V(k) 

where m is a superscript. 

Equations (37) become 
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£ I J E(j,k)j «00 ♦ ^ E n (j,n + 2) - (39) 
where j runs from 1 to n+1. 

Nov the Vs may be calculated from equations of the form (35) for each m 
but the amount of time spent doing so is prohibitive. Recursion formulas 
may be developed for the Vs by noting that Equations (33) lead directly to 

R(J) " e a j h when M(j) - 

and 

8(J) -e a J h cos ( Pj+1 h) l when M(j) -1 

R(j+ a) . e «^» sin ( p a+ih) | and M(>1) ■ 2 (1,0) 

Now from Equations (35) and (bO) it is obvious that 
m+1 m 

v(i) - v(i) - 1 (la) 

m+1 m , m 
V(2) - V(2) eT - V(2)* R(l) (1*2) 

1 

V(2) - 1 (13) 

The recursion formulas for V(3) and V( U ) are obtained with the help of 
well known trigonometric identities as 

m+1 . 

V(3) - eV" 1 (cos(P 3 (m-l)h) cos - sin (^(ra-ljh) sin ^h) 



1 

V(3) - 1 



m m 

V(3)* R(2) - V(h)*R(3) (WO 



and 



ITH-1 . 

Y(h) - eV" 1 (sin(f5 3 (m-l)h) cos p^h + oos(p 3 (m-l)h)sin p 3 h) 
m m 

- V(M*R(2) + V(3)*R(3) («) 

1 

V(ii) « 

The generalised recursion relationships are 
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m+1 
V(l) - 1 

m+1 m 
V(j) - V(j)^R(j-l) when M(j-l) - 

m+1 m m 
or V(j) - V(j)*R(;)-l) - V(j+l)*R(;j) (i|6) 

m+1 m » IwhaiM(j-l) - 

V (j+1) • V(j+l)*R(,j-l) + V(J)*r(J) |andM(j) - 2 

ra+1 

V(n+2) « -y(m+l) 
and the values for m » 1 are 
V(l) - 1 

V(j) - 1 when M(j-l) - 0, 1 
or (li7) 
V( 5) - when M( j-1) - 2 

V(n+2) * -y(l) 

The program calculates the sums on m in Equations (39) somewhat as before. 
The value of ra is set to one and the initial values of the Vs are calcu- 
lated from Equations (Ul) using M( j) to control the choice of equations. 
The Vs are stored and the Es are calculated using equation (38) and entered 
directly into the sum on m of equations (39) as the first term. The value 
of m is then set to two and Equations (Ij6) used to calculate the new Vs 
with M(j) as the control. The Es are calculated and entered into their sum 
as before. The program proceeds as for m * 2 until the sums are complete. 
As in the early part of this paper only the terras on and to the right of 
the diagonal need be calculated step-by-step. 

The same subprogram is used to calculate the xs of Equations (39) as with 
Equation (25). The coefficients of Equation (1) and form (2) are then 
calculated using 
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Co - x(l) 



C. x(j)l) 2 ♦ x(j+2) 



^j+1 " arctar 



x(j+l) / 




The method of calculating the two variance terms is discussed in 
F. A. Wilier* s book . 
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Gentlemen: 



I would like to talk to you today about a new computer applica- 
tion, Automated Design Engineering. 

A. D. E. is the use of the computer in the design of non- 
prototype products, the type of engineering commonly referred to as 
custom engineering, application engineering, or product engineering. 
Specifically, A. D. E. can benefit your company through significantly 
reducing your engineering lead time, increasing your engineering 
productivity , and decreasing your engineering costs. These and the 
other advantages of the application can, in turn, result in an improved 
competitive position for your company in your industry. 

First, I would like to tell you what A. D. E. is. Second, I 
will discuss where A. D. E. applies. Third, I will spend a little time 
on how A. D. E. works. Fourth, I would like to cover in more detail 
the advantages of A. D. E. 

Computers have been used, in the past, in the design of many 
products for industry. Computers have been utilized successfully in 
the design of circuits, missiles, motors, transformers, and telephone 
equipment, to name just a few. In each case, the use of the computer 
has resulted in great s avings and great increases in efficiency and 
productivity for the companies involved. Now a newly developed 
computer application, Automated Design Engineering, vastly increases 
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the power of the computer as an engineering tool . With A. D. E. , the 
computer can now be used to solve more engineering problems on a 
wider range of producta 

Conceptually, A. D. E. involves the storing of design logic 
in a computer so that the computer can accept customer orders as 
input and automatically generate complete designs. The input to the 
system would be the same type of orders you now receive from your 
customers. The completed design would contain the same type of 
information normally given by your engineers to your manufacturing 
department; such things as product characteristics, part numbers, 
assembly numbers, bills of material, purchases parts list, drawing 
numbers and so forth. 

The principles of A. D. E. can best be described by using an 
example. An A. D. E. system for a company that manufactures 
pumps, for instance, would receive customer requirements i n the 
form of orders or requests for bids . In this case, the customer needs 
a pump to pump 150° Fahrenheit carbonic acid at the rate of 100 gal- 
lons per minute with a head of 40 feet and he wants the pump motor 
wired for 220 and 440 volts. These requirements would be entered into 
the computer and the completed design would be printed out. 

These requirements would be entered into the computer, would 
be processed by the computer with the aid of the stored design logic, 
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and the completed design would be produced. The completed design 
information can be in many different forms . . . such as part numbers, 
assembly numbers, drawing numbers, manufacturing instructions, 
etc. In this example the completed design information consists of 
the pump frame number, FQ 6; the model number CXR; the suction 
and discharge pipe diameter, 3 and 4 inches, respectively; the im- 
peller diameter, 8 3/4 inches; the motor speed, 1750 RPM; etc. , 
etc. In short, the completed design would consist of the information 
needed by manufacturing to build the product. 

Now that we have discussed what A. D. E. is, where does it 
apply ? Automated Design Engineering applies primarily to "cus- 
tom" or "product" design which can be defined as "non-prototype 
custom design variations of a standard product line to meet your 
customers' requirements on a continuing basis. " Even though some 
of the tools, techniques, and methods available through A. D. E. are 
applicable to prototype design, the major impact of A. D. E. will be 
in the custom or product design area. 

There are many companies in industry today which produce 
custom designed variations of a standard product line in response to 
customer orders or requests for bids. "Custom" or "product" de- 
sign is common in such products as pumps, motors, generators, 



switch gear, transformers, electrical measuring equipment, indus- 
trial furnaces, switchboards, engineering and research instruments, 
heat exchangers, steam turbines, conveyors, and air compress ers, 
to name just a few. 

Let us look at the problem of the "custom" or "product" 
engineer. In industry today, customer orders for products will 
come to the design engineer. His problem is three-fold: 

First, he must determine what to build to satisfy customer 
requirements. 

Second, he must translate these customer requirements 
into workable parts and assemblies and their drawings. 

Third, he must prepare the complete paperwork for manu- 
facturing. It has been determined through actual experience in 
industry that these problems can be solved on a computer with 
great savings of time and money and great increases of efficiency 
and productivity. 

Now to help you visualize what a system would look like in 
your company, let us look at a typical Automated Design Engineering 
System in operation. A customer order coming in would go tb the 
Engineering Department where two functions would be performed. 
First, the order would be edited to insure the completeness and 
validity of the customer order. Second, any necessary pre-engineer- 
ing would be performed. Once the order editing and pre -engineering 
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has been completed, the customer order would go to the key punch- 
ing room where the customer requirements would be transcribed 
onto cards and fed into the computer. 

The computer is programmed with design logic and has 
available to it tables, standards, and reference information. The 
output of the system would be completed design which would then 
be reviewed by the engineering department and passed on to manufac- 
turing. Now that we have covered what A. D. E. is, and where it 
applies, and what an operating system would look like, how does an 
Automated Design Engineering Syste m actually work ? The key to 
Automated Design Engineering is the ability to capture the design 
logic by which customer requirements are translated into product 
specifications. 

We have developed new tools, new techniques and new 
methodologies which will assist you in developing an A. D. E. system 
for your company. One of the most important tools, and the heart of 
this new system, is Decision Tables , a technique for capturing the 
design logic that comprises custom or product design. The ability 
to capture design logic is the key to A. D. E. Before we discuss 
Decision Tables, let us clarify what we mean by design logic. 
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What is design logic. . . ? Design logic is the complex 
decision-making process through which the engineer proceeds in 
designing a product to meet a particular customer's requirements. 
The engineer reads the customer order and then, through the medium 
of design logic, proceeds to design and select the various parts and 
assemblies required to satisfy the customer's requirements. The 
result is the completed design of the product. 

Let us look at a sample of design logic and how this design 
logic can be captured using decision tables. We will use as our 
example the armature for a voltmeter. If we were to ask an 
engineer working for a company that manufactures voltmeters how 
he designs the armature for a voltmeter, he might say: "Well, if 
thercustomer requires DC service for a speed application and asks 
for a single phase instrument with villivolt rating, then I know I have 
to use a moving coil. lf_ two windings are needed, which frequently 
is the case, then the part number will be -12526A. Drawing number 
A26A will be used in this case. If the rating value specified is be- 
tween 76 and 200 millivolts and we need a moving coil, the main 
winding will use Aluminum wire 16 mils in diameter. The number of 
turns and the number of layers will come from these two formulas 
which relate to each other. Based on the case we're using, we'll 
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have 13 turns on the main winding and one layer. Again, I have to 
look up the drawing number . . . it's 012526- 1A. Now the coil will 
also need a damper winding which for a 76 to 200 millivolt rating and 
a scale size specified as four inches, will take 8 mil Copper wire and 
half the number of turns on the main winding which is also shown on 
drawing number 012 526-1 A. " This complex decision-making pro- 
cess can be captured by use of a new tool called decision tables. 

A decision table consists of four quadrants. In the top two 
quadrants are the customer requirement names and values. In 
the lower two quadrants are product specifications, names and values. 
This decision table is a summary of part of the data given to us 
by the voltmeter design engineer. In the upper left hand quadrant 
is a list of customer requirement names. . . service, application, rating 
units and number of phases. Values which the customer might specify 
for these requirements are given in the upper right hand quadrant. . . 
DC, speed, millivolt, 1, etc. In the lower right hand quadrant are 
values for these product specifications moving coil, inductive, 

1 plus the number of phases, etc. The decision table is read as 
follows: If the service is DC, and if the application is temperature, 
then the type of armature needed is amoving coil, and tables number 

2 is the table to which we should go next. Rule 4 would be read as follows. 
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If the service is AC, and if the rating units are milliamps, and if the 
number of phases is 1, then the type of armature required is inductive, 
and the number of windings is calculated as being ONE PLUS the number 
of phases, and the next table is number 2. Each different set of conditions 
which can exist is called a rule . In this case there are six feasible 
sets of conditions or six rules. Once the decision table is established, 
programmed and entered into the computer, the computer will auto- 
matically select which rule applies to each individual customer order. 

It takes hundreds of decision tables to capture and store the 
design logic for an entire product or product line. When all of the 
required decision talbes have been programmed and placed in the 
computer memory, customer orders can be entered into the computer 
design in a logical stepwise manner. 

Now that we have covered what A. D. E. is, where it applies, 
and have gone into a little detail on how it works, let us turn our attention 
to what an A. D. E. system can do for you, or the advantages of an 
A. D. E. system. We have already indicated that you can reduce 
your engineering costs and at the same time, obtain increased productivity 
and efficiency from your engineering force. Overall lead times can, 
in addition, be significantly reduced through a reduction of design 
engineering time, materials procurement time and manufacturing lead 
time. The design time, for example, can be reduced from weeks or 
even months to a matter of minutes. This reduction in design time 
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will allow you to give faster response to customer orders and requests 
for bids. Faster response to bids, on the one hand, can have a significant 
effect on your profits if there is a high correlation, as exists in many 
industries, between speed of your response and acceptance of your 
bids. Faster response to orders, on the other hand, can result in 
faster deliveries and improved customer relations. And furthermore, 
Automated Design Engineering can provide a significant increase in 
your business activity by allowing you to respond to more bids if you 
now find yourself unable to respond to all of the requests for bids that 
are made to your company -- due to the lack of available time in your 
engineering department. In addition, faster response to bids can have 
a significant effect on your profits if there is a high correlation, as exists 
in many industries, between speed of your response and acceptance 
of your bids. Also, bid costs can be reduced because A. D. E. can 
reduce the out-of-pocket cost incurred when bids are designed but not 
won, since the engineering cost per bid can be greatly reduced. 

Material savings are another potential advantage of Automated 
Design Engineering. By calculating exactly how much material is required 
for each job, A. D. E. can reduce overdesign and waste. Another way 
in which material savings can be realized is through reducing the number 
of parts in inventory with the same specifications but different part 
numbers. As you develop your A. D. E. system, these duplications 
will become readily apparent and can be eliminated. The consistent 
use of the best design practices eliminates the proliferation of methods 

a 
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and materials and insures a design commensurate with requirements. 
All too often today, the pressures of competition and of constantly 
shrinking delivery times tend to result in over -design or in picking 
a design which is perhaps more expensive than the specifications 
actually call for and the result is loss of profit. Bid and order cost- 
ing can be incorporated into your Automated Design Engineering 
system to allow you to not only design but also to price both labor 
and material for your design. Through the automatic generation of 
engineering paperwork much of the clerical burden which wastes so 
much of the time and talent of engineers today can be alleviated. 

Automated Design Engineering also provides an improved 
error-checking facility. The sources of errors in engineering today 
are many. These errors can creep in through the customer order, 
through the sales engineers, through engineering errors or material 
list preparation. Errors can result in the wrong material, too much 
material, or too little material being available at assembly points. 
They result in wrong design and in over-design. With an Automated 
Design Engineering System these functions are performed automatical- 
ly. Error checking capabilities are built into the computer to catch 
and eliminate such miscalculations, which are the real problem in 
industry today. Greater management control is possible because 
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management policies can be incorporated in an Automated Design 
Engineering System with the design logic, to insure that management 
policies on engineering are carried forward. 

Now that we have discussed the advantages of A. D. E. , I 
would like to show you how we can assist you in developing and im- 
plementing an Automated Design Engineering system in your company. 
IBM has available tools, techniques , and methods to help you in 
developing your system. There are two steps in the development of 
an operating Automated Design Engineering system. These are the 
Survey and the Implementation . The survey represents the first 
step. Its purpose is to determine the applicability of the A. D. E. 
system to your product lines. It will provide you with an analysis 
of your present system. It will determine the requirements of the 
new system and will give you a preliminary design of the new system 
as it can be applied specifically to your company. Lastly, upon com- 
pleting the survey, you will be able to measure the advantages that 
will accrue to your company through the use of Automated Design 
Engineering. The new tools have been developed to prepare an 
analysis of your present design engineering operation from four re- 
lated but different standpoints : Those of time, cost, accuracy, and 
o perations . After the survey is completed you are ready for the 
next step, implementation . The purpose of the implementation phase 
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is to develop a completely tested and operational Automated Design 
Engineering system. In order to do this, an analysis of your cus- 
tomer specifications, and an analysis of your product structure 
must be prepared. The techniques and tools which IBM will provide 
will also enab le you to capture the design logic of your product line 
onto decision tables, to perform the detailed systems design, the 
programming, the testing, the conversion and finally the initial 
operation of your Automated Design Engineering system. These tools, 
techniques and methods are tested and proven. They are available to you 
in the form of printed material. For example, the A. D. E. General 
Information Manual will introduce you to the survey and implementa- 
tion phases of this new system. In addition, we have prepared a de- 
tailed reference manual to provide your engineers with the "How To 
Do It" information necessary to develop an A. D. E. system. 
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I. INTRODUCTION 

It is no secret that the 1620 is not designed for excessive 
data handling and the IBM Corporation definitely is not marketing it 
for commercial uses. In fact it is a small scientific computer and 
as its many users will attest, it performs with excellence in this role. 
Despite the machine's fine record in the field, its inability to handle 
excessive data (especially the basic 1620 with paper tape) produces 
one distinct problem area which can present itself during the selling 
phase or after the installation of the system. It is known as machine 
utilization. 

Present day engineering or scientific establishments do have 
considerable, if - not abundant, amounts of commerical (accounting) type 
of work which could be automated and thus increase the utilization of a 
computer. This problem of the 1620 utilization is not as predominant 
in a large data processing oriented company as it most probably has a 
commercial computer system. In such a company the 1620 is being looked 
upon for the engineering or research and development departments. A 
large company would presumably not expect high utilization on a machine 
which would be operated in an open shop atmosphere. Their objectives 
would be more intangible, such as the release of engineers for more 
creative type work and the increase in speed of routine computations. 

On the other hand, in a small engineering account, where the 
rental cost of $1,600 per month represents a large investment, the subject 
of machine utilization is of great concern. Even though a two or three 
hour run a day can justify the installation of a 1620 in such an account 
there still is the feeling that more should be gotten for this kind of 
money. It is for such an establishment that this paper is primarily 
intended. 

The company for which these programs were developed has a 
bi-weekly payroll for two offices, one of ninety and the other of forty 
employees. It takes about two days (16 hours) of the computer's time 
every other week to produce the complete payrolls and labor cost distri** 
but ion. The operation of these programs since January 1, 1963 has helped 
to increase the average utilization of the computer to 40 hours per week. 
As for economic value of these programs, it cannot be measured in dollars 
and cents. It places valuable up-to-date information about the cost, 
progress and estimation of projects in the fingertips of management in a 
matter of a few hours. Neat reports may be produced in a moment's notice. 
It may be said that these programs in addition to increasing machine 
utilization do contribute in improving management operation. 
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II. SYSTEMS ANALYSIS AND DESIGN 

In trying to analyze and successfully design a payroll and 
labor cost distribution program package the following must be taken 
into account: the memory capacity of the computer, the presence or 
absence of special features, the type of input-output and the method 
of operation of the shop including personnel availability. The concern 
where these Payroll and Labor Cost Distribution Program Packages were 
developed employed the basic 1620 with a 20k memory, without any special 
features, and with paper tape input-output hardware. This concern, a 
consulting engineering outfit, operates on a semi-open shop basis without 
a specialized machine operator. One cannot start any simpler than that 
nor can one have any more handicaps as far as the 1620 is concerned. 

Of the three major languages of the 1620, machine, SPS, and 
Fortran, the latter was chosen by the authors because the programs could 
very easily be maintained and revised according to the needs of the 
company. In addition, Fortran offers a lot easier means of programming 
and debugging and to the best of the author's knowledge such a task would 
be a first. It must be noted here that all of the following programs 
(see exceptions later in this paper) were originally written in Fortran 
with Format and later changed to UTO Fortran. For those unfamiliar with 
UTO Fortran, it was developed by Mr. E. Stewart Lee and Mr. James A. Field 
of the University of Toronto, Ontario. UTO in general is similar to Fortran 
with Format except that it brings the program origin down to 06950, is 
slightly faster, the subroutines are about 10 per cent shorter, the use 
of EXECUTE PROCEDURE n statements (similar to CALL statements) facilitate 
programming and save memory when properly used, and has quite a flexible 
input format (input data does not have to be in strict accordance with the 
format statement). 

One question that had to be answered early in the planning stages 
was the availability of checks for the 1620 typewriter. Cardboard checks 
were out because they necessitated a typewriter RPQ. Paper checks were 
chosen and designed by IBM. A picture of a check appears in the appendix 
of this paper. 

The check is of the high-low design with the stub under the 
check proper. The perforation of the right can be placed any where. 
When designing a pay check the following hints should be considered: 

(1) When designing the check, it is advisable to make the 
spacing between the printed lines a multiple of two or three. This means 
that at program time during the carriage control operations, the type- 
writer can be run in double or triple space mode and thereby save a 
considerable number of carriage returns. 
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(2) Hint (1) may be accomplished with statements of the 
form 1000 FORMAT (//) or if using UTO Fortran with the CONTROL 102 
statement. In UTO the statement CONTROL }_08 may be used for tabulation. 

(3) The check should be designed so that variable alphameric 
data appears only on the left edge. Note the example check where the 
two dates (XX/XX/XX) and the name are all left adjusted. This is very 
important in a situation such as the following. The employees name is 
part of the master record and must be the first data in that record. 
For example, the name and year to date gross, FICA, and withholding 
taxes would read in as: 

ACCEPT TAPE 100, YGR, YFICA, YFWT 
100 FORMAT (15H EMPLOYEE NAME F9.2,F7.2,F8.2) 

The variable alphameric name is stored in the format statement 
itself. When saying PRINT 100, the name alone will be printed out. If 
the name were not left-adjusted in the master record and in the Format 
statement, it would be impossible to get it out alone. Another technic 
that could be used is to carry the name all by itself in a separate 
master record. 

The payroll system for which these program packages were 
designed is based on the following criteria: 

(a) Employees are to be paid bi-weekly with time being 
reported on two time sheets (the programs may easily be changed to 
any other pay system). 

(b) Reported time is to be reported by each employee by 
operation (also referred as activity herein) for each project. 

(c) There are three types of paying systems (pay codes): 
straight salary (code 1), straight time (code 2) and time and a half 
for overtime (code 3). 

(d) There may not be more than 130 employees. This company 
uses 20 activities only. If need be this number may be increased at 
the expense of running time • 

(e) Three master records are to be kept: (l)the payroll master 
tape containing the name, number, rate, pay code, dependents, year-to-date 
gross, FICA, tax withholdings and expenses, quarter-to-date gross, FICA, 
tax withholdings, social security number, hospitalization, bonds, savings, 
life insurance, pension, and miscellaneous deductions; (2) the project 
master tape containing the project number and todate totals on direct labor, 
expenses, and indirect labor (sick leave and vacation); (3) the activity 
Master tape containing to-date totals for each activity for each project. 
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As mentioned in the introduction the biggest handicap in 
this operation is that of handling data. Because of this, changes 
to any of the master tapes should be kept to a minimum as the slightest 
revision requires the reproduction of the entire tape. Master tape edit 
programs may be written to facilitate tape revisions, especially massive 
ones such as the zeroing of the registers at the end of a quarter or the 
year. 

The other case where data handling enters the picture is in 
the data. This is discussed in the following section. 



III. THE PROGRAMS 

Below is a general description of the programs in the Payroll 
and Labor Cost Distribution Programs. 

A. Data Preparation 

Considerable thought was given to the input data since this 
was the start of a chain of events which would produce all of the pay- 
roll and labor distribution reports. As part of this paper's appendix 
one will find the time sheet used for reporting work time and expenses. 
Since this data form the basis to both packages it is typed into the 
Data Preparation program* The sequence of this data is composed of a 
minimum of three records per employee. The first contains the employee's 
number and the remaining contain the activity, number of hours or 
personal expenses and the number of the project. A negative activity 
denotes the end of a time sheet. For a bi-weekly payroll two negative 
activities are required. The program immediately punches on tape the 
accepted data in a similar format. This is the Labor Distribution 
Data Tape. Then it proceeds to compute the regular and overtime hours 
and personal expenses on a weekly basis which is stored in memory for 
each employee. This information will later produce the Payroll Data 
Tape. To do this the Payroll Master Tape containing the employees' 
pay codes is read into the program at the start of the run storing each 
employee's number and pay code. It should be mentioned here that 
activities 18, 19 and 20 are reserved by the programs for vacation, 
sick leave and personal expenses. 

One bit of logic that might be helpful in the program is a 
technique for being able to go into memory and make corrections after 
all data is typed in and while the payroll data tape is still stored. 
This could be accomplished by assigning a storage cell number to each 
employee as he is entered into memory for the first time. This storage 
cell number would then represent the position in all four arrays (man 
number, regular hours, overtime hours, expenses) where this particular 
employee's data could be found. The number, along with a memory map of 
the arrays would then give the exact location of any data in question. 
A portion of such a table is included in the appendix. If an error were 
made in the typing of data for the say ninth employee the table will show 
that this employee's number, regular hours, overtime hours and personal 
expenses would be found in addresses 19750, 18450, 17150, and 15850 
respectively (see table in appendix). 
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It will be advantageous in the Data Preparation Program to 
have a dump of grand totals at the end showing regular hours, overtime 
hours and personal expenses. These totals should balance back to an 
adding machine tape which should have been taken when auditing the 
time and expense sheets. The totaling routine should be programmed so 
that it may be "branched to M after any corrections to the memory table 
have been made. At this point if everything balances, the Labor 
Distribution Data Tape may be removed from the punch and the punching 
routine of the Payroll Data Tape initialized. 

This marks the end of the combined operation between the 
Payroll and Labor Cost Distribution Programs. 

B. Payroll Register 

The first step in the Payroll Program Package is to obtain 
the payroll register showing all incomes and deductions for all employees 
of the current period. This, the Payroll Register Program . necessitates 
as input the Payroll Master Tape and the Payroll Data Tape as well as the 
starting check number and period ending date which are entered via the 
typewriter. The payroll data is stored in four arrays - employee number, 
regular hours, overtime hours and expenses. The Payroll Master Tape is 
then read in, the arrays are searched for a match, and the pertinent data 
from the master and detail are combined to calculate the net pay for each 
individual. Two lines of printed output appear for each employee showing 
the standard payroll register items. Accumulative totals are kept for 
each item and at the end of the report these totals are crossfooted for a 
check against the accumulated net pay. 

In this program as well as in the next there is the problem of 
half adjusting and truncating a product. Let it be assumed that a straight 
time employee waged at $4,375 per hour has worked 48.50 hours. Hence 

PAY - 4.375 x 48.50 * 212.18750 

For obvious reasons the employee should be paid $212.19. Fortran 
however will carry it as 2121875003. When several such products are 
accumulated there is going to be a disagreement between the actual 
sums and those of Fortran. This may be remedied by the following 
statement : 

PAY ■ RATE x HOURS + 0.005 + 100000.0 - 100000.0 

which will be executed as follows step by step 

PAY • 4.3750000 x 48.500000 - 212.18750 
PAY « 212.18750 + 0.0050000 « 212.19250 
PAY « 212.19250 + 100000.00 « 100212.19 
PAY « 100212.19 - 100000.00 « 212.19000 

When completed the payroll register should be spot ehecked by 
th<* ^personnel responsible to insure complete accuracy. There is no punched 
output to this program. 
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C. Payroll Checks 

As soon as the payroll register is checked and accepted the 
Payroll Checks Program may be used to produce the paychecks. Basically 
this program is similar to the previous one using the same input (instead 
of the check number the check date is entered). In addition to the checks 
this program updates the Payroll Master Tape concurrently with the checks 
and at the end it produces an Adjusted Pay Rates tape. Such a rate is merely 
the employee's gross pay divided by the number of hours worked. This will 
affect only the straight salary and time and a half employees. For example, 
a straight salary employee who is always paid for 40 hours and has a rate 
of $4,375 has worked 48.50 hours. To give a more complete picture in the 
distribution reports all of the 48.50 hours should be considered. Thus 
this employee's adjusted pay rate is 

Adjusted Pay Rate - 4.375 x 40.00/48.50 - 3.6082474 

Had the employee being paid time and a half for overtime then his adjusted 
pay rate would be 

Adjusted Pay Rate - (4.375 x 40,00 + (4.375 x 1.5) x 8.5)/48.50 - 4.7996134 

This adjusted rate may later be used to distribute dollar increments of the 
48.50 hours over the various projects on which the employee worked. This 
rate can be calculated at the same time as the employee's check and stored 
in memory. Since the Payroll Master tape is updated during the production 
of the checks the Adjusted Pay Rates will have to wait until the end of the 
program to be dumped on tape. These rates are later used with the labor 
cost distribution programs. 

D « Payroll Deductions Register 

Because of memory and output limitations not all deductions are 
itemized in the payroll register. These deductions, bonds, savings, 
hospitalization, etc. are itemized for each employee by the Payroll 
Deductions Register Program utilizing the just updated Payroll Master Tape. 

E » Labor Distribution Data Sorting 

Having completed the payroll phase of the operation we may 
proceed with the labor cost distribution phase. 

It was stated under the Data Preparation Program that data 
read directly into the program from the time sheets were divided into 
the Labor Distribution Data Tape and the Payroll Data Tape. Since it is 
possible for an employee to report time under several activities for several 
projects time data will have to be sorted in one way or another. Such a sort 
would be easy for a card machine. For a tape machine though it is a 
different story. Data could be sorted manually before the Data Preparation 
Program run but this would defeat payroll automation. As a matter of fact it 
would be a rather involved affair. Also mentioned earlier was the Project 
Master Tape containing all projects of the company with their respective costs. 
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It is possible for the operator to type an erroneous project number 

that is not included in the Project Master Tape or that a new project 

had been added without the roaster being updated. For these reasons two 

programs were developed in machine language to utilize memory space. 

These are the Pro ject Completeness Test and The Labor Distribution |^ 

Data Sorting programs # *^~ y 

The Project Completeness Test program first reads the Project 
Master Tape and stores the project numbers only; then it reads the Labor 
Distribution Data Tape and compares each project number with the table 
previously stored in memory. If a match does not occur the project 
number on the Labor Distribution Data Tape is printed on the typewriter. 
At the end of the run the number of projects and number of records are 
typed. The number of records refers to the records of the Labor Distri- 
bution Data Tape containing an activity, hours or expenses and a project 
number. When this test has been made the Labor Distribution Data Tape 
may be corrected if need be through a tape correction utility program. 

The Labor Distribution Data Sorting Program does a little more 
than just sort data. First it counts records to insure that no un- 
authorized records have been dropped or added while correcting the Labor 
Distribution Data Tape. Second it sorts data. The program inserts the 
employee number to each of his records, changes each record from an 
alphameric mode to numeric mode dropping all decimals and arranging them 
in the format shown below: 

XXXXX XXXKXXXXXX + 
1 2 3 4 

where 1: project number 
2: employee number 
3: activity number 
4: hours or personal expenses 

Each record contains a total of 16 characters including the record mark. 
These records are stored in memory beginning in location 19999 and working 
on downwards. The sorting of the records begins when all records have been 
stored. It is a simple replaceable sort routine which starts at the top of 
memory and keeps comparing two adjacent storage cells (records). If the 
upper is less than or equal to the lower, no replace takes place. If the 
upper is greater than the lower, the two records are interchanged, a switch 
is set and the program steps down one record to compare the next two. This 
continues until one complete pass on the records has been made. At this 
point, the switch is checked. If it is on (meaning at least one record 
out of sequence), it is turned off and the comparing routine is repeated. 
When the switch is finally off after a complete pass, the internal sort 
is considered terminated. It should be emphasized here that during the 
compression phase (changing data from alphameric to numeric) where data 
is edited in the input area and placed in a cell, the comparing field can 
be set up as one field and still get three breakdown levels - activity, 
within man, within project - providing the field is arranged as described 
in the previous paragraph. 
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Third, at the completion of the sorting routine the program 
will commence reading the adjusted pay rates which were produced at 
the end of the payroll check production. The rates and employee numbers 
are converted into numeric characters and the memory cells are searched 
for an employee number match. When the match occurs the rate is multiplied 
by the corresponding hours (no multiplication is performed for activity 
20 - personal expenses) and the hours are now replaced by the cost. Each 
employee number of the Adjusted Pay Rates tape searches the entire memory 
map as it is possible for one employee to appear in several records. 

Fourth, when all hours have been convered into costs the 
program will type the total direct and indirect labor cost which should 
agree with that of the payroll register. It is obvious that the adjusted 
pay rates of personnel other than those of straight time will not be the 
same with their original base pay. These adjusted pay rates may contain 
as many as eight significant digits. This problem is similar to that of 
the payroll register and payroll checks of rounding or truncating products. 
To compensate for this each product (rate times hours) is half adjusted to 
the nearest penny and all truncated parts are accumulated to form the 
truncation or rounding error. This which may be negative is assigned to 
project one, administration, to activity one, general, and to employee 
zero, fictitious. Hence the gross payroll figure plus the rounding error 
of this program should equal that of the payroll register. 

Fifth and last the fields of each record are changed back to 
alphameric and are punched on tape to produce the Sorted Labor Distri- 
bution Data Tape. 

F. Labor Cost Distribution 

This program is essentially a listing with some accumulating 
of the sorted labor data tape. This tape is now in project sequence. 
Within each project are the records of each man who worked on that 
project during the last pay period. Within each man are the operations 
which he performed on that project as well as the distributed money. 
There are actually three phases to this program which give three separate 
reports. 

The first report shows by employee number all the employees 
who worked on each project this current period. All the activities 
for each employee are accumulated and this accumulated money total is 
segregated into direct and indirect labor and personal expenses. 

The second report of this program is merely a listing of the 
totals lines for each project. It gives a neat, condensed, summary 
report of what the projects did this period. This is accomplished by 
storing project number, total direct labor, total expenses and total 
indirect labor for each project while producing the first report. 
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The third report is the project status report showing all 
projects of the company whether work was done or not on all of them. 
For this the Project Master Tape is used. As the master is read in the ,*0 
project number from the tape searches the memory (see previous paragraph) \J 
for a match. If a match does occur the figures are combined and the 
master is updated. If no match occurs (no work done on this project 
this week) the figures remain the same. There may be printout options 
at this point. One which the user would probably want is a project-to- 
date report showing three totals for each project. Line 1 shows project 
totals to-last period, line 2 shows project totals of this period (zeroes 
may be present) and line 3 shows the updated project totals (sum of lines 
1 and 2). 

A point of interest to this program as well as to the last of 
this series is that of the grand totals appearing at the end of the 
various reports. Fortran handles only eight significan figures thus 
limiting totals to $999,999.99. If the elements of a total added to 
more than a million dollars the machine computed total would be a 
truncated figure sometimes several dollars off the actual sum. This 
problem may be eliminated using Fortran in the following manner. Let 
SUM represent the accumulative value of an item and COST the value of 
any one element of the same item so that 

SUM - Z C0ST 

Before incrementing SUM by COST, COST is compared with the difference 
between 999,999.99 and SUM. If the cost is smaller than the difference 
the addition is performed. If not a carry-over factor is incremented 
by one and the sum reinitialized by the difference between COST and the 
first difference. Thus in UTO Fortran a procedure could be set up as 
follows : 

BEGIN PROCEDURE 100 
DIF - 999999.99 - COST 
IF(DIF-COST) 111, 110, 110 

110 SUM - SUM + COST 
RETURN 100 

111 KARRY - KARRY + 1 

SUM - COST - DIF - 0.01 
END PROCEDURE 100 

where KARRY indicates the millionth carry-over factor. Note the penny 
(0.01) in the last statement of the procedure to take care of the 
difference between a million and 999999.99. 
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G. Labor Cost Breakdown - Bi weekly 

This program merely rearranges data from the Sorted Labor 
Distribution Data Tape. That is, it reads in all the details for a 
project number, accumulates the amounts by activity number, and punches 
out a new tape, the Labor Cost Breakdown Tape, which contains the project 
number followed by each of the operation numbers and their respective 
monies. This is done for the entire input tape, so that the output tape 
contains all the project numbers followed by all the activity numbers 
and amounts. This program generates a new activity code number for each 
project Activity 21 representing an accumulative total of all the money 
from all the other operations on that project. 

There is also an optional print-out to this program. Any or 
all projects may be printed showing money spent by each employee under 
each activity for a project. Because of space limitations and the extend 
of the arrays only eight employees are to be included on a page. Thus 
several pages may be needed for one project. 

H. Labor Cost Breakdown - Accumulative 

This is the last of the Labor Cost Distribution series of 
programs. As one may recall the third report of the Labor Cost Distri- 
bution Program updates a Project Master Tape showing labor breakdown 
by direct labor, personal expenses and indirect labor. This last 
program gives the labor breakdown by operation within project. Because 
there are 21 possible operation under each project, the Activity Master 
Tape is actually composed of three reels of tape. The first carries all 
the projects with activities 1 through 7. The second aarries all projects 
with activities 8 through 14 and the third carries all projects with 
activities 15 through 21. As a reminder activity 21 represents the 
total cost of activities 1 through 20 inclusive. 

This program is executed in three passes. At the beginning 
of each pass the Labor Cost Breakdown Tape is loaded and the program 
selects and stores activities 1 - 7, 8 - 14 and 15 - 21 depending whether 
it is pass one, two or three. Having stored this information the 
corresponding Activity Master Tape is passed printing the breakdown 
report and punching the updated Activity Master Tape. The report may 
take a similar form to '-that of the project status mentioned under the 
Labor Cost Distribution Program. The company that developed these 
programs however chose to eliminate from this report the top line 
showing project-to-last period totals. 

This final report is an extremely comprehensive breakdown 
of the allocation of funds on the various projects. When a project 
is completed, the final master records become an invaluable reference 
to be used in the future for estimating costs of similar jobs and in 
a way to predict or forecast job progress and man power alocation. 
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TV. TIMING CONSIDERATIONS AND LIMITATIONS 

Below is a list of limitations to these programs: 

(a) The number of employees is limited to 130 in any one 
run. Slight variations in paying systems and in state or local tax- 
ations may increase or decrease available memory and thus introduce 
new limitations on the number of employees. 

(b) In the Labor Distribution Data Sorting Program there 
is room for 1048 records. 

(c) An increase in the number of activities will increase 
the number of passes in the Labor Cost Breakdown Accumulative Program. 

(d) The following timings are based on 90 employees working 
on 35 projects at any one time and reporting their time in about 500 
records. The total number of projects being reported on the distribution 
and breakdown reports are 70. It should be kept in mind that these 
figures of time are for a bi-weekly payroll. 

1. Time and expense sheets are turned in by Friday night. 

2. Time and expense sheets are audited Monday morning. 

3. Revisions and adjustments to the master 



tapes 1.00 hours 

4. Data Preparation Program 2.50 " 

5. Payroll Register Program 0.75 " 

6. Payroll Checks Program 2.25 M 

7. Payroll Deductions Register Program 0.50 " 

8. Labor Distribution Data Sorting ................ 0.75 M 

9. Labor Cost Distribution Program 

Report #1 ...1.00 M 

Report #2 0.25 " 

Report #3 (project status) 0.75 " 

10. Labor Cost Breakdown - Bi weekly 

without any print-out 0.50 ,f 

(allow 5 minutes per project 
print -outs) 

11. Labor Cost Breakdown - Accumulative 2.25 " 

12. Average minimum total •••• 12.50 M 
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The secret to a smooth and fast run is to avoid last minute 
changes, and to train the employees to report their time and expenses 
properly thus eliminating corrections and revisions while a $1600 a 
month machine is running doing nothing but tape corrections. 



V. SUMMARY 

It was the purpose of this paper to present guide lines and 
documentation for implementing a payroll and labor distribution application 
on a basic paper tape 1620. It takes but a little imagination to develop 
additional utility programs to handle information developed by these 
programs to prepare special reports such as Employee Taxable Wage 
Reports, W-2 Form Information (to the best of the author's knowledge 
there are no W-2 forms that will fit the basic 1620 typewriter), Sick 
Leave and Vacation Tally Reports and a number of others. With slight 
if any at all modification the distribution programs may be used to 
handle other distributions unrelated to the payroll. 

These programs have supplied their developers with up-to- 
the-minute information on the financial status of all their projects - 
information which in the past was delinquent and incomplete. It is 
the authors* belief that such programs as these presented in this 
paper although not money making will provide a user with invaluable 
service • 



VI. APPENDIX 

Included in the appendix to this paper is a general flow 
chart of both program packages and input-output samples of their most 
important phases. 

The back side of the time sheet is not used in either of 
the packages. In the register of payroll deductions the last column 
represents the horizontal sums of the other columns. Note that the 
grand total is obtained twide; one from the vertical and once from the 
horizontal sums. 

Two projects of the labor distribution report are shown on 
thi» same page. Actually they would appear each on a separate sheet 
the carriage being controlled by the program. Totals of the project 
summary are obtained at the end but are not shown in the sample. The 
project status at the bottom shows a sub-total and a total. The 
developer of the programs wanted to divide the projects into two groups 
each having a subtotal. The subtotal shown is that of the second group 
of projects. 
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The bi-weekly Labor Cost Breakdown report shows a partial 
listing of project No. 61151, There are only 16 of the 20 activities 
shown. The number across the top indicate the employee numbers. Had 
there been more than eight employees a second or even more sheets would 
have been needed. In such a case the column headed ACC would indicate 
the accumulative totals to that sheet. The last display shows a partial 
Accumulative Labor Cost Breakdown report. The numbers 15 - 21 across 
the top indicate the activity numbers. There are actually two more 
sheets to this report for activities 1-7 and 8-14. Note that activity 
21 is the sum of the first 20 activities. 
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ERDMAN a ANTHONY 

DAILY TIME REPORT 


LOCATIONt 

Zocfajfer, A/. Y 


NAME: 




WEEK ENDING: 

// -/5-63 


TITLE: 




OPER. 






OAILY 


HOURS 






1 TOTAL 


PROJECT 


NO. 
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T 


W 
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NO. 
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\ HEREBY CERTIFY that the tim« reported above is a true 
statement of my working time for this period. 


and complete 

I.U! 


Emolovti 










A and 
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DAILY TIME REPORT 


LOCATION 
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1 HEREBY CERTIFY that the time reporte 
statement of my working time for this t 


d above is true or 
>eriod. 

Rd - 


id complete 

.1.1 


y J Empoye^ /Signature API 



(ATTACH RECEIPTS WHEN OBTAINABLE FOR ALL ITEMS) 
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EXPENSE ACCOUNT |j 
2 Weeks ?/? /2 
Ending Ajf * " 

!■ a? v*m pSSBSmBtSOSBSSSm SSBS83S P* HOB W B 



Rochester, New York 




SHEET 



(Please print) 
Statement of 



CKDBY 



{/£N£TH IA/. 3gOl>VAJ 



Mailing Address 



Date 



Brief statement of itinerary, etc. , 
showing nature of expense items listed 
Day| in columns right 



MI 



Mileage 

and 
Transp. 



Lodg- 
ing 



Meals 



Other 
Itemize) 



Daily 
Total 



S 



M 



4.00 



6. SO 7. SO 



/AOO 



Job topocfior) 6/i '5/, 



4.0O 



4.0O 



w 



F. 



M 



T 



Job /vspetfroo 60/ 5/. 



4.DD 



3.0O 



loo 



w 



Jupp/'fS 



Z2-/4. 



Z.23 



2.2 b 



F 



I HEREBY CERTIFY that this is a true and 
complete statement of authorized expenses 
incurred by me in connection with company 
business for this payroll period. 



JZDD 



J0.5O 



2.23 



Employee 



Approved^ 

fry 



11/22/63 



E ROMAN AND ANTHONY, CONSULTING ENGINEERS 
PAYROLL. MASTER- JAPE 



O 







1 



Pay 



2. ISO 



Bonds 



.00 



Employee ^/Vamd [... 



JOHN J. JOHNS 



25.00 



.00 



.00 



1 



Income 



720.00 



26.12 



Gross 



129.60 



15.92 



4230.00 153.45 761 .40 



State 



95.29 



Persona/ 
£/pe/7S. 



.00 



ee*/s+<rrs 



5oc. 6ec. •#* 


Pension 


ufe 


/v.y.5. 
&iso(? 






■ ■ - — - — - — — — - — — " — 


120 30 1404 


.00 


.00 


.60 



2""4~.T88~3~ "'.00' 
KENETH W. BROWN 
95 1 3346 



-.oo~--~aa" 

5.50 .60 



-.THr"n3*0~.16 .00 1B5.0B "28.56' 
7705.92 174.00 1064.21 164.22 43 8.29 



1.00 



3 3_.?50_4_ _ J 00_ .00 5^40 _,01_2 1261.00 .00_J52.10 21.61 

"MARY'BENSON " " 6909.45 174.00 "813.16 108.81 515.68 
126 34 9929 .00 .00 .60 



299___ J ogg_o___ J pg___ J po___ J p5___^x)().o ^oo .,.0.0 ...oo ^00. 
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TABUS OF ADDRESSES 



The addresses shown in the table represent the address of the left roost position 
of each variable. 



ENTRY 
NUMB* 


MAN 
NUMB. 


REG. 
HOURS 


OVT. 

HOURS 


PERS. 
EXP. 


ENTRY 
NUMB. 


MAN 
NUMB. 


REG. 

HOURS 


OVT. 

HOURS 


r:rs. 

EXP. 


1 


19830 


18530 


17230 


15930 


.23 


19610 


18310 


17010 


15710 


2 


19820 


18520 


17220 


15920 


24 


19600 


18300 


17000 


15700 


3 


19810 


18510 


17210 


15910 


25 


19590 


18290 


16990 


15690 


4 


19800 


18500 


17200 


15900 


26 


19580 


18280 


16980 


15680 


5 


19790 


18490 


17190 


15890 


27 


19570 


18270 


16970 


15670 


6 


19780 


18480 


17180 


15880 


28 


19560 


18260 


16960 


15660 


7 


19770 


18470 


17170 


15870 


29 


19550 


18250 


16950 


15650 


8 


19760 


18460 


17160 


15860 


30 


19540 


18240 


16940 


15640 


9 


19750 


18450 


17150 


15850 


31 


19530 


18230 


16930 


15630 


10 


19740 


18440 


17140 


15840 


32 


19520 


18220 


16920 


15620 


11 


19730 


18430 


17130 


15830 


33 


19510 


18210 


16910 


15610 


12 


19720 


18420 


17120 


15820 


34 


19500 


18200 


16900 


15600 


13 


19710 


18410 


17110 


15810 


3:5 


19490 


18190 


16890 


15590 


14 


19700 


18400 


17100 


15800 


36 


19480 


18180 


16880 


15580 


15 


19690 


18390 


17090 


15790 


37 


19470 


18170 


16870 


15570 


16 


19680 


18380 


17080 


15780 


38 


19460 


18160 


16860 


15560 


17 


19670 


18370 


17070 


15770 


39 


19450 


18150 


16850 


15550 


18 


19660 


18360 


17060 


15760 


40 


19440 


18140 


16840 


15540 


19 


19650 


18350 


17050 


15750 


41 


19430 


18130 


16830 


15530 


20 


19640 


18340 


17040 


15740 


42 


19420 


18120 


16820 


15520 


21 


19630 


18330 


17030 


15730 


43 


19410 


18110 


16810 


15510 


22 


19620 


16320 


17020 


15720 


44 


19400 


18100 


16800 


15500 
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JGROSS. 



.EW.L 



DED 



OVT.H 



OVT.P 



NET 



FICA 



SWT 



EXP 



JOHN J. JOHNS 11858. 4443.75 

._-_l__.8Q.OQ_10.00__180.00___33-75__2l3-75--_7.75__38.if8___5.29-._-25.60 



.00__136.63 



T^rtirvrr^RoWN~rf859T 



~8o4b.9~6~ 



2_-_8Q.gg_25.gg___35.g4 .gg__335.g4 .gg__46_.27___7.14 6_.1g___31_.28__306.8i_ 



RARY" BENSON fl86o; 
3 80.00 5.00 260.00 



7169.45 

.00 260.00 .00 28.08 3.I3 10.01 6.34 __225.12 



REG.H REG.P GROSS FWT 

OVT.H OVT.P FICA 



SWT 



DED 



EXP 



NET 



IOh.jyo nS.Oi ^._808_.79 11?_.83 V ___4_U71 668. 56_ 

^.00 33.75 7.75 15.56 37.62 
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ERDMAN AND ANTHONY 

CONSULTING ENGINEERS 
ROCHESTER, N. Y. 



DATE 

PAY 

TO THE 
ORDER 
OF 



11/29/63 VM CHECK NO 11858 



JOHN J. JOHNS 



PAYROLL ACCOUNT 

EAST AVE. BRANCH 

TO GENESEE VALLEY UNION TRUST CO. 

ROCHESTER. N. Y. 




DOLLnKS ICENTS 



136.63 



s/6A//? razzes 



STATEMENT OF EARNINGS AND DEDUCTIONS 



REGULAR 


OVERTIME 


RATE 


REGULAR 


OVERTIME 


GROSS 






HRS. 


HRS. 


PAY 


PAY 


PAY 




EXPENSES 


80.00 


10.00 


2.250 


180.00 


33.75 


213.75 




.00 



FED. TAX 


F.I.C.A. 


STATE TAX 


HOSP. INS. 


BONDS 


SAVINGS 


PENSION 


N.Y. DISA. 


MISC. 


NET PAY 


38.48 


7.75 


5.29 


.00 


.00 


25.00 


.00 


.60 


.00 


136.63 



YEAR -TO-DATE TOTALS 



PER. END. DATE 


GROSS PAY 


FED. TAX 


F.I.C.A. 


STATE TAX 




PENSION 


11/22/63 


Wf3.75 


799.88 


161.20 


100.58 







__________________________ 

ranWTNDTN(_~FRTDAY~ 

_QY|MBER_22 t _1963 



-0 



Bonds savings 


HOSPITAL- 
LIZATION 


MISC 
DEDUCTIONS 


DISABILITY 


GROUP 
INSURANCE 


I 


1 


JOHN J. JOHNS T. 
•00 25.00 


.00 


.00 


.60 


*Q0 


25.60 _ 














KENETH W". BROWN 2. 
.00 .00 


.po__ 


_ .00 . 


.60 


5.* 50 


6.10 


~~~RARY~BENSON IV 
.00 .00 


9.1*0 


.01 


.60 


.00 


_ _ _ ._ ... _. ... 
10.01 




.00 25.00 


9.40 


.01 


1.80 


5.50 


M.71 


M.71 
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10 
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_ _ E ROMAN AND, ANTHONY, CONSULT I NG ENG.! JhlEERS 

PROJECT N O . 1 . 

_ 9/27/63 __ _Q 

LABOR EXPENSES, _____ S_ +V, TOTALS 

-.06 ,00 _PP -.0 6 

1 _ 1 80_, 00 ,00 _ PP_ _ 1 80,. 00 

37 3 6P_,PP_ _ _ ,00 .,00 3 60_, 00 

52 L_7__25 5.70* .00 152 .95 

6?_____ 33_*51 ,op _.00 ___39.,51 

63 _ J32_85 .Q0_ _.PP__ J32_,85 

97 115.15 „pp_ _00 115.15 



31h_.J0 5 ,7P„ _, 00 _ _ 9_8p_, 40 380_._f0. 



ERDMAN_ AND, ANTH&NY^_. CONSULT I NG . ENGJ Nf-ERS 

PROJECT NO, 1114,. 

-.5/27/63 

LABOR EX PEN S E S _ _ _ S_ +V_ , TOTALS. 

1 3 50.00 00 _0P 15PJD0 

_ _ _ _ 2 2_ 2J 6_.PP . 00. 00 21 _6_, 00. 

35 _*2. 25_ _ _ _ .00 _,00 if2_,25_ 

3§ 27.00 .00 27.00 5A.0P. 

45 25.00 _.0P .,00 .2 5..PP. 

6p_ JP7_.j25__ _.pp _.pp_ JP7.,_25. 

80 *L_li ,00 _00_ 8.14 



775_. 6^_ .00 27_.PP 802 , 64. 80 2.64 



Z4 



ERDMAN. AND. ANTHONY , CONSULT! NG_ ENGJ NEJERS. 



_ 3/27/63 





LABOR, 


EXPENSES 




TOTALS 


1 . 


97^.70 


5.70 


.00 


980.40 


2. 


444.01 


.00 


.00 


444.01 


_1 1 1 if. 


775_.64___ 


.00 


27_..00 


802.64 


1115. 


117.00 


'.00 


.00 


117.00 


__JJ52. 


J4.P_6_.8_8___ 


45.14 


_.pp 


J 452,02 


2201 . 


778.82 


22.37 


140.00 


941.19 


2202. 


140.02 


.00 


.00 


140.02 


2204. 


591.50 


.00 


270.73 


862.23 


„_2_205_. 


JP3_8_.3_6___ 


__.po 


.___77_.PP____ 


__.JJJ5_.3_6_ 


2206. 


295.25 


.00 


.00 


295.25 


2208., 


325.00 


.00 


44.00 


. __ 365^.00 


2212. 


_ _ 23.75 


1.86 


.00 


_ 25.61 


2213. 


1292.64 


.00 


53.25 


1345.89 


___22J4. 


33.18,_Q8___ 


___5_1_6J7___. 


_.00 


___ 3-834.25 _ 


2215. 


257.38 


.00 


.27.00 


284.38 


2218. 


.nn 


Q 


nn 


«/ • ^* • - 


2221 * 


165.19 


.5.55 


-00 


170.74 


__.2222_, 


180^73— 


_.-00____ 


_,.00____ 


J 80. 73 


2223-. 


32..55 


.00 


.00 


-32.55 


...3102. _ . 


21.25- 


....00 . 


. _50.63 


71 .88 


_ .310 7-. _ . 


729_78-__ 


_ 102.59- 


. __._40_.63 ___ 


.873 . 00 .... 


, 3108,. 


4.25 


.00 


.00 


4.25 


_ 59105. 


357.16 __ 


11 .40 


. _ _-.oo 


368.56 


60107. 


1 1 1 . 80 _ _ 


.00 _„_ 


. ...00. 


_.. 111.80 


.60151. 


1 666.80 _ 


... 61 .39- 


10.00 


... 1738.19 


_ 61102 


1 556.76 .... 


.00 


Jj3_4_-0_0 


1720.76 



01 fc, 



P£CU£CT STATUS 

"'" " ~~"ERDMAN AND ANTHONY, CONSULT] NG" ENG i NEERS 

10/11/ 63 " 







LABOR 


EXPENSES 


s.+u 


TOTALS 




2203. 




538.20 
.00 
538.20 


.00 
,00 
.00 


.00 
.00" 
.00 _ 


538.20 
.00 

538.20 


... . ^v^.j* 


62208. 




2206.13 


172.94 


128.50 


2507.57 








.00 

2206.13 


,00 
172.94 


.00 

128.50 


.00 

2507.57 




62213. 


— 


4293.42 
.00 

4293.42 


544.00 
.DO 
544.00 


* _ «_ . 

.00 

180.89 


5013.31 
.00 

_______ 5018.31 


' 


3101 . 


— 


313.31 
.00 

313.31 


44.45 
.00 

44.45 


.00 

" v\ 

.00 
.00 


357-76 
.00 
357.76 


- - 


_3J0;2_.__ 




2053.95 

_ « • s\ 

1384.89 
3*68.84 


160 ,.24 
27.9- 
' 183.16 


50.63 

9.00 ~ 
59.63 


2264.82 

— — . . . r> * 

1 421 .81 
3686.63 





3103 v 




1179.89 
.00 
1179.89 


74.33 
.00 

74.33 


76.00 
.00 
76.00 


1330.22 
.00 

I33O.22 




31.0_4_.___ . 




61.88 

.00 

61 .83 


.00 
.00 
.00 


.00 

.00 

.00 


61 .88 

.00 

61 .38 


-•— - 






760.47 


193.40 


55.93 


1014.80 








.00 
760.47 


.00 

198,40 


.00" 
55.93 _ 


.00 
1014.30 


-~ ; ' - 


3106. 




H838.OI 
,00" 
11838.01 


2915.55 
.00 

2915.55 


112.54 

"" "TOO" 
112.54 


14866.10 
o 00 
14866.10 




3107. 





746 23 

r/y 

796.66 


102.59 

102.59 


40.63 
40.63 


889 50 

C p 2 

939-38 




3108. 




4.25 


.00 


.00 


4.25 






— 


- " .00 
4.25 „ 


.00 
.00 


" .00 

.00 


"~~ .00 
4.25 




3 _ _62?PA.____ 
-7 




2484.90 
48.00 

2532.90 


24.15 
.00 
24.15 


214.74 _ 

.00 
214.74 


2723.79 
48.00 
2771.79 






.1 
1 


520352.50 
17667.45 
538019.95 


46879.50 
745.53 
47625.03 


58380.27 
564. 1 1 

58944.38 _ 


1 625612.27 
18977.09 

.__l_..644589.36 


1 625612.27 
18977.09 
: V*MQ 




1 
1 


556023.36 
18998.65 
575022.51 


43152.42 

745.53 
48897.95 


60543,30 
564.11" 
61112.41 


1 664724.58 
20308.29 
1 685032.87 


] 664724.58 
20308.29 
1 685032.87 



V 



a 



o 



ERDMA_N_ MP. _ MTHO_NY^_ CONSUL TJ NG. JE_NGJ_N_E_E_R_S_ 



_PRC J_ECJ_ N0,_ 61 J5J. 



6/21/63 





2 




I -> 


91 


56 


n 

V 


o 


n 


TAT 
1 V 1 




1 

I 


50.00 


.00 


46. 12 

~V • 1 »- 


14.25 


.00 


-00 


00 

. V V 


no 

. V V 


1 1 n 37 


1 

i 

lift }1 

1 LVOV. 


2 


.00 


.00 


.00 


.00 


7 .22 


.00 


.00 


.00 


7 99 


i 

7 99 


3 


.00 


.00 


4.8^ 


.00 


.00 


.00 


00 




4 85 


r • Op 




.00 


.00 


.00 


.00 


00 


-00 


00 


on 


nn 

. vv 


nn 

t uu _ 
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INTRODUCTION 



The management of ACF Industries, Incorporated realized that manual methods for 
accessing personnel information were inadequate both with respect to speed and 
accuracy. Many requests for summarized information either could not be answered 
or became a costly venture in man hours alone. As an attempt to remedy this prob- 
lem, a punched card system was designed to hold personnel information in coded 
form. This card system proved inadequate in speed, type and amount of information 
stored, and updating capability. With this system as a background, SPIRE was 
developed . 

SPIRE is an acronym for "System for Personnel Information REtrieval . 11 This system 
was designed for and implemented on a magnetic tape 60K IBM 1620 and has been 
operational since February, 1963. Since that date, no revisions have been neces- 
sary in the type or format of the information. 

CONTENT 

The information content of SPIRE was specifically designed for ACF's use. Personnel 
information was carefully analyzed and organized to determine what information had 
been recorded by the previous system and what part of this information had been use- 
ful. By categorizing the useful information, it was possible to devise coding techniques 
which made each item unique in the system and thus recoverable. 

The information in SPIRE is broken into the following five categories: 

1 . Personal — name, date of birth, sex, date of hire, marital status, social security 
number, etc. 

2. Previous Experience — job code, position title, time on fob, extent of super- 
vision, etc. 

3. ACF Experience — complete file on all pertinent happenings at ACF with the 
employees status (job title, grade, salary, department, etc.) recorded at the 
time of the action. 

4. Education — includes all completed levels of education and the extent, major, 
minor, and university, if applicable, for each level. 

5. Military Experience — contains the same information as the previous experience 
category with the branch of service added. 

To obtain this information in coded form for all employees, supervisors interviewed 
their subordinates and coded their education and experience backgrounds. The ACF 
experience involved so much detail that it was necessary to code each part of the 
employee's experience using the permanent personnel records. As new employees are 
hired, their backgrounds are coded and added to the master file. 
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III. APPLICATIONS 

Although SPIRE was developed as a general purpose information system that could be 
used for the generation of virtually any report involving personnel information, one 
application appears to be the most notable, that of in-plant recruiting. Prior to the 
development of SPIRE, to fill a vacant position from within required the personnel in 
Salary Administration to have an intimate knowledge of the background of all salaried 
personnel . This, of course, became a virtually impossible task as the number of em- 
ployees increased. Qualified people could easily have been overlooked, and the 
preparation of readable resumes for the managers became a monumental task. SPIRE 
now assures that no qualified individual is overlooked and automatically prepares 
their resumes . This resume, as printed, contains virtually no coded information, thus 
providing management with a concise, complete, easily read and understood document 
on any salaried employee. 

Of course, qualification is a matter of definition. A professional analyst takes the 
vacant position and determines the background necessary to qualify an individual . The 
analyst specifies acceptable work experience and the minimum years of experience in 
related fields. He also specifies the acceptable areas of education and the minimum 
acceptable level . From the years of experience and the education level specified, a 
minimum number of "points" are computed by the system which will qualify an individual 
for the vacant position. Individuals may accrue the necessary points in any combination 
of related experience and education. Thus, if an individual had a masters degree in an 
applicable field and only a bachelors was required, he could qualify for the position 
with less experience than that specified. The equating of education and experience 
used in the SPIRE system is based on the ACF position rating plan. The analyst also 
specifies the acceptable salary grades for the search. This keeps senior men from "drop- 
ping out" for Junior fobs for which they would be qualified. In addition to the precise- 
ness with which this search for qualified employees is conducted, a considerable savings 
in both time and money exists. 

The normal time to complete a search for any given position varies between six and 
eight minutes. This same fob accomplished by manual methods would take about six 
man hours. 

Another application which has saved a considerable number of man days is the gener- 
ation of the quarterly Merit Review Report. This report is by department and contains 
pertinent information concerning past salary increases for all employees who are to be 
reviewed during the next quarter. The normal time to complete this run is about 15 
to 20 minutes as opposed to approximately 32 man hours. 

Other current applications of SPIRE include: creation of summary reports on salary 
increases for any time period, skills inventories, projections for salary budgeting, 
preparation of data for salary surveys, current salary status reports, salary distribution 
studies, and the generation of resumes for specific purposes. 



IV. COMPUTER TECHNIQUES 



SPIRE is centered about a master information tape. The information stored on this tape 
is first punched into cards. Each card contains all the pertinent information about a 
single happening in the employees past. Every card entering the system is identified 
by the category number and the employee's permanent number. The remainder of the 
card contains identifying codes for the data as well as alphameric descriptions. 

The master tape is blocked by employee. All the information about any employee forms 
a single tape record. Within this record the information is held in card images sorted 
by category and date where applicable. This means that the maximum number of records 
to be read or bypassed will be the number of employees in the system. The employee's 
current salary grade is carried as part of the first card image so that for many applications 
a decision can be made as to whether or not a particular record is of interest to the search 
before actually performing a resume scan. This is particularly useful in the recruitment 
appl ication . 

The accuracy of the information contained in the master file is, of course, quite critical. 
To provide this accuracy, a great deal of time and effort has been placed in the program 
used to update the file. A total of 47 different errors are uniquely detectable in the up- 
dating data. This program allows total resume replacement or deleting, individual card 
replacement or deletion within any resume, and the addition of any new information to 
any given resume. 

In general the customary way the master file is used in supplying information, other than 
complete resumes, is to scan the master tape and write a summary tape of the required 
information. The information on this summary tape may then be sorted or interrogated 
in any prescribed manner. Several of these codes that prepare summary tapes have been 
written. In the majority of special requests, one or more of these summary tapes can 
provide the necessary data. In these instances it is usually possible to write a program 
using ACF's magnetic tape FORCOM system which will interrogate the summary tape 
and prepare the necessary report. This provides a quick, economical way to satisfy 
demands made on the system . 

V. CONCLUSIONS 

SPIRE has proved a complete and useful system. It has resulted in large cost savings 
and has been well accepted by management. This acceptance is born out by the con- 
stant increase in requests for information from this system. To date there has been no 
request that could not be satisfied because the data had not been included in the system. 

SPIRE was designed specifically for ACF's needs by ACF personnel and, as such, it may 
or may not be applicable to any other organization. This system is operational and in 
fact is far exceeding the original expectations in speed and versatility. 
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A Computer Program for the Calculation of 
Prime Implicants from a List of Boolean Minterms 



I Introduction 

Boolean algebra has proven to be a most useful tool in the logical 
design of digital circuits. In general, the method employed involves 
four phases: 



1- definition of the problem (verbal statements) 

2- translation of verbal statements into Boolean equations 

3- simplification of the equations 

4- implementation of the simplified equations 

Step 3 - Boolean simplification - can be done in a number of ways. One 
method is attributed to Quine. ^ The program to be discussed here enables 
the 1620 computer to simplify Boolean functions in a manner very similar 
to that proposed by Quine, obtaining as a result all the prime implicants 
of the original function. 

2 

II The Quine Method 

To simplify a Boolean function by the Quine method, it is necessary 
to express the function as a sum of minterms. The N- letter minterms are 
then systematically compared with each other two at a time, and all pairs 
differing in the state of only one letter are combined according to the 
identity: 



where A represents the letter eliminated and X is all other letters in 
the minterms involved. 

After all possible comparisons have been made, there is in general 
a list of terms having N-l letters. If any minterm did not combine at 
all during this process, it is a prime implicant. Minterms are checked 
as they are found to combine, to facilitate spotting of "non- combiners". 

Comparisons continue among the N-l letter terms, with the restriction 
that only terms containing the same letters can possibly combine. The 
N-l letter terms fall in as many as N different sets, and further combination 
is possible only within these sets. 



X A + X A - 



X 




lot 



Successful combination of N-l letter terms produces N-2 letter terms. 
Again, uncombinable terms are prime implicants. 

This process is continued until no further combinations are possible. 
Uncombinable terms at any level are prime implicants. The Quine method 
then goes on to select the simplest set of prime implicants that is 
logically equivalent to the original function. This last phase is not 
handled by the program to be discussed - its output is simply a list of 
all prime implicants. 

The reader is referred to Reference 2 for a much more detailed 
description of the Quine method, along with numerous examples. 

III The Program - General 

The program to be discussed was written by Mr. H. Huhta of the General 
Electric Co. as part of his M.S. thesis project at Union College. The 
required inputs are: 

1: KV - the number of variables (cannot exceed 9) 

2: KN - the number of minterms in the function (cannot exceed 60) 

3: the list of minterms, coded in octal form (see IV A) 

The resulting outputs are the prime implicants, coded as described in IV D. 
The FORTRAN-0 language was used. The following paragraphs point up 
various aspects of the program organization. A detailed example is included, 
together with Veitch diagram verification of its correctness. 

The Union 1620 has only 20 K memory. If both KV and KN are large, the 
program may not fit. To insure the ability of the machine to cope with any 
9-variable problem, a second program was written to perform a single level 
of reduction. Thus the output, with 9- letter minterms, would be 8- letter 
terms which would then be reentered to produce 7- letter terms, etc. The 
second program will not be discussed further in this paper. 

The first program has been shown to handle an 8-variable, 24-minterm 
problem successfully. The exact limits are not yet known. 

IV Program Details 

(A) Input Coding 

Input minterms are coded as fixed-point octal numbers. Each 
octal digit specifies the state (i.e. - complemented or not) of three 
variables, hence three octal places will suffice for a 9-variable 
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problem. For example the 5- letter minterm A B C D E would be entered 
as 25, using the conventional scheme of representing an uncomplemented 
letter by 1 and a complemented letter by - but converting to octal (25) 
rather than decimal (21). 

(B) Combinability 

Terms at any level are combinable if: 

1- they contain the same letters 

2- all letters but one are in the same state 

Satisfaction of requirement 2 is detected by a two-step process. First, 
terms must differ in only one octal place, and second, the differing 
octal digits must differ in only one binary place. The fact that the 
second of these steps can be carried out at the octal level follows 
directly from the relationship between the binary and octal number 
systems. A table L (I) contains the necessary information. This 
43-entry table relates Reference Octal Digit M and Compared Octal Digit 
N. It is addressed by the number N + 5M, and contents of each address 
is either 4,2,1 or — to indicate the weight of the differing binary 
place. means: "These two octal digits do not differ in just one 
binary place, hence they are not combinable. 11 

Terms at any level in the process are always arranged in ascending 
order prior to comparison. Thus if M is larger than N, it means that 
the terms being compared differ in at least one other octal position - 
hence they do not combine. The entries for all cases of M^N are there- 
fore in the table. 

(C) Set Identification 

As mentioned in II, comparisons at all levels beyond the minterm 
level take place only within sets . The program must therefore generate 
a set identification number ( S I N ) to identify each combined pair. 
The SIN increment generated by any particular comparison is simply 
the octal weight of the letter eliminated. As terms pregress through 
successive levels, an up-to-date S I N is maintained by cumulative 
addition of the weight of each eliminated variable. 

(D) Example - to determine the prime implicants of the function: 

f » m 3 + m 4 + m 5 + m ? + m 9 + m 1() + m^ + 



MO 



(1) List minterms in ascending order in octal system- 

m 3 —» 0003 



m 



4 



0004 ^ 



m c — > 0005 
m 7 — * 0007 
m 9 0011 
m 10 — » 0012 

m 12 ~* 0014 
m 13 — > 0015 

(4 digits are shown, because this is the way FORTRAN represents 
fixed point numbers internally. Leading 0's need not be typed.) 

(2) Investigate all possible combinations. Start with 3, and compare 
each term with all higher- numbered ones. 

ex: 3 - 011 three binary dif ferences , \ not combinable 
4 = 100 



3 = 011 two binary dif ferences t .not combinable 
5 = 101 



3 = 011 combinable ; the most significant binary place 

7 - 111 is the one that combines 
The decision as to whether or not a pair of octal digits is combinable 
is made by referring to the table L (I) , as discussed in IV-B. In the 
3-7 comparison just described, the most significant place is the one 
that combines, so address 22 (7+5X3= 22) of the table contains the 
number 4 - the octal weight of the eliminated variable . 4 now becomes 
the S I N of the combined term. It is attached in the 5th place of the 
lower-valued of the two combined minterms, and the result stored to 
await processing at the next level (5 places are now possible because 
the combined term is given a floating point FORTRAN name) . 

This process continues until all possible combinations are tried. 
Minterms involved in one or more comparisons are tagged by putting a 
9 in the 4th (most significant) place. If any minterms have not entered 
into any comparisons when the processing of level N is complete, they 
are printed out. They are prime implicants. 
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The complete level 4 results appear below: 



T-r, -l «-4 Q 1 
JLIl 1 UXel 1 

Level 4 


Level 4 


Tffln c f A y*t*o /I 

to Level 3 


fl otnhi tip i" i on of 




7UUJ 


£.000 ^ 


000^ 9nH 0007 


0004 


9004 


10004 


0004 and 0005 


0005 


9005 


100004 


0004 and 0014 


0007 


9007 


20005 


0005 and 0007 


0011 


9011 


100005 


0005 and 0015 


0012 


0012 


40011 


0011 and 0015 


0014 


9014 


10014 


0014 and 0015 


0015 


9015 







At this point, 0012 (minterm ten) is outputted - it is a prime 
implicant. Output coding will be discussed in paragraph (3) 

Processing of level 3 starts by arranging the transferred terms 
into sets having identical S I N's, and putting the terms in each 
set into ascending order. The result is: 



10004 
10014 



20005 



40003 
40011 



100004 
100005 



:> 
;> 



S I N - 1 

S I N - 2 

S I N * 4 

S I N - 10 



Digit-by-digit comparisons are now made within each set, just as 
before. Since the S I N's of two terms being compared are the same, 
comparison need be made only of the four least significant digits. 
(The S I N is actually dropped before comparison, by returning to use 
of fixed-point variables.) 

The terms with S I N ■ 1 combine, since they differ only in the 
second least significant octal place, and and 1 combine (Table L (I)) 
to eliminate the variable of weight 1. The new S I N is therefore 10 
(generated at level 3) + 1 (generated at level 4) « 11. The term 
transferred is therefore 110004. Note that the SIN increment of 10 
indicates that the variable eliminated had weight eight. 



There is only one term having S I N = 2, hence 20005 is outputted 
at the conclusion of level 3 processing. Similarily, the two terms 
with S I N = 4 do not combine, hence both 40003 and 40011 are outputted. 



The terms with S I N = 10 combine, since 4 and 5 combine to 
eliminate the variable of weight 1. The new SIN would be 11, and 
the transferred term 110004 - identical to that derived from the 
SIN =1 set. The program avoids this duplication by a means that 
will be discussed in paragraph 6. 

Since only one term was transferred to level 2, it is of course 
a prime implicant. The problem is now complete, except for the 
decoding of the computer output. There are five prime implicants, 
so far identified as follows: 

At level 4: minterm 00012 (note that the S I N of all minterms 



At level 2: term 110004 

(3) Output Format 

Three fixed point numbers are printed out for each prime 
implicant. They are: 

1- the minterm from which the P.I. was derived; 

2- the final S I N, to tell which letters have been 
eliminated; 

3- the level, to tell how many letters will appear in 
the P.I. (this is a check on the others; it contains 
no new information) 

If two P.I.'s come from the same set, the SIN and the level 
will be printed only once. 

The print-out for the example will therefore appear as 



All three of these terms are prime implicants. 



is 0) 



At level 3: terms 20005 
40003 
40011 



follows: 




12 




5 



2 



3 



3 



11 



7 



(4) Output Decoding 

The output format of (3) can be translated to literal form as 
follows (using the first prime implicant to illustrate) : 

1- write letters for variables: A B C D 

2- translate octal minterm to binary: 12 — > 1010 

3- translate S I N to binary: 0000 

Now each in binary SIN indicates a letter present in the prime 
implicant (in this case, all four). The binary minterm now indicates 
complemented variables by 0, true variables by 1. The prime implicant 
is therefore A B C D, which is minterm ten, in the usual decimal 
notation. The output level -4- verifies the fact that the term should 
contain four letters, which it does. Similarly, the other prime 
implicants decode as follows: 



A B C D 



10 1 
0010 J 



A B D 



11"] 
10 J * 



A C D 



11 
4 



1 1~T 
10 J * 



A C D 



4 

11 



10 0"] 
10 1 J * 



B C 



(5) Check of Results 

The solution of this problem by Veitch diagram is shown below. 
It is readily verified that the computer solution is correct. 

Prime implicants: A B C D (™±q) 



B 



1 

n 




* 


1 

H 






1 

7 


i. 




it 


i. 


/ 


s 


i. 


X 


* 



D 



B C (m, + m_ + 
4 5 

A B D (m 5 + m 7 ) 
A C D (m^ + m^) 
AC D (m 9 + m 13 ) 



m 12 + m 13 } 



nit 



(6) Elimination of Duplication 

It is inherent in the Quine system that identical combined 
terms will arise from different series of comparisons. For 
example, a term A C could arise by eliminating B from ABC and 
ABC, or by eliminating D from A C D and A C D. The program 
eliminates such duplication at all levels by the following 
procedure : 

When a combinable pair is found, the old SIN and the weight 
of the variable just eliminated are compared. The term is 
transferred to the next level of comparison only if the old 
S I N is less than the weight of the variable just eliminated. 
For example, in the problem just illustrated, the two terms 
with S I N 8 1 combined to eliminate the variable of weight 10 
(decimal eight ) . Since 1<10, the term was transferred. 
Later on, the two terms with S I N = 10 combined to eliminate 
the variable of weight 1 (which would have produced the same 
term). This term was not transferred because 10 > 1. 
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SUMMARY 



Critical speed, stresses, and bearing reactions can be calculated on a 
digital computer for a general two bearing shaft, which can have each 
bearing located anywhere on the shaft, any number of cross sections, 
variable loading, and any length. 

Besides eliminating the tedious labor of the calculations, the computer 
provides flexibility. A number of calculations may be rapidly made to 
optimise design or to check out application variations on a standard de- 
sign. 

Basis of the Calculations 

In general, the critical speed calculation is based on Raylelgh's method 
using the following formula 



where N c is the first critical speed, rpm 

W g is the load at the s position excluding external loads without 

mass, such as belt pull, pounds. 
y s is the static deflection of the load at the s position due to 

the loads V with the loads reversed between bearings. 

The static deflections and bending moments are found by using the same 
principles of mechanics which are the basis of graphical determination 
of beam deflections. Numerical integration is used to relate load, shear, 
bending moment, slope, and deflection. This integration, although mm* 
erical, gives exact areas as long as the loads and moments of inertia 
are constant over the increments of shaft length. 

The bearing reactions are found by taking moments about each bearing 
equal to zero in accordance with statics. 

The shaft bending tensile stress, torsional shear stress, and combined 
shear stress are calculated according to the usual equations of strength 
of materials. 

Method 



The method is illustrated in detail by an example calculation. Briefly, 
it consists of dividing the shaft into increments of length, determining 
the load and shaft moments of inertia in each increment and the computer 
calculates critical speed, stresses and bearing reactions. If there is 
only shaft weight in a number of increments, the O.D. and I.D. of the shaft 
are sufficient input data and the computer will calculate the load and 
moments of Inertia. 

Other Factors in Critical Speed 

There are other factors affecting critical speed, such as, shear deflection, 
gyroscopic effects, bearing length, and bearing flexibility. Bearing support 
and oil film deflection may be accounted for by adding the bearing deflection 
to the deflections used in calculating critical speed. 
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IHTftOPUCTIOK 



Critical speed, stresses, end bearing reactions for a general shaft may - 
be calculated on a digital computer by utilising numerical integration. v 
The general shaft has each of two bearings located anywhere on the shaft, 
any number of cross sections, and any number of loads. 

Besides eliminating the tedious lebor of the calculations, the computer 
offers flexibility. It allows the calculation of a number of variations 
in order to optimise design or evaluate application variations of standard 
designs. 



THEORY 

l*yioinh's Bnorgy Method (4) and ffl 1 



Ray 1 sigh's energy method may be used to determine lowest natural frequency 
of a system. The method consists of assuming a configuration for the system 
which will approximate the maximum amplitude of the fundamental mode or low* 
est natural frequency. Based on this configuration, the maximum potential 
energy is calculated at maximum displacement (zero velocity) and equated to 
the msxlmum kinetic energy at the system equilibrium position (maximum velo* 
city). This relationship may then be solved for the lowest natural fre- 
quency. 

The dynamic deflection curve at maximum displacement for laterel displacement 
of a beam may be assumed to be very close to the static deflection curve, 
when the configuration due to static loading of a beam with its own loads is 
used for calculating the fundamental frequency by ftayletgh'a method, the fre- 
quency is within the accuracy required for most engineering-calculations. 

Since Rayleigh's method consists essentially of equating inertia forces of 
the rotating masses to the restoring forces of the shaft, external forces 
which do not have mass, such aa belt pull, must not bs included. Eayleigh 
found that the system vlbretes in a manner which makes the frequency a minimum. 
This means that the inertia forces of the masses must act in directions causing 
maximum deflections. In a two bearing shaft the static loads between the bear* 
Ings must be reversed to produce maximum deflection and minimum frequency. 

The potential energy Involved in Rayleigh's method is the work done on the 
beam in moving its losds from the equilibrium position to the dynamic do* 
flection curve at maximum deflection. (This is assumed the same as the static 
deflection curve without, external loads and with loads reversed between 
bearings.) 
Thus: 



i/2£lv.| 



The maximum velocity of the weight V, with amplitude ys in harmonic motion 
is pys * Hence, the maximum kinetic energy of the system is 



2 g m* 



1. Mumbers in parenthesis designate references at end of paper. 

Hi 
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The symbols in the previous aquations are: 



P.I. 


is 


potential energy, 




is 


the load st the s position, s * 1,2,3,— 


y$ 


is 


the deflection of the load at the s position. 


K.E. 


is 


the kinetic energy. 


g 


is 


the acceleration of gravity. 


P 


is 


the natural circular frequency. 



The law of conservation of energy requires that the potential energy equal 
the kinetic energy in undamped systems. Solving for p* gives the fol loving 
result: 

p 2 - s / S |w «v.j 
V £l».y. 2 i 

The first critical speed of a shaft may be considered to occur when the 
circular frequency (G>) of shaft rotation equals the natural frequency of 
lateral vibration of the rotors (considered as a beam). The rotating shaft 
is deflected by the small unbalance forces in the rotating system, which 
cannot be eliminated entirely in actual rotors. These exciting forces have 
a frequency in a lateral plane equal to the circular frequency of rotation. 
When the exciting forces have the same frequency as the natural frequency 
of lateral vibration, resonance occurs, which usually causes undesireabie 
excessive vibration of the machine. Thus: 



N c is the first critical speed, rpm. 

W t is the load at the s pos i t ion exc luding external loads 
without mass, such as belt pull-, pounds. 

y $ is the static deflection of the load at the s position 

<W g ) with the loads reversed between the bearings, inches. 

Deflection 

The static deflections of the shaft system may be determined by utilising 
the principles of mechanics which are the basis of graphical determination 
of beam deflections as shown in Figure 1,(13), (10), and (3), 

Bearing Reactions 

The bearing reactions are found by taking moments about each bearing equal 
to sero in accordance with statics. 
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Simply Supported Beam 
t is the load per unit length- H ||f | 1 1 1 



MI1IIHLI 



tt, is the bearing reaction at the left end 
converted to load per unit length. 

r 2 is the hearing reaction at the right end 
converted to load per unit length. 

w is the load per unit length including 
the bearing reactions. 



Load Diagram 




-A Hdx 



V is the shear. 



V«^w dx * Area under the w curve; and * w 



V Shear Diagram 



M is the bending moment. 

M«Cv dx 4- external moments; and M = V 



M Moment Diagram 



M 



E is the modulus of elasticity. 

I 4 is the moment of inertia of the bean 
cross- sectional area with respect to 
the axis in the neutral plane. 




Area under the V curve -Y M ext# " ! j , 

(Most shafts have zero external bending moments.) ~1 r~ajc 




Hkdx 



•e- is the area under the (M/EI^) curve. 
S + Cg^CM/EI^) dx = 9- ; and ||U |M_ 



s*e>-c 



s » 



and 

S is the slope and C s is the integration 
constant. 




H kdx 



S is the area under the 6- curve. $ 
y + C y * £ S dx -j^-Cg) dx*S-C s x; and |2*$~0 B 
y*S-C s x-C y ; and.^C«^U%; bo ||»e- 
y is the deflection, and Cy is the constant. 




FIGURE 1 - BASIS OF 3-RAPHICAL DETERMINATION OF DEFLECTIONS 
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The shaft banding tensile a trot , torsional sheer stress and combined shear 
strata ara calculated according to tba usual equations of Strang th of 
materials (13). 



S 



bs 



where all valuas ara at taction s of tha shaft: 



a bs 

*• 

c 

s t. 
T . 

c 

Sc. 



is banding tent 11c stress, pal 

is banding moment, in lbs. 

is tha radius to tha outer shaft surface divided by 
the moment of inertia of the cross sectional area 
vlth respect to radial axis, in. "3. 

is tha shaft torsional shear stress, psi 

is tha torque in tha shaft, in. - lbs. 

is tha radius to tha outer shaft surface divided by 
the moment of inertia of the cross sectional area 
with respect to tha longitudinal axis, in.* 3 

is the shaft maximum eonbinad shear stress, psi 



The integrations shown in Figure 1 may be performed by numerical means 
with a computer, numerical integration techniques may be approximate, 
(8) and (15). However, this numerical integration is exact if the load 
and the moment of inertia are constant over the increment of x, (A x) 
and the following formulas are used. 

Taylor's Formula (8) 

From calculus Taylor's formula with remainder follows: 

*<•> ♦ JL14 (x . «) + ♦ t f!^ (x . „»♦ f <^i . » ♦ i 

1 ' « . (n ♦I): 



f(*) 



In Figure 2 if x • b then (x - a) • Ax and if x • c then (x - a) • Ax 



3t 



S - '/2 



s ♦ Vt 




Ax is the increment of length along the beam, 
s is tha station line in the center of Ax. 
Subscript s designates the value at the station line. 
Subscript s • ^designates the value at the left of Ax. 
Subscript s + '^designates the value at the right of Ax. 



FIGURE 2 
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Shear Formulas 

V ' ■ dV ■ v «■ (constant over &x)» w e ; V" « 

s 8 

Bending Moment Formulas 

m' « dM » V; M" « V 1 - w } M< 3 * « 
dx* 8 

Slope Formulas 

dS » M n and since S « 9 - C g ' ? dS « d£ 
dx EI Z dx dx 

where EX is constant over A x 

(4) 

- d© « JJ_ ; 9" « j£ • J/_ ; 9 % ' m M" • v $ ; 9 • 

dx Bis BI z KI z KI z "IT 



Deflection Formulas 



dy » S and since yiz & - C g x - Cy d£ » d S - C 

ds dx dx 



d_S - S + C 8 » 9 therefore S' * • i 8" - d * " JL 

dx EI^ 

& W » 8" » _V_ ; - 9 < 3 > - w. , S (S) - 



EI 



zs 



s 



£S *- / 



■zs - 4 

384 8 
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Error 



Therefore, the only error introduced in the integration is due to the approxi- 
mat ion which is necessary to keep the loading and moment of inertia constant 
over the increment of x, ( Ax). If desired, the amount of error can be 
determined by approximating the load vith high values and the moment of in* 
ertia with low values for one calculation. A second calculation la then 
made by approximating the load with low values and the moment of inertia 
with high values. The total error will be less than the difference between 
results of the two calculations. The error may be reduced by making further 
calculations with smaller increments of x. This is not required for most 
engineering calculations with the usual safety factors* 

MET HOP 

The calculation method is illustrated by using the shaft shown in Figure 3a, 
This is a compressor rotor mounted on a stepped shaft with two bearings. Be* 
cause the loads need to be regarded differently, the calculations are done in 

two parts. 

Example Calculation - Fart 1 

First, the bearing reactions and stresses are determined including external 
loads without mass and with loads in their normal direction (downward) as 
tabulated in Table 1. 

The input data for Fart 1 of the calculations is described below i 



Column 



9 
11 
12 



Symbol Description 

&x Tha. increment of lengths 

Station numbers of bearing locations, 

a Station line number in the center of Ax, 

Ji B Load per unit length including external loads 

without mass in normal d ire c t ion . Concentrated 
loads must be distributed over the increment Ax 
Actual shaft and rotor loads are more often dis 
trlbuted than concentrated, 

C/I,f 3 See Fage 6. 

T. See Fage 6. 

See Fage 6, 



8 

The calculations made for Fart 1 are described below s 



Column 



Symbol 
ra 2 



mi 



/*3 



/ 



Description 

Moments about bearing Z station line as follows 
m m Jl [ s(of bearing 2) - s] calculated at 

e» A 

each station line, lb/in x station number. 
Moments about bearing 1 station line as follows 
°1 • ^s [»(of bearing 1) - *1 calculated 

at each station line., lb/in x station number. 

Fage 6 



123456789 10 















a* Rotors 



m 




d, Combined Shear Stress. Part 1 

x 



-10,000 



e. 3ending Moment. Part 1 





001 



f. Deflection. Part 2 
FIGrURE 3 - EXAMPLE CALCULATION DIAGRAMS 
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SS&SBL 



r, 



Description 
searing reaction 1, lbs/in 



•(at bear lag l)»s(at bearing 2) 
Bearing rcectiee 1, iba. • r^ ( & x) 



• 


r 2 


letting reaction 2, lbs /in 

•(at bearing 2) - » (at bearing 1) 




*2 


Bearing reaction 2, lbs. R 2 • *2 < A *) 


5 


*a 


Load ear unit length including bearing reactieni 
w § • t 9 except at bearing 1, w § */ # ♦ r^ 

and at bearing 2, w^ » ^ ♦ r^ 


6 

7 


* 

*» +Vi 


See Fagot 6 and 7. 
See Fagee e and 7» 


8 


*. 


See Pages o and 7* 


10 


"b. 


See Fage 6 


13 


s t. 


See Page 6 


14 


S c. 


See fage 6 



Rxample Calculations r fart ft 

Second, the critical speed it determined axe led lag external loads without mass 
and with the leada between bearings reverted (upward) aa tabulated in Table 2* 
Note the shrink fit of the compressor rotor and the coupling hub nay be con* 
sldered to increase the shaft diameter when calculating critical speed, (14). 

The input data for Fart 2 of the calculation is described below: 



Col 



Ax The Increment of length. 

Station numbers of bearing locations, 
s Station line numbers in the center of Ax. 



i. 



Load par unit length excluding external loads 
without mass, Concentrated loads must be dis* 

tributed over the increment A x. Actual shaft 
and rotor loads are more often distributed than 
concentrated. 



12 



llyg BXyg at station s. See Fags 5, 
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The calculation* aade for Part 2 ara described below i 



3 
4 



7 
8 
10 
11 
12 
13 



14 



•» ♦ V* 



Determination of y* froa S g 



Sea tart I, (O 
Saa Part 1. 
Saa Part 1. 
Sao Part I. , 

Saa Part 1. The indue ion of bearing reaction* 
■ay be quae t toned a iace they are external loaea 
without meat, However, It nay be shown that 
the inclueion of bearing reaction* ia necessary 
to make the integration V • J* wdx correct along 
toe length of the ehaft. 

Load, excluding external loaea at atation auabet 
a, Ibo, W # - Jt B ( Ax). 

See Pagan 6 and 7. 

See Page* 6 and 7. 

Sea Pago* 5, 6 and 7. 

See Page* 5, 6, and 7. 

Deflection* See following, 

Abaoluto value of coltmn 6 -x column 12 at 
atation s, 

Abaoluto value of col tan 13 x eoiuaxi 12 oft 
atation ». 

Pirat critical Speed. Son Page 3* 



It wa* shown (Page 5) that 

y • o • C *• • C 
j w a y 

Since x la the distance along the fceaa: x m a ( Ckx) 



Let i 4 * the atati&k nunbor at bearing 1 

S b • the atation nunoer at bearing 2 

& ■ the C at bearing 1 
a a 

& b • the S g at bearing 2 
Since the deflection y • at bearing 1 and bearing 2 then 
& a - C g S a ( A x) + Cy and © b • C s « b ( Ax) ♦ Cy 

126 
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Solving these two equations for C g and Cy and substituting into the equation 
for y then 



y t " S s * < § b^a * 8 « *b> < & b m *a> a 
I«b • «a > 

This ia the sane aa 



(•b • »a> 



^ - K ~ K ~ <^b lia> < »-.»«> 

Hence y it the difference in value between the value of S g and the straight line 
&^ which goea through each bearing aa shown in Figure 4. 



s 



line - *a+ " b ;^ < x-x«) 
(*b * x a) 




bearing 2 



runmi 4 



it u ii i j ii \i i 



Use of a Computer 



Beeeuae of the tedious labor involved in these general shaft calculations, the 
use of a computer ia neeeaaery. The input data for the computer computet lena 
differs slightly from that described in the example, which illustrated the method. 

The input data to the computer ia aa follows t 

1. Identification 

2. A x. 

3. Bearing location station numbers. 

4. External loada without maea and location a tat ion numbers. 

5. Horsepower transmitted with start and end station numbers. 

6. rpm 

7. Modulue of elasticity. fmgm 12 



In addition either of the following; 

3. Loads excluding shaft weight 
and external loada withour mass. 

9* Shaft O.D. with atart and and 
station number a. 

10. Shaft I.D. vlth start and and 
station numbers. 



8* Loada excluding external loada 
without mess. 

IV C § which la 1/2 shaft 0.9. 



J 



Values such aa 5, 9, and 10 which nay ba tha seme for a number of stations sra 
oaalar to pnt into tha computer aa a alnglo valna fron start station number to 
and a tat ion niaabor. External loada, 4, ara put into tha computer asperate ly so 
that thay nay ha used for calculations eiaillar to Fart 1 of tha example and 
excluded from calculatlona alnillar to Fart 2 of tha example. Xf tha shaft has 
an unnaual aoction so that c, and r cannot ha calculated from shaft 0.9 and 
I.D. than 8' 9; 10» and 11* ara tha Input data of that station instoad of 8, 9 
and 10. 

Tha calculations dona by tha computer ara similar to those in tha example except 
<<l $ (4ua to shaft weight), */^g t » <*/ Ip t » t # and IZ^ are calculated from £n»ut 

data 5, 6, 7, 8 and 9. the computer uses either input data 4 ♦ 8 ♦ (due to 

shaft wt) or input data 4 ♦ 8* as A § in calculations similar to Port 1. Zn 

calculations similar to Fart 2, the computer uees input data 8 ♦ J $ (due to shaft 

wt) or input data 8', with the loads between hearings reversed ,ae JL % , 

The principal advantages of these general shaft caleulatlona ara Hated belovt 

1. The calculatlona apply to any two bearing shaft with each bearing 
located anywhere on the shaft with any variable cross sections 
and loading and with any length. 



2. Boundary conditions are evaluated in tha calculatlona by having 

the first load at station 1 and using zero def lection at each 
bearing. 

3. Designs may ba optimised by calculating a number of variationa by 

computer and selecting the best. 

4. Variations in standard design, such aa a number of belt pulls, may 

be checked by making a number of calculatlona by computer. 

The accuracy of tha calculatlona are affected by the following) 

1. Having the load constant over the increment of length (Ax), 

2. Saving the ahaft cross section moments of inertia constant over 

the increment of length (Ax). 

3. Assuming the static deflection curve la the dynamic deflection In 

calculating critical speed. 
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4. Concentrating the lead in the increment of length In calculating 

critical speed. 

5. Uaing Ray la ighr. 1 a method to calculate critical speed. 

The error due to I and 2 preceding may be evaluated aa indicated on Page 8. 

The error due to 4 nay be minimised by taking small increments of length (An). 

The errors due to 3 and 5 nay be reduced by recalculating the critical apead 
uaing the inertia forces aa loads (4) 

Then calculating a second critical speed value 

Subscript 1 refers to values in the first calculation and subscript 1 refers to 
values in the second calculation. 

Koto that W §1 Is need in the calculation of » c 2 since the kinetic energy is a 
function of ease while potential energy is a function of the inertia force* 

QPm T v c tfT« frWM Iff « 

There are other factors affecting critical speed seek as shear deflection, (3) 
end (5), gyroscopic effects, (2), (3), end (6), bearing length (3), and bearing 
flexibility (3). Bearing aupport and oil file flexibility nay be accounted for 
by adding the bearing deflection to the y s in the calculation of critical apeed. 
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TABLE I - EXAMPLE CALCULATION - PART 1 
(Including external loads without mass and vlth normal direction of loads) 
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x » 4 


inches 


- Bearing No 1 at e • 


4 - Bearing No 2 at 
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/ JU 


26300 


0.08 


1120 
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U. AO 
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! 
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0.08 


1120 
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-8 
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1128 
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! 

1 
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-4 


-4 


i 
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> 
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i 
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448 




i 
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+8 
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! 

r. 
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! 

•67 2 1 

! 
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A / Of 
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8.5 i 






i 
] 


■ 


64 
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t 


; 

■ 

i 


! 


i 

! 
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-2 


♦2 


i 

+10 ] 
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1 

.) 

i 
i 


1 

-4/0 : 


J 
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O.'IU j 

i 


?i | 
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1 i 

1 




56 


-112 1 




. f 
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} 
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10.0 
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4-28 


| +84 ! 

i 
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i 
1 

j 


i 
j 


-28 


0.40 j 

i 


»l 
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i 
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i ; 


! 

J 


| I 




i 

i 


> 

1 
t 


f 

i 








TOTAL 


-252! -1384 


\ -376! 


! I 


T 

j. 




i 


— f 

1 

—4.. 


1 




1 





r x m + 346 #/in. r 2 - - 94 #/in. 
« 1384 lbs R 2 « -376 lbs. 
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TABLE 2 - EXAMPLE CALCULATION - FA&T 2 
(Excluding external loads without aass and with load* reverted between bearings) 
&x « 4 Inches - Bearing So. 1 at s ■ 4 - Bearing No-. 2 at s » £ 
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# 
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I 


-70 
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-4 
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CRITICAL SPEED* STRESS* AND BEARING REACTION CALCULATIONS 



C FOR A GENERAL SHAFT USING NUMERICAL INTEGRATION 

C ********** 

C PROGRAM DOES NOT INCLUDE A TRACE ROUTINE 



C NO SWITCH SETTINGS 



C ***** PART A-l ***** 

C **** DIMENSION STATEMENT 

DIMENSION A(102) »b(102) »C(102) 
C **** INITIALIZE ARRAYS AND SUMS LOCATIONS 
50 DO 51 K»l»102 



A{K)=0.0 
6(K)=0.Q 
51 C(K)=0.0 



SUMA=0«0 

suMa=o»o 



SUMC*0.0 



** READ IN EXTERNAL LOADS 



60 READ t I0*Wl»Il»W2>I2»W3tl3»W4»I4 
K*I1+1 



A(K)=W1 

K=I2+1 

A(K)=W2 

K«I3+1 

A(K)=W3 

K*I4+1 

A(K)=W4 

IF ( I0)60»60»100 

**** READ NO. 1 INPUT DATA 



100 READ »ISIZE,ICAL,ITYP»M,S1.S2»DX*HP,RPM 



PUNCH»ISIZEt : ICAL»ITYP»M»Sl»S2.HP.RPM,DX /33 



C **** CALCULATE TORQUE. *T.. ♦THE TURNING MOMENT ' 
T=63025/*HP/RPM 

M=M+1 ._ 

Sl=Sl+l. ........ _ 

S2=S2+1. ^ _ __ _ _ ^ _ _ 

C **** READ =2 INPUT DATA 

110 READ ' >0D > X I D » W * X I X , x IP , J ,L _____ ,._..,.,„,„ 

J=J+1 

L=L+1 

A(1)»0.0 ., _ _ 

IF (OD) 114*114»115 
C **** TO DET.IF IX»W ARE GIVEN 

115 IF (W)lll, 111*112 ^ _ _ _ _ _ „.„ 

C **** CALCULATE SHAFT WEIGHT W=PI/4(OD**2-ID**2)*0.285*DX 

111 W*0.7853982*(0D**2-XID**2 )*0«285*DX 

112 IF (XIX)113, 113*120 _ _ ............._.,.„ 

C **** CALCULATE THE MOMENT OF INERTIA OF THE CROSS SECTIONAL AREA » IX 
C WITH RESPECT TO THE RADIAL AXIS — I X= ( 0D**4-ID**4) *P I /64 

113 XIX=<OD**4-XID**4)*0. 049087385 

C **** CALCULATE THE MOMENT OF INERTIA OF THE CROSS SECTIONAL AREA » IP 
C WITH RESPECT TO THE LONGITUDINAL AXIS — IX=IY»SO IP»IXMY*2*IX 

XIP=2.*XIX _ 
C **** CALCULATE C/IX(S)»THE RADIUS TO THE OUTER SHAFT SURFACE DIVIDED 
C BY THE MOMENT OF INERTIA OF THE CROSS SECTIONAL AREA W.R.TO 

C RADIAL AXIS,IN.**-3 

120 CDIX=(OD/2. J/XIX 

C **** CALCULATE C/IP(S)»THE RADIUS TO THE OTHER SHAFT SURFACE DIVIDED 
C BY THE MOMENT OF INERTIA OF^THE CROSS SECTIONAL AREA W.R. TO 

C THE LONGITUDINAL AXIS»IN**-3 

CDIPMOD/2. )/XIP 

GO TO 116 



C **** CALCULATE LS ( A < K ) +LS ) »LOAD PER INCREMENT INCLUDING EXTERNAL ^ 
C LOADS WITHOUT MASS IN NORMAL DIRECTION. CONCENTRATED LOADS MUST 



C BE DISTRIBUTED OVER THE INCREMENT DELTA X. ACTUAL SHAFT AND 

C ROTOR LOADS ARE MORE OFTEN DISTRIBUTED THAN CONCENTRATED 

C **** DIRECTION OF LOAD » DOWN * - 
114 CJIX=G.O 

CDIP=0.0 
116 DO 130 I=J,L 

PJNCH»CDIX»OD»CDIP,XID 

A( I )=-(A( I )+w ) 

C **** CALCULATE M2 ( B ( J ) ) "MOMENTS ABOUT BEARING 2 STATION LINE 

C M2=LS*((S OF BEARING 2)-S) WHERE S=STATION NO. SO THAT M2 IS 

C CALCULATED AT EACH STATION LINE 

S = I 

B ( I ) sA ( I )*<S2-S) 

C **** CALCULATE Ml ( C ( J )) "MOMENTS ABOUT BEARING 1 STATION LINE 

C M1=LS*((S OF BEARING 1)-S) WHERE S=STATION NO, SO THAT Ml IS 

C CALCULATED AT EACH STATION LINE 

C ( I ) =A ( I ) * ( Sl-S ) 

SUMC=SUMC+C ( I ) 

SUMB=SUMB+B( I ) 
130 SUMAsSUMA+A ( I ) 

IF (L-M) 110»140>140 
140 I=S1 

A ( I ) »A ( I ) +SUMB/ ( S1-S2 ) 

SUMA=SUMA+SUMB/ (S1-S2 ) 

I=S2 

A ( I ) »A( I ) +SUMC/ ( S2-S1 ) 
SUMA«SUMA+SUMC/(S2~S1 ) 

C **** PUNCH OUTPUT DATA S»W»M2»M1 

DO 15C I*1>M 
S1«I-1 

150 P UNCH »S1»A(I)*B(I)»C(I) 
C **** PUNCH TOTALS OF M2»M1»W 
PUNCH 9 T » SUMd 9 SUMC > SUMA 



GO TO 50 
END 

C ***** PART A-2 **** 

C **** DIMENSION STATEMENT 

DIMENSION A<102) »P(102) *X(102) 
C **** INITIALIZE ARRAYS AND SUMS LOCATIONS 

100 DO 101 K=l >102 
P(K)=0.0 

101 X(K)=0.G 
SUMV=0.0 
SUMD=0.0 
SUME=0.0 

READ >ISIZE»ICAL»ITYP.M 
READ »Sl»S2tHP»RPM 

READ »DX 

PUNCH » ISIZE> ICAL> I TYP »M*S1 »S2 »HP tRPM.DX 
C **** STORE C/IX IN X(I) ARRAY 
C **** STORE C/IP IN P(I) ARRAY 

M = M+1 

DO 102 I=1»M 

READ ,CDIX»OD»CDIP »XID 
X( I )=CDIX 

102 P(I)=CDIP 

DO 103 1 = 1 »M 

103 READ » S > A ( I ) >Z»ZU 
READ »T,SUMB,SUMC»SUMA 
S = 0»0 

V1=0.0 
DX»0.0 

DO 120 1 = 1 »M 

C **** CALCULATE SHEAR A STATION (S+l/2) V ( S+l/2 ) = V ( S-l/2 ) +WS*DX 

V2=V1+A ( I ) *DX 

C **** CALCULATE bENDING MOMENT AT STATION (S+l/2) 



c 



M ( S+l/2 ) sM ( S-l/2 )+V ( 5-1/2 ) *DX+1/2*WS*DX**2 



D2*=D1+0.5*DX*( V1 + V2 ) 



C 



**** CALCULATE BENDING MOMENT AT STATION (S) M(S)=E(I) 



C 



-M(S)=M( S-l/2 )+l/2*V( S-l/2 ) *DX+ 1 / b*WS*UX**2 



El=Dl+0. 12i?*DX*( 3.*V1+V2 ) 

C **** CALCULATE bENDING TtNSILE STRESS>PSI S6S=MS*C/IX 

SB=E1*X( I ) 

IF (SB) 104»105»105 

104 SB=-SB 

C **** CALCULATE SHmFT TORSIONAL SHEAR STRESS»PSI STS=T*C/IP 

105 ST=T*P ( I ) 

C **** CALCULATE THE SHAFT MAX. COMdINED SHEAR STRESS, PSI 

C SCS=< ST**2+(SbS/2 )**2 )**.5 

SC=( ST**2+Q.2 5*SB**2 )**0»5 
C **** PUNCH OUTPUT DATA V ( S+l/2) »M ( S+ 1/2 ) »MS»C/IX»SBS»T,C/IP»STS»SC 

S=I-1 

IF ( X ( I ) ) 110»110»111 

110 Q=0. 

111 PUNCH »S»£1 »X( I ) »Sb »Q »P ( I ) »ST » SC 
Q=T 

S=S+.5 

PUNCH»S» V2 ,D2 
V1=V2 
120 D1=U2 

GO TO 100 

END 

C REARRANGE THE OUTPUT DATA OF 62-0015Al,2 INTO TABLE 
DIMENSION A(102) »B(102) »C(102) » T ( 1U2 ) » V ( 102 ) »D(102) 

100 READ, ISIZE,ICAL» ITYP,M 
READ,S1,S2»HP,RPM 
READ , DX 

PUNCH»ISIZE, 1CAL, I T YP »M ♦ S 1 > S2 » HP , RPM ,DX 
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C **** REAu»C/lX»OD»C/IP»Iu 

DO 110 I=ltK 
110 READ,R1 »R2»R3 »R4 
C **** READ S»WS»M2>M1 

DO 120 I = 1»K 

READ » S >R1,R2,R3 

PUNCH»S »R1»R2,R3 

S=S+.5 
120 PUNCH»S 

S = S-S 

REA0tR4»SUM6»SUMC»SUMA 

PUNChWS »SUMA»SUMB»SuMC 

R1=SUMB/ (S1-S2 ) 

R2=SUMC/ (S2-S1 ) 

READ, I SIZE* ICAL» I TYP»M 

READ,S1»S2,HP,RPM 

READtDX 

C **** A( I )=C/IP»B( I )=STS»C( I )*SCS 
DO 130 1 = 1 >K 
READ,S1»SM»CDIX»SB 
READ»T ( I ) ,A< I ) >B( I ) ,C( I ) 
READ,S»V( I ) , D ( I ) 

130 PUNCH»S1 tSMtCDIXtSB 
DO 131 I=1»K 

S= I . 

S=S-. 5 

131 PUNCH»S»V< I ) »D( I ) 

s=s-s 

PUNCH»S»R1»R2 
DO 140 1 = 1 »K 
140 PUNCH »T(I)»A(I)»B(I)»C(I) 

GO TO 100 

END ^£ 



4 q ***** PART B-l ***** 

C **** DIMENSION STATEMENT 

DIMENSION A(102) »X(lu2) »Y(102) 
C **** INITIALIZE WORING AREAS 

90 DO 91 I«l»102 
A(I)=0.0 

X( I ) =0.0 

91 Y(I)=0.0 
V1=0.G 
D1=0.0 
F1=0.0 
G1=0.0 
SUMB=G. 
SUMC=0. 

C **** READ INPUT DATA NO. 1 

READ » ISIZtt ICAL» I TYP»M*S1»S2»HP»RPM»DX 

PUNCH t IS I ZE» ICAL, I TYP »M,S1 ,S2 »HP ,RPM ,DX 
C **** READ C/IX (CulX) 

M=M+1 

S1=S1+1. 

S2=S2+1. 

DO 101 I = 1»M 

READ »X ( I ) >0l> »CP >CX 

IF ( X ( I ) ) 101»101»100 

100 X( I ) = (OD/2 . ) /X ( I ) 

101 CONTINUE 

C **** READ STATI0N»LS»M2*M1 
DO 105 .1 = 1 »M 
READ »S»A(I),B,C 
S=S+1. 

C **** CALCULATE LS=LOAD PER INCREMENT EXCLUDING EXTERNAL LOADS 

C WITHOUT MASS AND WITH LOADS REVERSED BETWEEN dEARINGS. 

C CONCENTRATED LOADS MUST Bt DISTRIgUTED OVtR THE INCREMENT DX. 



c 



ACTUAL ShAFT AN J ROTOR LOADS ARE MORE OFTEN D I STR I bUTED THAN 



C 



CONCENTRATED. 



IF CS-52) 102> 104*104 

102 IF (S-Sl ) 104»104»103 

103 A ( I )=-A( I ) 
B = -6 

C = -C 

104 SUMb=SUMb+6 
PUNCH»I»A(I) 

105 SuMC=SUMC+C 
READ ♦T»b»C»SA 
K = S1 

b=A(K)-6/(Sl-S2) 

A(K)=B+SUMB/(S1-S2) 

PUNCH»K»A(K) 

K=S2 

C=A(K)-C/(S2-S1) 
A(K) =C+SUMC/ ( S2-S1 ) 
PUNCH, K>A(K) 
DO 170 K=1»M 

C **** CALCULATE SHEAR AT STATION (S+l/2) V ( S+l /2 ) »V ( S-l/2 ) +WS*DX 

V2=V1+A (X ) *DX 

C **** CALCULATE bEMDING MOMENT AT STATION (S+l/2) 

c M(S+l/2)=M(S-l/2)+V(S-l/2)*DX+l/2*WS*DX**2 

D2=D 1+0 • b*0X* ( V1 + V2 ) 
C **** CALCULATE THE AREA UNDER THE (M/EIZ) CURVE AT STATION (S+l/2) — 
c — F (S+l/2 >=F< S-l/2 )+<2M(S-l/2+ 1/2V( S-l/2 )+M( S+l/2) )*DX/(3*EIX) 

IF ( X ( K ) ) 11 1 1 10 * 120 
110 F2=F1 

G2=Gl+Fl*DX 
Y(K)=G1+.5*DX*F1 

GO TO 130 J 




120 F2=F1+(2.*D1 + 0.5*V1 + D2-)*DX/ ( 3 .*X ( K ) *30 • E06 ) 



C **** CALCULATE. THc AREA UNDER THE F CURVE AT STATION (S+l/2) 

C G(S+l/2)=G(S-l/2)+DX*(F(S-l/2)+( 12M(S-l/2 ) +DX ( 3V ( S-l /2 ) 

C +V(S+l/2) ) )*DX/24EI ) 

G2=G1+DX* <Fl+< 12.*D1+DX* ( 3.*V1+V2 ) ) *DX / ( 24 .*X ( K ) *30 • 0E06 ) ) 

C **** CALCULATE THci AREA UNDER THE F CURVE AT STATION S 

c G =rG(S-l/2)+l/2DX(Fl+(M(S-l/2) + l/48*DX(7V(S-l/2>+V(S+l/2) 

Y(K)=G1+.5*DX*(F1+(D1+.020833 33*DX*( 7.*V1+V2 ) )*DX/ < 4 . *X ( K ) *30 . E06 ) ) 

C ♦*** CALCULATE DEFLECTION AT STATION S GS1=G1»G=GS2 

130 PUNCH »V1»V2»D1»D2»F1»F2»G1»G2>Y(K.) 

V1=V2 

D1=D2 

F1=F2 

G1»G2 
170 CONTINUE 

K = S1 

G1«Y(K) 

A(K)=B 

K=S2 

G2=Y(K) 

A(K)=C 

C **** CALCULATE ABS. SUM OF WS*YS — WS*WS*DX 
SUMP=0«0 
SUMK=0«0 
DO 240 K=1»M 
S = K 

Y(K) =Y(K)-G1-(G2-G1 )*(S-S1 ) / ( S2-S1 ) 

SUMP=SUMP+ ( A ( K. ) *0X ) *Y ( K ) 

SUMIOSUMK + ( A ( K ) *DX ) * Y ( K. ) *Y ( K ) 

PUNCH » SUMP ♦ SUMK » Y ( K ) 

IF (SUMP)210»220»220 
210 SUMP=-SUMP 
220 IF (SUMK)23C»240»240 
230 SUMK=-SUMK 



240 CONTINUE 

C **** CALCULATE THE FIRST CRITICAL SPEED » RPM 
C NC*187.7*(SUM OF WY/SUM OF WY**2)**.5 

XNC=187.7*(SUMP/SUMK)**.5 
C **** PUNCH OUTPUT — YS AND NC 

R1=SUM6/ (S1-S2 ) 

R2=SUMC/ (S2-S1 ) 

PUNCH»R1>R2»XNC 

DO 250 1=1 »M 

K=I-1 
250 PUNCH»K»Y(I) 

GO TO 90 

END 
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INTRODUCTION 



An ever present problem in engineering is the need to represent the three dimensional array 
of data by a mathematical equation. Many physical quantities can be described only in 
three-space. Computer storage limitations prevent table look-up of a large number of these 
quantities, expecially if their range of values is extensive. M-151 uses nth order multiple 
interpolation together with extensive transformation of variables to describe most three 
dimensional functions. The method used has the advantage of separating the three dimensional 
characteristics of the array into components described in two dimensions, hence easy to de- 
fine by the user. Since many functions (exponentials, sinusoidals, etc.) are not best described 
by the polynomials used in interpolation, the first approximation to the function might be sadly 
deficient. Better successive approximations to the function can be easily made using the 
appropriate output option accompanied by transformation of variables. An error analysis 
assists the programmer in choosing the best fit of a function, and in the case of second and 
third order interpolation, the best fit is selected from a number of possible fits. 



ANALYSIS 



Let us represent a three dimensional array by x = f |Y(y), Z(zJ] , where Y(y) and Z(z) are 
polynomials in y and z, respectively. Let us limit the order of these polynomials to n and 
m so that 

Y(y) = a + a iy -f • • • +a n y n = £ a.y 



j = 



m 

Z(z) = b n + b,z + • • • + b z m = Y b.z' (2) 
1 m z_ _ i 

At the ith value of y, 
m 

z. = 2 b.. z i 
' i = o 11 



Define: 



x = 



n n m 



f( Y ,z) = £ V = £ £ yv (3) 

1 i=0 i=0 j=0 1 



m 



Then, a. = V b..z' (4) 
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Note that x is a smooth function of y and z with (n + l)(m + 1) coefficients and hence re- 
quiring the same number of known members of the array. In general, an (n + l)(m + 1) order 
matrix would have to be solved to obtain values for the unknown coefficients. 

Referring to equations (3) and (4), note that both a. and x vary in the same way with changes 

in z. Thus, by studying the behavior of x with changing z, the behavior of a. with changing 

z is also known. By holding y constant, a. is a function of z alone. If values of a. are 

known at each of (m + 1) values of z, the b.. values in equation (4) are defined by the 

appropriate matrix solution. If z is held constant at z. and if (n + 1) values for x = g(y, z«J 

are known, (n + 1) a., coefficients are determined. Define a restriction that allows easy 

determination of the a., and b.. coefficients: y takes on, at most, (n + 1) values; z takes on, 

'11' 

at most, (m + 1) values; and the value of x is defined for all allowed values of y and z. At 
particular z., 



which is a form that allows solution for (n + 1) values of a.., when the (n + 1) values of y 
are used. 1 



x.(y, z.) = a . + a ]jy 



n 



(5), 




l =a 01 +a ll>' + '" +a nl y 



n 



n 



x 





X m+1 a 0,m+l +a l,m+l y 



7 

/ 
/ 

/ 

/ 



/ 



/ 



/ 



/ 



/ 



z 



Z m+1 



FIGURE 1 
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Choose the (m + 1) values of z successively, then an array of a *( z j) values results: 
i = 0, 1,2, • • •, n; j = 0, 1, 2, • • •, m. Fit the (m + 1) values of a. by 



m 



a.(z) = b n . + b..z+- • -+b .z (6), 
i Oi li mi 

with the result that the b.. coefficients in equation (3) are evaluated. Equation (3) defines 

a function of y and z that passes through (n + l)(m + 1) data points. Furthermore, the function 
of y and z is smooth, and if planes parallel to the y or z axis and perpendicular to the yz 
plane are chosen, they intersect the x = g(y, z) surface in lines described by polynomials 
of order n and m, respectively. For practical purposes, n and m may be selected by choosing 
their values dependent on the study of the z = constant plane intersections and the y = 
constant plane intersections with the surface x = g(y, z). Note that experimental data are 
often determined most easily by holding all variables constant except one. This procedure 
guarantees compliance with the restrictions imposed regarding values of x at each combination 
of allowable y and z. 

The usefulness of the assumed fit given by equation (3) is limited when the surface to be fit is 
best described using exponential, trigonometric, fractional exponent and hyperbolic functions. 
Although it is true that most useful engineering functions can be described by infinite series 
approximated by the form given in equation (3), it would be better to use the appropriate 
functional transformation directly. For instance, assume that a surface is known to be of the 
form x = f(y) = |sin(7y + 3) J + a^yj . Let x^ = f(y^)andx2 = f^) an d define 

g(y) = sin(7y + 3), x'j = Xj/gfy), x£ = x^ 9 ^ and ' ' n 9 enera '/ x * = x /g(y)- Then fit 
x 1 = Oq + a^y by a first order polynomial using x^ and x^ as data points. 

Having determined the values of Oq and a^, x is defined by x'g(y). The expression for x is 
valid not only at y^ and y^i but over the entire region in which g(y) is a suitable transfor- 
mation. It should be noted that the data points x^ and ^ could have been fitted by 

x = a n + a.y without considering g(y). However, only the two data points and the points 
-1 

where (7y + 3) = sin (1) would be satisfied by x = Oq + a^y • 

Theoretically, provided the matrix solution is exact, the polynomial fit given by equation (3) 
should fit exactly the (n + l)(m + 1) primary data points. If the function to be described is of 

n m 

the form x = g(y) £ b..z'y', a study of the primary data points will not give 

i=0 (=0 M 

information concerning the suitability of g(y) as the appropriate transformation of variables, 

since all primary data points are fitted, independent of the form of g(y). Hence, secondary 

data points not used in the determination of the b.. coefficients must be studied to arrive at 

M 

g(y). M-151 allows secondary data points for either arbitrary y or z, but not both arbitrary. 
Either y or z must be an "allowed" value, in the sense that defines primary data points. The 
other coordinate, z or y respectively, may be arbitrary. 
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To help decrease the importance of the user in deciding on allowed values of y and z for 
primary data points, over-definition of the problem is allowed in the cases of second and 
third order fitting. As many as seven allowed values may be chosen and the program will 
use all combinations of seven points taken three and four at a time, respectively for second 
and third order fits. 



SURFACE FITTING TECHNIQUE 

1 . Select the dependent variable, x 

2. Plot families of curves 

a. y = constant plotted on the xy plane 

b. z = constant plotted on the xy plane 



n . m 

i 



3. Assume n and m for fitting y and z such that x = £ a.y and x = ])T b.z' 

i=0 ' j=0 1 

on the xz and xy planes, respectively. If there is uncertainty concerning the best 
value of n and m, all possible values should be used in separate cases, using the 
average error and maximum error features of the code to decide between cases. It 
is not crucial that the very best case be used, if transformation of variables can be 
used to force a best fit. For instance, x = sin(7y + 3) is fit best in certain regions 
2 

by x = A a, y w '^ nout transformation of variables, but with transformation of 
i = ' 

variables x/sin(7y + 3) = x 1 = 1 is fit exactly by a zero order fit: x' = a^ = 1. 

In the interest of developing a technique independent of curve fitting experience, it 

2 

is better to use x = y a.y , provided that the region of definition will remain 

i = 1 

small. However, any curve fitting experience at the user's disposal should be used., 

4. To effect more accurate fits, use the output option in which transformed data is out- 
put ready for further transformations. This option prints the transformed values of x 
(for example, x/g(y) is printed if the transformation is x' = x/g(y) and x is printed 
if there is no transformation), the fit of the transformed variable, the difference in- 
cluding sign between x' and its fit, and the ratio x'/(fit of x 1 ). Referring to figure 2 
x = f(y), the fit of x = g(y), and h(y) = x - (fit of x), 



g(y) 




FIGURE 2 



The transformation that effects an exact solution is f(y) = h(y) + g(y). Since h(y) 
can be approximated by a sinusoidal with appropriate amplitude and period, the 
sinusoidal transformation is made on all data points using the program to make the 
change: x 1 = x - h(y). Upon rerunning this case with the transformation included, 
the printout will include x 1 , the fit of x', x 1 - fit of x 1 , and the ratio x'/(fit of x 1 ). 
Once again, h'(y) = x' - fit of x' = f'(y) - g'(y). To determine our fit, x 1 = 
h'(y) + g'(y) and x = x 1 + h(y) = h(y) + h '(y) + g ! (y). Correspondingly, the ratio 
x/(fit of x) may be used by plotting this ratio and describing it just as h(y) was de- 
scribed. Code capability allows as many as ten transformations to be made. It is 
necessary, in order to achieve an improvement in the fit, that the transformation 
of variables be such that the original data points will not be seriously shifted by the 
transformation. As an example, consider the translation without distortion of the 
fit of the array to more closely approximate secondary data points in the array. If 
this translation is fit once again, the fit of the translated array will be parallel to the 
original fit with no improvement over the original fit in regard to secondary data points 
in the array. 



TRANSFORMATION OF VARIABLES 

The utility of transformation of variables was demonstrated in the section on surface fitting 
technique. Basically, there are two types of transformations available in this code: func- 
tional transformations that deal with the manipulation of complete generated functions of x, 
y and z; and transformations of individual variables x, y and z to form the functions used 
in functional transformation. Functional transformations available are listed in Table 1 
and transformation of variables are listed in Table 2. 
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TABLE 1 

FUNCTIONAL TRANSFORMATIONS 



NVT NMODE = 1 NMODE = 2 

1 None None 

2 f(x) f(x) 

3 xg(y) xh(z) 

4 x + g(y) x+h(z) 

5 x/g(y) x/h(z) 

6 x[g(y) + k(y)] x[h(z)+l(z)] 

7 x/[g(y) + k(y)] x/[h(z) + l(z)] 

8 xg(y)h(z) x/[g(y)h (z)] 

9 x[g(y) + h(z)] x/[g(y) + h(z)] 
10 xg(y)/h(z) xh(z)/g(y) 
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TABLE 2 

TRANSFORMATION OF VARIABLES 
Let v take on the value x, y or z. v shall be transformed. 



TRANSFORMATION 



INVERSE TRANSFORMATION 



NT = 1 
None 



None 



NT = 2 

. cv + d 
a + be 



In [fei 



I cv + d\ 
be ) - a 



-d /c 



NT = 3 

i c 
a + bv 



"(a + bv C ) - a ] 
. b J 



NT = 4 

a + b sin (cv + d) 



sin" 1 | + b S '" ^ CV + ^ " a ^ |- d 



NT = 5 

4-U dv 
a + be 



(a + bc^ V ) - a 



d In c 



NT = 6 

a + b In (cv + d) 



exp 



[a + b In (cv + d)] - a 1 
b 



* 



(Table 2, Continued) 



NT = 7 



a + b In (v + d) 
c 



NT = 8 
a + b In 

NT = 9 
b 



v + d 



a + 



(v + c) Q 

NT = 10 

a + b sinh(cv + d) 



NT = 11 

a + b cos(cv + d) 

NT = 12 

a + b tan ^ (cv + d) 

NT = 13 

a + b cosh(cv + d) 

with the restriction that (cv + d) > 1 



Ir l \ 

a + b In (v + d)J - a 



'V 



- d 



[ a + b,n 7+t] ' a 



ce 





b 


b 










- a 


(v + c) d 





- d 



1/d 



- c 



1 + 



[a + b sinh(cv +. d)] - a 
b 



1/2 



[a + b sinh(cv + d)] - a 

+ b 




cos 



1 / [a + b cos(cv + d)].~ a 



- d 



tan (^Btan-V + d) 3 -a), d 



± In 



- 1 + 



[a + b cosh(cv + d)] q j 2 



1/2 



[a + b cosh(cv + d)| - a 
b 



- d 
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INPUT PARAMETERS 



AT, BT, CT, DT — These are the coefficients a, b, c, and d defined in Table 2. Referring 
to Table 1, these coefficients apply to each single function transformation, but apply 
only to g(y) if NMODE = 1 or h(z) if NMODE = 2, when a transformation is used 
that involves two distinct changes of variables. (An example is NVT = 6 and 
NMODE = 1 from Table 1, where g(y) and k(y) are distinct changes of variables 
having distinct characteristic coefficients.) 

AT2, BT2, CT2, DT2 — These are the coefficients a, b, c, and d defined in Table 2. Refer- 
ring to Table 2, these coefficients apply only when a two function transformation is 
used. When NMODE = 1, only k(y) or h(z) is described; and when NMODE = 2, 
only I (z) or g(y) is described. 

DC — x values defined as primary (see analysis). Values of x are allowed only for (YC, ZC) 
coordinate pairs. To fill input cards, nonvalid values of DC are set equal to 1 .000E-9. 
The order in which DC values are placed on cards is defined by the order in which ZC 
values are read. YC is held constant for any DC input card. 

DY — x values defined as secondary (see analysis). Arbitrary values of z are allowed for 

each YC value. To fill input cards, nonvalid values of DY are set equal to 1 .000E-9„ 
The order in which DY values are placed on cards is defined by the order in which Z 
values are read. YC is held constant for any DY card and packets of DY and Z cards 
are read in the same order as YC values were read. 

DZ — x values defined as secondary (see analysis). Arbitrary values of y are allowed for 

each ZC value. To fill input cards, nonvalid values of DZ are set equal to 1 .000E-9. 
The order in which DZ values are placed on cards is defined by the order in which Y 
values are read. ZC is held constant for any DZ card and packets of DZ and Y cards 
are read in the same order as ZC values were read. 

KO — output selector: KO = 1 - only the b.. coefficients are printed; KO = 2 - b.. co- 

1 

efficients, x, y, z, fit of x, x - (fit of x), average per centage error, and maximum 
per centage error are printed; KO = 3 - same as KO = 2 with the addition of the 
transformed x = x 1 , fit of x', x 1 - (fit of x 1 ), and x'/(fit of x') printed. Table 3 shows 
a fourth order function fitted with a first order approximation using KO = 1,2, and 3. 

MO — order of polynomial Z defined in equation (2), where MO = m. MO may be as 

large as ten, but calculational error is such that it is suggested that MO be kept less 
than seven . 

MC — (MO + 1) except for second and third order fits where MC may be as large as seven. 
There are MC values of ZC. MC < 12 

NC — (NO + 1) except for second and third order fits where NC may be as large as seven. 
There are NC values of YC. NC < 12 



IS3 
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NC2 — the number of secondary values of x to be read in for each value of YC. There are 
NC2 arbitrary z values. NC2 is < 18. 

NC3 — the number of secondary values of x to be read in for each value of ZC. There are 
NC3 arbitrary y values. NC3 is < 18. 

NE — operator that governs the type of error analysis usetdhto distinguish between possible 

fits when there are multiple fits allowed. Consequently, NE is useful only for second 
and third order fits where over-definition is possible (see analysis). 

NE TYPE FIT 

1 numerical average 

2 product 

3 numerical average X product 

4 averages/maximum error 

5 product with minimum error stripped out 

6 average with maximum stripped out 

7 product with maximum stripped out 

8 product with maximum and minimum stripped out 

NMODE — an operator, together with NVT, that is used to define transformations in Table 1 . 

NO — order of polynomial Y defined in equation (1), where NO = n. NO may be as large 
as ten, but calculational error is such that it is suggested that NO be kept less than 

seven . 

NPASS — number of transformations allowed sequencial ly . Referring to the section on tech- 
niques, NPASS takes on the value 2 if x" transformations are made. 

NT — defined in Table 2 and treated similar to coefficients AT, BT, CT and DT. 

NT2 - defined in the same way as NT but treated similar to coefficients AT2, BT2, CT2 and 
DT2. 

NVT — an operator, together with NMODE, that is used to define transformations in Table 1 . 

Y — an arbitrary value of y which, together with a ZC, defines a secondary point, x. 

YC — the value of y used to define primary data points in the DC array and secondary data 
points in the DY array. 
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Z — an arbitrary value of z which, together with a YC, defines a secondary point, x. 

ZC — the value of z used to define primary data points in the DC array and secondary data 
points in the DZ array. 

INPUT TECHNIQUE 

Cards 3 and 4 below are read in repetitively in packages consisting of one card 3 followed by 
one card 4 provided MC ^ 6 or by two card 4s when MC > 6. NC of these packages are read 
in. Card 5 is repeated MC times. Cards 6, 7 and 8 form a package that is read in NC times. 
If NC2 = 7 cards seven and eight are not included in the package. If < NC2 < 6, fust one 
card 7 and card 8 follow NC2. If 6 < NC2 ^ 12, the sequence is 6, 7, 8, 7, 8. If 12 < NC2 < 18, 
the sequence is 6, 7, 8, 7, 8, 7, 8. Cards 9, 10 and 11 form a package that is read in MC times. 
This package is handled the same as the 6, 7 and 8 package. Cards 13, 14, 15 and 16 form a pack- 
age that is read in NPASS times. If NVT (I) < 6, cards 14 and 15 are omitted. 



Card Description Format 

1 55 symbol comments card 55H 

2 NO, MO, NC, MC, NE, KO 614 

3 YC(I) E8.4 

4 DC(I,J), DC(I,J + 1), DC(l,J + 2), DC(l,J + 3), 6E8.4 
DC(I, J + 4), DC(I, J + 5) 

5 ZC(J) E8.4 

6 NC2 14 

7 DY(I,J), DY(I,J+ 1), DY(l,J + 2), DY(I,J + 3), 6E8.4 
DY(I, J + 4), DY(I, J + 5) 

8 Z(I,J), Z(I,J+ 1), Z(l,J + 2), Z(I,J + 3), Z(l,J + 4), 6E3.4 
Z(l, J + 5) 

9 NC3 14 

10 DZ(J, I), DZ(J +1,1), DZ(J + 2, 1), DZ(J + 3, 1), 6E3.4 
DZ(J + 4,I), DZ(J + 5,I) 

11 Y(J,I), Y(J+ 1,1), Y(J + 2,I), Y(J+3,I), Y(J + 4,I), 6E8.4 
Y(J + 5,I) 

12 NPASS 14 

13 NVT(I), NMODE(I), NT(I) 314 

14 NT2(I) 14 

15 AT2(I), BT2(I), CT2(I), DT2(I) 4E8.4 

16 AT(I), BT(I), CT(I), DT(I) 4E8.4 
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KO - I. 



XIHR66 OIN6NSIONAI SURFACE FIT CODE , M-151 
N-151 DEMONSTRATION PROKlEMt OROERMXl 

CcEFF ICItNTS OF THE POWER SERIFS I N 1 



fll 1 



1 ) 

2 1 



20.95776-01 

10.9032fc-04 



I 7. d064t-04 
-65. >496t-08 



XTH'IH DIMENSIONAL SUrfFAtfc FIT CuOF, 
>• - 1 j i OlMONSTRAT IUN PK031E*, OROtR*Ul 

CufcFF|CU;US OF THE PuKER StRlFS IN I 



HI I, I) * -20. >577 L -01 

BI 1, 2) = U. 80326-04 



H ( 2, 1) * 1 7.H0646-04 

e ( 2 , 2 I = -6 j. 94 961 -08 



TArtUlAIt.O 
X 



CAK.ULITED 



niFFtamtt 



50.000Ce-00 27.0000F*0;> <}?.ll00fc-02 62.1100F-02 



.OOOOt.- )> 

50.00006*02 32.3000E-01 32.2J9JF-01 10.0000E-OH 
lO.OOOOtO? 27.0OO0t*02 82.11006-02 82. 10.0000i-08 
17.8?00fc-0l 17.8899E-01 20.0Q00L-08 
>0.9100t-02 11. 3i30E-Ot-?2.62CKf:- 



50.00006*02 
50.00006-00 30.00001*02 



33.0000E*02 10.1*001-01 14 . 4950E -0 1 -4 3 . 55031 - (V 



37.00006*02 11.>OOOE-Ol 1 8. 6844E -01 -67 . 8447E-02 
41.00006*02 14.7l00t-0l 22. 8 7 J8E-0 1-8 1 . 6 3866-02 
47.00006*02 23.80O0E-OI 29. t 57*6-01-53.57956-02 
49. 00006*02 29.09006-01 31. 2526C-01-21 .62656-02 
10.00006*02 30.00006*02 •«>. 54006-02 94. 7 34 76 -02-51 . 947eE-03 
33.0000E*02 97.2700E-02 10. 73596-01-10. 0895E-02 
37.00006*02 10.90006-01 1 2. 4 1926-01- 15. 1926E-02 
41. 00006*02 1?. 37006-01 14. 1 0256-91-17. 3256E-02 
47. 00006*02 I5.5600t-0l 16.62756-01-10.67526-02 
49.00006*02 1 7.04006-01 17.46916-01-42.91726-0 3 



AVfRAGE P6RC6NTAGI: EMW * . 122214«66*02 
MAXIMUM PERCENTAGE ERROR ■ .969107706*02 



XTHKEE OINCNSIONAl SUKFACE FIT CODE, M-151 
M-151 DEMONSTRATION PROELCH, UROfcRMXl 

COEFFICIENTS OF THE PUM6R SLRIES IN I 



HI 1, 1) 
BI I i 2) 



BI 2, II 
Hi 2t 2> 



-20. 9577t-0l 
10.8032t-04 



1 7.«064c-04 
-6>. )4 96t-06 



THli OPTION OUTPUTS THE Fl/llO* I 04 AlTtP.NATF L1NF^>-- 



ORIGINAL OA T A 



I A 'tl)l AT! l) LALCul A T t 



f IFF, *CF 



TMANSFOKMFO OAT A 

TAFUJLAtflJ CALCULATtO 

XT XC Xf-Xt 



50. OOOul -00 
82. 11 001. -02 

32. 30006-01 
10.00001*02 
82. 11001-02 

17.89006-01 
50. OOOOt -00 



I0.0000t*02 



27. 
82. 
50. 
32, 
27. 
62. 
50. 
17. 
30. 
33. 
37. 
41. 
47. 
49. 
30. 
33. 
37. 
41. 
47. 
49. 



00006*02 
I 100E-02 
0000E*02 
299*6-01 
0000F*02 
10996-02 
0000t*02 
88996-01 
00006*02 
00006*02 
00006*02 
0000t*02 
00006*02 
00006*02 
00006*02 
00006*02 
00006*02 
00006*02 
00006*02 
00006*02 



fl?. UOOt-02 
.0000E.-99 
32. 30006-01 
t0.0000t-08 
82. UOOf-02 
10.OOOOE-08 
17.8900E-01 
20.0000E-08 
90.9100E-02 
10. I400t-01 
11.90006-01 
14. 7100F-01 
23.80006-01 
29.0900E-01 
89.54006-02 
97.27006-02 
10.90006-01 
12.. 37006-01 
I5.5600t-0l 
17.04006-01 



82.1 1006- 
10.OOOOF- 
32.299 Jt- 
99.999 it- 
62.109 7E- 
99.99996- 
17.869)6- 
99.999JE- 
11.35306- 
14.4950F- 
18.A8446- 
22.87386- 
29.15796- 
31.25266- 
?4. 73476- 
10.7359E- 
12.4192E- 
14.10256- 
16.62756- 
17.46916- 



02 
01 

01 11. 

02 

02 10. 

02 

01 20. 
02 

0I-i2. 
01-43. 
01-67. 
01-81 . 
01-03. 

01- 21. 

02- 51, 
01-10, 
01-15, 
01-17, 
01-10, 
01-42, 



0000t-08 

00006-08 

6204t-02 
55CH6-02 
fc-02 
618ot-0» 
5 7")jE-02 
626^E-02 
9478E-03 
0895C-02 
1926E-02 
3256E-02 
6752E-02 
-J172E-0) 



AVERAGE PERCENTAGE ERROR « .122214866*02 
MAXIMUM PERCENTAGE ERROR - .363107706*0? 



TABLE 3 



I ' ! 
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SAMPLE PROBLEM 
1 

A practical problem involving the fit of hydrogen thermal conductivity as a function of 
temperature and pressure was met with success using M-151 . The technique of regression 
analysis proved unwieldy in solving the same problem; lack of success being attributed to 
inexperience in picking suitable form for the fit plus the inherent limitations of regression 
analysis, such as inaccuracies resulting from the need to solve for large numbers of co- 
efficients simultaneously. Some of the data points used are presented in Table 4. Referring 
to the section on Surface Fitting Technique, thermal conductivity was made the indepen- 
dent variable, pressure = y and temperature - z. The families of curves were not plotted 
since it was decided to choose many combinations of n and m and study the resulting error 
analysis. Table 5 presents input used when n = m = 4. Step 4 in the technique section 
was not used extensively since good results were achieved with the 4x4 order fit. Because 
of the volume of output with the 4 x 4 case, only the error analysis will be presented to- 
gether with the errors from other cases in Table 6. Table 3 presents examples of output 
available . 

TABLE 4 

THERMAL CONDUCTIVITY OF HYDROGEN 
(X10 6 ) 



TEMPERATURE, °R 





2700 


3100 


3560 


3900 


4300 


4600 


4800 


5000 


50 


; .8211 


.9425 


1.094 


1.312 


1.687 


2.166 


2.626 


3.230 


100 


.8211 


.9352 


1.073 


1.262 


1.571 


1.947 


2.301 


2,760 


200 


.8211 


.9277 


1.041 ; 


1.212 


1.451 


1.721 


1.965 


2.273 


300 


.8211 


.9264 


1.047 


1.203 


1.430 


1.682 


1.907 


2,191 


400 


.8211 


.9250 


1.043 


1.194 


1 .409 


1.643 


1.850 


2,109 


500 


.8211 


.9237 


1.039 


1.184 


1.388 


1,604 


1.793 


2.026 


700 


.8211 


.9210 


1.031 


1.166 


1.345 


1.526 


1.678 


1.862 


1000 


.8211 


.9204 


1.029 


1.159 


1.327 


1.491 - 


1.627 


1,789 



151 
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TABLE 5 



M-15I DEMONSTRATION PROBLEM, 0RDER=4X4 



4 


4 5 


5 2 


3 








50. 














.8211 


1.094 


1.687 


2.626 


3.23 


1.000E- 


9 


100. 














.8211 
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TABLE 6 
ERROR ANALYSIS 



Order (n x m) 

1 x 1 
1 x2 
2x2 
2x3 
3x2 
3x3 
4x4 
5x5 
6x6 
7x7 



Average % Error Maximum % Error 



12.22 

7.755 

8.681 

2.645 

9.439 

2.653 

.04456 

.07827 

2.201 

51 .29 



36.31 

18.20 

18.82 

13.25 

21.38 

9.839 

.3445 

.7006 

38.67 
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MAXIMUM LIKELIHOOD RESOLUTION OF TWO 
MIXED NORMAL DISTRIBUTIONS 



Reimut Wette 

Biomathematics Unit, Division of Research, 
The University of Texas M.D. Anderson Hospital and 
Tumor Institute, Houston, Texas 



Abstract ; 

Samples exhibiting particular deviations from an assumed normal 
distribution can, in certain situations, be interpreted on the 
assumption that the parent distribution consists of two normally 
distributed subpoputations , with different means and/or variances, 
and mixed in a finite proportion. The statistical problem is, then, 
to estimate the respective parameters. In the general case, where 
no simplifying: assumptions (e.g. equality of variances) can be 
made, five parameters have to be estimated, viz. two means, two 
variances (or standard deviations), and the proportion of mixture. 
The author's estimation method of choice was the maximum likeli- 
hood scoring system. The derived estimator is of the iterative type, 
and modified insofar as the information matrix is estimated as the 
dyadic square of the score vector. Numerical execution of this 
iterative estimation procedure requires initial estimates of the 
five parameters; these are, at present, obtained manually from a 
normit graph of the data. A generalized distance from origin of the 
score vector (i.e. the quadratic form of the estimated variance- 
covariance matrix) is used as a criterion to exit from the iteration 
cycles. 

The procedure was programmed for the IBM 1620 (with 1311) in 
FORTRAN II-D. The main program consists of two parts: 1.- Routines 
for control transfer set-up, standard form data input or transfer 
to either of two optional non-standard data input subprograms. 
2,- The resolution procedure proper, routines for optional inter- 
mediate and other output, and transfer to two optional subprograms. 



The program provides, at present, the following I/O (punched cards) 
and test options (under parameter card control): 



200 variates or classes. 

(2) Grouping of ungrouped raw data; class inter- 
val and grouping range computed from initial 
estimates. 

(3) Grouping of fixed- interval valued ungrouped 
data. 

Output and test options : 

(1) Output of grouped data frequency distribution. 

(2) Output of expected frequency and cumulative 
frequency distribution together with observed 
distribution (grouped data only). 

(3) Chi-squared test for normality against s&ewness 



Input options 2 and 3 and output options 2 and 3 are subprograms 
called as LOCALS on a 40k machine. A data trailer card controls 
exit to MONCAL or return linkage to part 1 of the main program. 
Maximum core storage used is 35>604 cores (Monitor I, modif. 2, 
all-in-core subroutines), program package (2 main/link, 4 subprograms) 
disk storage is 435 sectors. Interation cycle running time is about 
2.7 seconds per point (variate or frequency class). 

Convergence in this iterative estimation procedure depends on 
the goodness of the initial estimates, on sample size and structure. 
Experience gathered so far indicates that convergence, if attainable, 
is comparatively fast (less than 10 iterations). Convergence could 
not be attained in a few instances of small and ill-behaved samples, ( 
of size around 20 and less, from which reliable initial estimdes 
could not be obtained. 



Data input : 



(1 ) Ungrouped or grouped data not exceeding 



(4) 



and kurtosis (cumulant test), with programmed 
bypass of resolution peocedure when not 
significant. 

Intermediate iteration output of estimated 
information and variance-covariance matrix, 
parameter estimates and corrections, and 
convergence criterion. 



The author gratefully acknowledges the extensive assistance of 
Mr. Lawrence E. ITewton, Jr., in programming and testing different 
versions of the method in the Computer Science Laboratory of the 
M.D. Anderson Hospital and Tumor Institute. 



Note : Methodological and programming details will be published 
elsewhere. 



EVALUATION OF TWO METHODS OF FINDING SIGNIFICANT 



CONTRIBUTORS IN MULTIPLE REGRESSION ANALYSIS 



1/ 



M a Jo Garber and Richard H. Hill 



With the advent of computers in recent years tremendous strides have been 
made in the ability to reduce to manageable numbers vast accumulations of data. 
In evaluation of experimental results the major emphasis is now on punching and 
proofing the observed values with the assurance that a generalized computer 
program is probably on the shelf ready for use. In the majority of cases this 
is true, but thoughts and concepts have expanded along with the increase in 
ability to do arithmetic quickly and accurately. We are, of course, now find- 
ing problems even large scale computers cannot solve in a reasonable period of 
time. 

One such problem is multiple regression involving many independent vari- 
ables* A number of years ago the first of approximately 60 such problems was 
handed to the senior author One phase of the experiment dealt with growth in 
a greenhouse of citrus seedlings in 102 nonfumigated old citrus soils. A total 
of 35 measurements was made on each. These included chemical, physical, and 
biological properties of the soil, plant composition, and relative growth. 
Locating statistically the best single contributor and the best set of 33 of the 
3H independent variables is relatively a trivial operation. The best single 

2 

contributor is the variable with the highest coefficient of determination (r ) e 

x i v 

The best set of 33 (or n=^l) contributors is easily found by entering all n vari» 
ables into the regression equation, and deleting the variable with the smallest 



1/ 
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2 

absolute partial correlation (|r — - i ) , Mr 6 Hill (then at the University 

of California, Los Angeles) wrote the MISLE program for the 709, 

The best set of two, the best set of n-2, and all best sets between these 
extremes present a very different problem. The best set of two will not neces- 
sarily contain the best single contributor s and the best set of n-2 could contain 
the variable not included in the best set of n-1. For a statistically efficient 
analysis all best sets should be calculated. But this presents a formidable 
problem. For example, the final regression equation for the above problem con- 
tained 10 independent variables. Finding the best set of 10 among 34 variables 
can be quite a chore, even for a computer. The number of sets to be evaluated 
is 34!/(10!2*H) or 131,128,140, Assuming that a computer can solve these at the 
rate of one per second, over 36,424 hours are required for finding the best set 
of 10. 

A number of questions immediately come to mind; 

(1) Is the best set of 10 significantly better than the one 

selected by this procedure? 

(2) How does the best set of 10 compare statistically with 

the best set of 11, or any other best set? 

(3) As the correlation coefficients are statistics, and not 

parameters, is the additional labor worth doing? 

After some soul searching it was decided that initially the following three- 
phase procedure would be used? 

(I) MISLE Program (presently UCRBL 0052 9 1620 Library File No. 6.0,37) 
All independent variables are entered into the regression equation, the 

inverse matrix cf sums of squares and cross products being calculated in the 

2 ■ • . 

process. All b^/c^ are calculated, where b^ is the partial regression coef- 
ficient of y on the ith independent variable, and c** is the major diagonal 



IU 



* 

3 
2 

element for that variable in the inverse matrix e As Fisher (1) has shown b^/c^ is 
the variance (equivalent to a sum of squares with one degree of freedom) which 

would drop into the error sum of squares of y if the variable were deleted. 

2 

The variable with the smallest b /c is deleted, and the remaining n-1 
independent variables are entered into the regression equation. This procedure 
is continued until only one variable remains. 

An example of such analysis is illustrated in Table 1 in terms of the squares 
of partial correlation coefficients. Note that, in general, the simple correlation 
coefficients, the squares of which are shown in the second line of the table, do 
not reflect in any way the magnitudes of the squares of the partial correlation 
coefficients as variables are deleted. 

The leftmost 10 variables were found to be statistically significant contribu- 
tors (0.05 level of probability), and it should be noted that four of these showed 
nonsignificant simple correlation with the dependent variable. The magnitude of 
R 2 was 0.5387. 

Variable 34 has the largest simple correlation. Yet, in combination with 
other variables its contribution is relatively small. On the other hand, Variable 31 
which was not significantly directly associated with the dependent variable is rela- 
tively a large contributor in the presence of other variables, 

(2) Stepwise Regression (UCRBL 0018, 1620 Library File No. 6.0.031, a 

modification of the BIMD 9 developed at UCLA by the Biostatistics 
group) 

Here the procedure is to begin with the variable most highly correlated with y. 
From the remaining variables the one with the highest (absolute) partial correlation 
is selected for entry (assuming that it meets the F criterion), and the magnitudes 
of the two contributions are evaluated. Either both are retained or one is deleted 



t 



4 

(a parameter card entered with the data contains the minimum F value required to 
enter a variable into the regression equation and the minimum F value required for 
retention of a variable). This process of entering, evaluating, and deleting vari- 
ables continues until the criteria can no longer be met. 

For the problem being considered here nine variables were found to be 

2 

statistically significant contributors (R = ,5205) when grouped. Eight of these 
variables were also selected by the MISLE program. 
(3) MISLE Runs 

(a) The eleven variables selected by either or both of the above 
programs were fed into the computer for another MISLE evaluation. Ten of the vari- 
ables were found to be significant contributors, and the first variable of the 
eleven to be deleted was the one unique to the Program 18 run. There was thus no 
significant gain. 

(b) It was then decided to make another run, this time including 
the eleven above and certain additional variables. Examination of the F values 
(- 1) led to the choice of the two variables last deleted by Program 52 and the 
succeeding three which would have been entered by Program 18. Two of these were 
common to both runs. 

Of the fourteen variables entered for the final run 10 were found to be sig- 
nificant. Seven of these were among the variables found significant by both 

programs, two common, but in the noncritical area, and the tenth being a signifi- 

2 

cant contributor found by Program 18. The magnitude of R was .5433. 

The 10 variables statistically significant in this run were the ones con- 
cluded to be real contributors (paper by Martin, Harding, and Garber (2)) 6 

2 

The four curves are shown in Figure I. The curve of multiple R calculated 
by Program 52 (labeled 521) is almost smooth, running essentially horizontally 



until approximately six variables have been deleted. The drop off becomes quite 

steep as additional variables are deleted* It should be noted that the last remain 

ing variable is not the one most highly correlated with y. The latter 9 Variable 34 

was dropped after the 26th run. 

2 

The curve of multiple R calculated by Program 18 begins by entering 
Variable 34 then rising in a somewhat oscillatory manner to essentially match 
the 521 curve after 9 variables have been entered* 

The short spur labeled 5211a is the graphic result of the MISLE run with 
11 variables. 

The short section of curve (5211b) exhibits a sharper break than either of 
the first two* 

The remaining two curves in the figure are from the 521 run* The curve of 
2 

"Cumulated b /c" shows the addition to the residual sum of squares (all variables 
initially entered into the regression equation) as variables are deleted* The 
curve of "Residual Mean Squares" falls off from its initial value 9 then rises,. 
Recalling that the residual mean square is the quotient of the residual sum of 
squares by the residual degrees of freedom it is seen in this example that in the 
initial phases of variable deletion the denominator is increasing relatively more 
rapidly than the numerator* Only after approximately half the number of variables 
has been deleted does the curve begin to rise, 

A summary of the evaluation of 55 experiments is given in Table 2* For the 
left side of the table the F criterion was set at 1*00* For the right side of the 
table the 0,05 level of probability was the criterion* 

In eight of the 55 experiments Program No„ 18 was the better of the two 
programs. In four experiments Program No* 52 gave better results* In one experi- 
ment each of the programs provided information both significant and unique. 




In one experiment the order of entering variables indicated the advisability 
of a third run. The outcome was the set of four variables selected by both programs* 

For the remaining 41 experiments both programs gave the same results in the 
initial runs* The third run was unnecessary. 

In summary, there is no clear cut indication of the superiority of one method 
over the other. For our purposes at Riverside we will continue to evaluate large 
problems by both the deletion and the stepwise methods. 

LITERATURE CITED 

1. Fisher, Sir R. A. (1948). Statistical Methods for Research Workers, 

10th Edition, Hafner Publishing Company. 

2, Martin, J. P., R„ B, Harding, and M, J» Garber. (1961). Relation of Soil 

Properties and Plant Composition to Growth of Citrus Seedlings in 100 
Nonfumigated and Fumigated Old Citrus Soils. Soil Science 91(5): 317-323. 



no 



o 



6184 (3 


1 ( '> 


32 


31 


1 1 


6 


8 


2 


29 


34 


24 


20 


27 


22 


16 


17 


23 


19 


9 


15 


5 4 


30 18 


28 


3? 


13 


7 26 


3 


21 12 


14 


10 


25 


]2< 2 >11 


1 


7 





5 


4 





16 





1 


12 


11 











4 


3 


10 


4 2 


3 1 


1 





6 


1 1 


1 





2 


1 


6 


3 


y 8W)i3 


10 





3 


2 


1 


4 


5 


5 


7 


4 


1 


1 


1 


3 


2 








1 


1 





1 


























6184 


8 


13 


10 





3 


3 


1 


4 


5 


6 


7 


4 


1 


1 


1 


■5 


2 








1 


1 





1 

























6184 


8 


13 


10 


3 


3 


8 


1 


4 


5 


6 


7 


4 


1 


i 


1 


3 


2 


2 


1 


1 


1 





1 


1 





















6184 


8 


13 


10 


4 


5 


8 


1 


5 


5 


6 


8 


4 


1 


1 


1 


3 


2 


2 


1 


1 


1 





1 


1 




















6182 


8 


14 


10 


7 


6 


8 


3 


5 


5 


6 


7 


4 


1 


1 


1 


3 


2 


2 


1 


1 


1 





1 


1 



















6179 


9 


14 


10 


7 


5 


8 


3 


5 


5 


6 


8 


4 


2 


1 


1 


3 


3 


2 


1 


1 


1 


1 


1 


1 


















6175 


9 


14 


10 


7 


5 


8 


3 


5 


5 


6 


8 


4 


2 


3 


2 


3 


2 


2 


1 


1 


1 


1 


1 


1 

















6171 


9 


14 


10 


7 


7 


8 


3 


5 


5 


6 


8 


4 


2 


3 


2 


3 


2 


2 


1 


1 


1 


1 


1 


1 

















6164 


9 


15 


10 


7 


7 


10 


3 


6 


5 


6 


9 


5 


2 


4 


2 


3 


2 


2 


2 


1 1 


1 1 


1 






















6146 


10 


15 


12 


9 


7 


9 


5 


6 


6 


6 


8 


4 


1 


4 


2 


3 


2 


2 


2 


1 1 


1 1 
























6131 


11 


15 


13 


8 


7 


9 


5 


6 


6 


6 


8 


5 


1 


3 


2 


3 


2 


2 


2 


1 1 


1 























6120 


11 


15 


13 


8 


7 


9 


5 


6 


7 


6 


8 


5 


1 


4 


2 


2 


2 


2 


2 


1 1 


1 






















6104 


11 


15 


13 


9 


8 


9 


5 


5 


8 


6 


8 


5 


1 


3 


2 


2 


2 


2 


2 


1 1 


1 






















6067 


10 


15 


13 


8 


8 


10 


5 


6 


8 


6 


7 


5 


2 


4 


3 


2 


2 


2 


2 


1 I 
























6035 


10 


15 


13 


8 


10 


9 


5 


6 


8 


7 


7 


5 


1 


4 


3 


2 


2 


1 


1 



























6018 


9 


15 


13 


8 


10 


9 


5 


6 


8 


6 


8 


5 


2 


4 


3 


2 


2 


1 


1 


























5967 


10 


17 


12 


7 


12 


8 


7 


7 


6 


5 


8 


4 


3 


3 


2 


2 


1 































5948 


10 


17 


11 


7 


12 


8 


9 


6 


6 


6 


8 


5 


3 


4 


2 


2 


2 






























5876 


15 


16 


12 


7 


11 


7 


8 


7 


8 


5 


7 


4 


2 


3 


3 


1 
































5817 


14 


15 


1 1 


8 


11 


8 


8 


8 


9 


5 


6 


3 


2 


3 


2 


































5714 


12 


18 


11 


8 


10 


9 


8 


7 


7 


5 


4 


3 


2 







































5693 


12 


18 


13 


7 


10 


9 


7 


7 


7 


5 


4 


3 


i 






































5644 


11 


17 


14 


6 


10 


8 


7 


6 


6 


4 


4 


3 








































5489 


14 


17 


14 


10 


10 


8 


8 


7 


6 


5 


2 






















LEGEND 


















5387 


13 


17 


14 


10 


11 


9 


8 


6 


5 


4 












































5171 


12 


17 


17 


11 


11 


9 


7 


6 


3 














1. 


Variable number 


(Variables were deleted from right 


to 


left) 




5006 


16 


24 


21 


12 


9 


7 


7 


6 
















































4682 


12 


25 


19 


9 


9 


6 


5 


















2. 


Simple 


r 2 


(Two decimal digits) 
















4426 


14 


22 


16 


7 


6 


4 


























o 


























4194 


28 


19 


17 


6 


5 






















3. 


Multiple R y 


.ijk- 


(Four decimal digits) 














3866 


26 


15 


13 


7 






















































3434 


26 


16 


12 


























4. 


Partial r y i.jk-- 


(Two decimal digits) 














2565 


16 


16 




























































1154 


12 































































Table 1. Computer results (34 independent variables, 102 sets of 
observations) . 



Table 2„ Summary of evaluation of 55 experiments. 
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Summary of 41 additional experiments: 

No. of independent variables: 6 to 32 
No. of degrees of freedom: 24 to 152 

No. of significant contributors selected in common by 0052 and 0018: to 5 
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IBM 1620 ELECTRONIC CIRCUIT ANALYSIS 



The Electronic Circuit Analysis system is a series 
of programs to help electronic engineers greatly 
extend their design capabilities through the use of 
the IBM 1620 Data Processing System. Previous 
computer experience is not required. 

Electronic circuit performance requirements 
vary between industries and products. In some 
instances, the design objective is extremely high 
reliability; in others, it is a balance between cost, 
reliability, and consumer acceptability. The design 
effort includes determining the tolerance required 
of each individual circuit component in order to 
obtain an overall circuit performance within certain 
tolerances. Since the relationship between component 
tolerance and circuit tolerance is not a direct one, 
the analysis of the possible combinations of com- 
ponents is a time-consuming and tedious task to 
perform manually. 

The IBM 1620 Electronic Circuit Analysis system 
provides the user with a high degree of confidence 
in circuit design, since many variations of circuit 
performance can be "tested". The system is user- 
oriented, providing means for easy communication 
and complete user control during the solution of 
the problem. 

To use the system, the engineer starts with a 
schematic representing the equivalent circuit. He 
then numbers the nodes and branches of this net- 
work. A special input language is used to describe 
the network, its physical characteristics, and 
various user options such as type of output, type of 
analysis to be performed and circuit parameters to 
be varied. The input statements are then punched 
into cards in free form for entry into the 1620. 
The system interprets the input statements, de- 
termines the network topology, generates the 
network equations and performs the desired analysis. 
During the course of solution, the engineer can 
easily intervene at the console and modify system 



parameters, output desired, or many of the input 
language instructions by merely typing modified 
language statements on the console typewriter. Thus, 
the system allows either intervention for manual 
modification or completely automatic modification 
and solution. In this manner, the engineers can 
simulate their normal mode of design, which is 
experimental at the outset and iterative on a fixed 
design in the final stages. 

The system consists of four sections: the monitor, 
the AC program, the DC program and the transient 
analysis program. The monitor reads and interprets 
the input cards, obtains the program requested (AC, 
DC or transient) from the auxiliary 1311 Disk Storage 
Drive, and transfers control to that program. Upon 
completion of the analysis, control is returned to the 
monitor, which reads any other input cards. Programs 
to be run may be stacked one behind the other, and 
solutions developed without manual intervention. 

Some other features of the system are: 

1. Networks of up to 21 nodes and 60 branches 
may be processed. 

2. Nonlinear circuit elements are handled in the 
transient analysis through linear approximation. 
This is accomplished by automatic switching of 
circuit elements from one value to another in order 
to change the characteristic of the model being used. 

3. The numerical stability of the transient 
solution is not affected by the size of the time step 
selected. Thus, the user can select large time 
steps in order to quickly obtain an indication of the 
transient behavior. In addition, the steady state 
solution can be obtained by selection of a single, 
sufficiently large time step. 

The outputs available to the user are: 

• Node voltages 

• Branch currents 

• Branch voltages 

• Power dissipation 

• Switching times 
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• Sensitivity coefficients which indicate the 
relative effect of variations in the circuit 
elements on circuit output parameters 

• Standard deviation of the node voltages 

• "Worst case" analysis, which yields the max- 
imum and minimum variations of node voltages 
when the worst possible combination of circuit 
element tolerances occurs 

• Coefficients of the circuit equations generated 
by the System 

• Node incident matrix indicating the topology 
of the circuit that the system generated. 

A technique of solution is used in the AC/DC 
programs that substantially reduces the overall 
problem solution time when multiple solutions are 
required because of the varying of circuit parameters. 
The technique requires only a single nominal solution 
of the system equations. Subsequent solutions that 
would normally be required. when circuit parameters 
are varied are obtained through a matrix update 
technique of the nominal solution matrix. 



MACHINE CONFIGURATION 

The minimum system required is a 40K IBM 1620 
Data Processing System, with a 1311 Disk Storage 
Drive, a 1622 Card Read Punch, Automatic Divide 
and Indirect Addressing. 

PRECISION 

The standard floating point format of the IBM 1620 
is used throughout, which carries eight significant 
digits. This precision results in a solution for the 
largest problem accommodated that is well within 
acceptable design accuracy. 

PROGRAM LANGUAGE AND SYSTEM 

The program is written in FORTRAN II D with 
Monitor I. 
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ABSTRACT 



The Electronic Circuit Analysis Program for the IBM 1620 is an engineering 
tool to aid in the design of reliable electronic circuits. The system provides 
facilities for complete DC, AC and transient analysis of an electrical network. 
The user can easily determine the effect of parameter variations on circuit 
performance that might arise due to effects of environment, aging and 
production tolerances. The manual calculations required to analyze circuit 
performance are tedious and time consuming for all but the most trivial circuits. 
The use of the computer relieves the designer of this restriction and greatly 
expands his capability. 

FEATURES 

provides user with a high degree of confidence in circuit design since 
many variations of circuit performance can be "tested" which could not 
have been physically duplicated in any reasonable amount of time. 

System is user oriented, providing means for easy communication and 
complete user control over course of problem solution. 

Simple input language used to describe circuit topology and component 
characteristics. 



Automatic generation of network equations. 

Handles network of up to 21 nodes and 60 branches. 

Automatic modification of parameters or selective modification through 
console. 

Transient analysis of non-linear circuits through automatic piece -wise 
linear approximation. 

Worst- Case analysis of circuit design performed. 

Obtains sensitivity of circuit performance to various parameters. 

Numerical stability of transient solution not affected by size of time step 
selected. 

Voltage and current dependent switching of element values is provided 
in transient program.. 
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DESCRIPTION OF APPLICATION 

The problem of reliable electrical circuit performance faces many designers 
in various industries. It is especially critical in the Aerospace Industries 
where extremely high reliability is desired. In other industries (e. g. , electrical 
and communications) the problem is one of attaining a desired reliability that 
the consumer will tolerate. 

The question to be answered in design of these circuits is what is the tolerance 
required of each individual circuit component in order to obtain an overall 
cricuit performance within certain tolerances. The relationship between 
individual component tolerance and overall circuit performance is not a direct 
one. Thus, it cannot be said that to obtain a maximum of + 1% variation of a 
given circuit output parameter, all components must be at least 1% components 
(tolerances of + 1% from nominal value). In fact, it may be that some components 
will require a + . 01% tolerance others can be as high as + 10%. 

The analysis of the possible combinations of components is a tedious, if not 
impossible task to perfrom manually, and it is not practical to build and test a 
large number of variations of any given circuit. Thus, the digital computer 
provides a means for solving this widespread problem. 



Use 



The Electronic Circuit Analysis Program is written in Fortran II D, 
works under Monitor- 1 on the 1620 and resides on the 1311 disk. It 
is comprised of four programs: A main line program which interprets 
the input language and calls in the specific program requested, the AC 
program, the DC program and the Transient program. With the input 
language, the user describes his circuit, indicates the program he 
wants, the output he wants and which parameters he wants varied. The 
main line program determines the circuit topology, generates the cir- 
cuit equations, calls in the proper program which then provides the 
desired solution. 

MACHINE REQUIREMENTS 

The minimum requirements for the Electronic Circuit Analysis Program is a 
1620 with 40K storage, a 1311 disk drive, 1622 Reader Punch, Automatic Divide 
and Indirect Addressing. 

ANALYSIS PROCEDURE 

The procedure employed by the user to effectively utilize the IBM 1620 Electronic 
Circuit Analysis Program is outlined in Figure 1. The critical ingredient of the 
procedure is the generation of an appropriate equivalent circuit for the electron 
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devices in the circuit under examination (tubes, diodes, transistors). These 
devices must be replaced by an equivalent circuit consisting of passive elements, 
current sources and voltage sources. The use of the Electronic Circuit Analysis 
Program allows the user to incorporate a level of sophistication into the equiva- 
lent circuits that renders a hand analysis of the result totally impractical. For 
a discussion of equivalent circuits for electron devices, the reader is referred to: 

SheffLer, H. S. , Terry, F. R. , "Description and 
Comparison of Computer Methods of Circuit 
Analysis", North American Aviation, Inc. , Auto- 
netics Division, Manual EM-6839. 
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Figure 1. Analysis Procedure 
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I. Introduction 



The need for and usefulness of an integrated system of programs to do 
network analysis has already been shown. In order for the system to be 
most effective in carrying out its task, it must "speak" to the user in 
his own jargon. The circuit design engineer must be able to visualize 
a direct connection between the output which he obtains from the system 
and the input data presented to the system. 

For this reason the SNAP language program was developed. The input 
data which must be prepared by the engineer is as close as possible to 
the format the engineer would normally use to carry out hand analyses of 
his circuit. 

The input data could, perhaps, be handled in fixed- fields on data cards 
as FORTRAN E-D normally does. But the engineer's task seems over- 
bearing when one considers that a 20 node-60 branch circuit could contain 
over 470 separate items to be input. This in itself could discourage the 
use of the system by all but the most earnest engineers. The present 
input language does not require the engineer to learn something entirely 
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new to him, but instead allows him to use directly what he is already 
familiar with. The teaching of the language to prospective users is, 
to coin a phrase, a "SNAP. " 

n. Preliminaries 

The approach a person would take in order to prepare data for the SNAP 
language program would be, as mentioned above, very close to the normal 
hand analysis procedure. The engineer would prepare from a schematic 
diagram an equivalent circuit using the basic building blocks (i. e. , 
. resistors, capacitors, inductors, voltage and current sources, trans- 
conductances, and mutual inductances) in place of special devices such as 
transistors, diodes, vacuum-tubes, etc. . The engineer must decide upon 
element values including tolerances, if needed, for D. C. , alternate values 
in the case of switching in transient and phase angles and frequency for A. C. 
This particular phase — the gathering of actual and representative values — 
is usually the most difficult for any analysis. 

Ingenious uses of the so-called "standard-branch" can in many cases lead 
to a simplification in the input data. The standard branch is defined as 
having each end terminate at a separate node and the branch between these 
two nodes can be one of the combinations shown in Figure 1. 
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Y can be any of the following: 

1. Resistor or Conductor 

(expressed in ohms or mhos) 

2a. Capacitor - (farads) 

2b. Capacitor with leakage 

(farads and ohms, or mhos) 



3. Inductor (henries) 



Figure 1 
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The assumed positive directions are as shown in the figure and will appear 
positive in sign on the data cards when oriented as shown. Note that the 
assumed positive direction of branch current goes from the INITIAL node to 
the FINAL node. The specification as to which node is initial and which is 
final is important in many cases and is part of the input information. 



Once the engineer has represented his schematic in terms of the standard 
branch and has decided upon positive current directions, he then consecutively 
numbers the nodes and branches, as well as transconductances and switches. 

in 



The common node (usually ground) is referred to as node zero. This done, 
he is ready to write down the input data to be key-punched. 

EI. Specific Input Card Types 

The engineer is now prepared to make use of the specific input card types 
with which he can properly describe his circuit to the SNAP system. There 
are three types of SNAP language statements: COMMANDS, CONTROL and 
DATA. A COMMAND specifies what analysis program is to be called and 
when to call it. A DATA statement specifies the circuit configuration and all 
of the circuit parameters. The class of statements called CONTROL cards 
instruct SNAP on what to do, but in a more specific sense than a command. 

A. COMMANDS 

1. DCNODE 

This command specifies that a DC analysis is to be made 
on the data which follows it. 

2. ACNODE 

An AC analysis is to be made on the data which follows 
tnis card. 

3. TRANSIENT 

A transient analysis is to be performed. 
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4. MODIFY 

This command indicates that the data that follows modifies 
the previous analysis and that the previous analysis is to 
be either repeated or continued with the modified data 
incorporated. (Note: only 1 MODIFY command is allowed 
before any EXECUTE command). 

5. EXECUTE 

This command will cause execution of the commands 
preceding it which have not yet been executed. 

6. CONSOLE 

This command causes execution of the commands preceding 
it as with the EXECUTE command, but after execution is 
completed, will re-define the input device as the CONSOLE 
typewriter for manual entry. 

The EXECUTE command will not re-define the input device. 
The SNAP program initially expects input via the card 
reader and will continue to do so with EXECUTE commands 
until a CONSOLE card is read. Thereafter, EXECUTE 
commands will leave the console typewriter as the input 
device. 

7. READER 

This command is similar to the CONSOLE command except 

m 
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it defines the card-reader as the input device after 
execution of the preceding unexecuted commands. 

8. END 

This command defines the end of a SNAP run. It causes 
the program to return control to the 1620 Monitor I 
Supervisor. It would not be used until a user were completely 
through with circuit analysis. 

9. 1USER 

This command as well as the two following can serve to 
call special user routines which can be compiled on the 1620 
Monitor I, added to the SNAP package and called by use of 
these commands. 

10. 2USER 
U. 3USER 



Examples of commands are shown in Figure 2 and are summarized in Table I. 
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Figure 2. 



Command 




DCNODE 


Calls D. C. Analysis Program 


ACNODE 


Calls A. C. Analysis Program 


TRANSIENT 


Calls Transient Analysis 


MODIFY 


Denotes data to follow as modifying previous 
analysis and creates linkage to repeat pre- 
vious analysis with modified data 


EXECUTE 


Signals end of data for analysis called for 
and starts execution 


CONSOLE 


Same as execute, except after execution 
returns control to typewriter 


READER 


Same as execute, except after execution 
returns control to card reader 


END 


Signals end of SNAP jobs . 


1 USER 

2 USER 

3 USER 


These commands can 
call special user com- 
piled programs. 



Table 1. 
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B. DATA STATEMENTS 

There are two general classes of data statements. The first classification 
specifies data from the equivalent circuit such as circuit topology, 
element values and initial conditions. The second group specifies data 
which applies to the circuit in a more general way. In this classification 
lie the frequency to be used in an A C circuit, the time step to be used 
for a Transient solution, as well as solution criteria such as is needed 
for Transient. 

1. Equivalent Circuit Data Statements 

These data statements contain information of two general 
types. The first type identifies the kind of data statement 
and its serial number and the second type identifies the actual 
data by prefacing the data with a data label. 

a. Data Statement Card Identifiers 

1. B - identifies a branch card. 

Branch cards contain information about the individual branches 
within the circuit to be analyzed. 

2. T - identifies a bransistor card. 

A transistor card specifies a transconductance or current gain 
between two branches in the equivalent circuit. 
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3. M - identifes a mutual card. 

Mutual cards specify the values of mutual inductance which 
exists between two inductive branches. 

4. S - identifies a switch card. 

Switch cards contain information about the switches in a 
circuit to be analyzed by the transient program. 

b. Data Labels 

Twelve different types of data could possibly be 
found on a data; statement card. The following list 
associates the parameter with the particular identifier 
and indicates on which type of data statement card a 
given parameter may occur. 



1. 


Nodal data: 


"N" 


(B-cards only) 


2. 


Branch data: 


"B" 


(T-M-or S cards) 


3. 


Resistance data: 


""R" 


(B-cards only) 


4. 


Conductance data: 


"G" 


(B-cards only) 


5. 


Voltage data: 


"E" 


(B-cards only) 


6. 


Current data: 


MJMt 


(B-cards only) 


7. 


Inductance data: 


"L" 


(B-or-M cards) 


8. 


Capacitance data: 


"C" 


(B-cards only) 


9. 


Initial Capacitor Voltage: 


"EO" 


(B-cards only) 



10. Initial Inductor Current: 

11. Trans conductance data: 

12. Current Gain data: 
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"10" (B-cards onV) 

"GM" (T-cards only) 

"BETA" (T-cards only) 



Examples of Data statements are shown in Figure 3. 

Card 

columns 1 7 





TRANSIENT 


Bl 


N= (0,1), R= l.E + 6 


B2 


N = (1, 2), C = 5. 3E - 6, EO = 150. 


B3 


N = (2,5), G = .00125 


B4 


N = (5, 6), L = 5. 7E - 3, IO - . 0015 


Tl 


B = (1,3), BETA = 65. 


T2 


B = (2,3), GM = l. E -2 


SI 


B = 5, (7, 1, 9), OFF 




DELTA = . 5E - 6 




EXECUTE 




END 



Figure 3 
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General Data Statements 

These data statements contain only two items: a data label 
and a value. 

a. " FREQUENCY" 

This data card specifies the frequency to be used in 
an AC analysis. Frequency is specifies in cycles. 

Note: All of the following data cards are for use with 
the transient analysis program. 

b. "DELTA" 

Specifies the time increments to use in the transient 
analysis program. . This would normally be specified 
in units of seconds. 

c. "MAJOR" 

Specifies the number of time steps which should pass 
before an output should occur. The user may want to 
integrate with a small time step but may only want out- 
put at every tenth step. MAJOR is assumed to be 1 by 
the program unless otherwise specified. 
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d. "1 ERROR" 



Specified the maximum allowable error in the sub- 
stitution check. A check is made to see that all 
of the currents into each node sum to zero. Nor- 
mally, due to round- off, they will not sum to zero. 

1 ERROR allows the user to control the limit beyond 
which a message will be output indicating that the 
criterion has not been met. 1 ERROR is assumed 

. to be 1 milliamp unless otherwise specified. 

e. "2 ERROR" 

When the transient analysis program is searching 
for the exact time at which a switch should be 
actuated, 2 ERROR provides a tolerance on how exact 
this time should be found. Unless specified by a 

2 ERROR card, the program uses a value of . 001 
indicating that the switching time will be found to 
within. 1% of the time step. 

f. "3 ERROR" 

Accuracy in results can be increased by decreasing 
the time step under certain conditions of switching 




when inductors are involved in a circuit. 3 ERROR 
specifies the reduction to be made. Its value is 
assumed to be 1. by the program unless otherwise 
specified. This will inhibit this feature which is the 
normal procedure for most circuits. 

g. "INITIAL TIME" 

Specifies the starting time printout. DELTA will 
increment time from this base. This can be useful 
when re-analyzing a certain portion of a transient 
response when this portion occurs after time equals 
zero. The program assumes 0. unless otherwise 
specified. 

h. "FINISH" 

Specifies the time at which to terminate a transient 
analysis. This can allow the user to leave the com- 
puter unattended knowing that his job will not run 

endlessly. The value of FINSIH is initialized to 
49 

1. * 10 unless otherwise specified by a FINISH 
card. 
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i. "SHORT" v 

These two cards are available to the user 

j. "OPEN" 

to specify quantities used to calculate 

the Time = 0.0 voltages and currents. At time zero 

capacitors act as short circuits and inductors act as 

open circuits. These two values are substituted in 

place of the actual values when the transient analysis 

program is calculating these initial conditions. The 

Program assumes values of . 01 for short and 
7 

1.0 x 10 for open unless cards are input to the con- 
trary. If the user switches these values, that is, if 

7 

he sets short to 1.0 • 10 and open to .01 he can obtain 
a steady state solution corresponding to the state of 
the circuit at t = CO . 

Examples of the use of Data Cards are shown in the 
sample problems in APPENDIX A. 

The following Control Cards may be used with the SNAP Program (See 
TABLE III for a summary): 
1. "PUNCH" 

The SNAP Program is initialized to print the output data, 
but the entering of this control card will redefine the card 
punch as the output device. 
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2. "PRINT" 

This card is used if a user wants both PRINT and PUNCH 

as output options. The punch control card will inhibit printing 

unless this card is used. 

3. "SENSITIVITY" 

The D. C. analysis program is prepared to calculate the 
sensitivity of each node voltage with respect to a small change 
in each circuit parameter. In addition, the actual partial 
derivatives are output. This control card initiates this output. 

4. "WORST-CASE" 

When given the proper tolerance data, the D. C. Program will 

make use of the partial derivatives mentioned above to calculate 

the worst-case maximum and mimimum node voltages. This 

card allows the user to get this information. 

Note: the partial derivatives are necessary in order to calculate 

these worst-case conditions, therefore, the SENSITIVITY card 

is not required as input information. The program will automatically 

calculate the partial derivatives and sensitivities. 

5. "STADEV" 

The DC Analysis will calculate the standard deviation of each node 
voltage when the parameter tolerances are input and this card is 
entered. The sensitivities and partial derivatives are also required 
for these calculations so that the SENSITIVITY micro -data card is 
again not required. i q <f 




Control Card 


Use 


Analysis 


PUNCH 


Enables results to be punched. 
Inhibits Printing. 


ALL 


PRINT 


Used only when punch inhibits 
print and printed output is also 
desired. 


ALL 


SENSITIVITY 


Enables DC to perform sensi- 
tivity and partial derivative 
calculations. 


D.C. 


WORST-CASE 


Enables Worst- Case calculations 
as well as sensitivity. 


D.C. 


STADEV 


Enables Standard deviation cal- 
culations as well as sensitivity. 


D. C. 
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Equivalent Circuit of Single Stage Common Emitter Amplifier 
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N=( 1,2) ,R^6000.0 
N=(2,0) ,R^1000.0 
N=( 1,3) ,R=i2000.0 
N=(3,4) ,R=43000.0 
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EXECUTE 
EXECUTE PROGRAM. 
DC CALLED. | 
DETMA ! i ' j 

TOPOLOGICAL MATRIX DIMENSIONED 
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NODES 
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-1 
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-1 




1 - 4 0.1000006I7E 03-0. 16666667E-03H-0. 5000000j0£-03 0. I 

1 - 4-0. 1666666i7E-03 0. 32 1 79487E-02j 0. j -0. 20j5 1282 1E-02 

1 - 4-0.50G0C0Cj0E-03 0. 86153845E-01 0. 52 32 55811 E-G 3-0. 861177 1Q0E-0 1 
1-40. 1 -0. 88205 127E-0 1-0. 232558 14E-04 0. 90228382E-0 1 



EFCURR! K= 1 
LMVBM K= 1 

EQUIVALENT CURRENT VECTOR 



NODE NO. 



CURiRENT 



1 

2 
3 

4 I 
MISER 

NODAL SOLUTION MATRIX 



0.200OC00OE 04 
0.11897436E-02 
0.49969229E-01 
-0.5115I8973E-01 



NODE ROW NODE CGLS 



910 1 23456789J0 1 2 345678910 1 23456789 



01 2345 67890 I 23456 7 39|0 



'' |VJ ■■'■■'■■"'■'■■■^'■MiLi.Mm.^Biil!*, 



fflli1li^i^]MV^m L °' ^ ■'• L^1 ^^" mm ^■^ m '^ L " J ■^■'-^■^"^■ 



6789 01234 5 6 



1-4 0.99999590E-02-0. 14312440E-01 0.99640164E-02 0. 9 1912452E-02 
t 6i8 0il3?644^1E^Q2 5 0.82i32&ii2E 7 03; 2 . 8T00B%$1 E 2 OQ 5 7 . 19I547309E 7 02 
1 - 
1 - 



3 
4 

PMVBM 

NODE VOLTAGE SOLUTION VECTOR 



4 0.473438a8E-02-0.31369275E 04 



4 0. 134679d4E-02 0.80399203E 031 0.1364135I7E 01 0.30|664095E 02 



NODE NO- 



0, 
0. 



0, 



VOLjTAGE 



19999928E 



1 

2 I 

3 I 

* ! 

LMVBM |K» 2 

BRANCH VOLTAGE RECTOR 



27758129E 
il59|i434E 
2149|5623E 



02 
01 
02 
01 



BRANCH NO 



VjOLTAGE 



0.17224115E 02 

0.27758129E 01 

0.84084945E 01 

0.94|418712E 01 

0.21495623E 01 

0.62625059E 00 

0.19999928E 02 



COIL VjOLTAGE VECTOR 



BRANCH NO. 



VOLTAGE 



0.17224115E 02 
0.27758129E 01 
0.84084945E 01 
0.94j418712E 01 
0.21I495623E 01 
0.4^250604E-01 
-0.71I763992E-04 



EFCURRj K= 2 ! 
BRANCH CURRENT viECTOR 



•J 



BRANCH NO. 



CURRENT 



0.28|706858E-02 
0.27758129E-02 
0.42I042472E-02 
0.42j042457E-02 
0.42!991246E-02 
0.94I873034E-04 
-C.7lj763993E-02 



COIL CURRENT VECTOR 



BRANCH NO 



0.19925266E 04 0.18I317456E 04 



910 1234567 



2 1 2 3 4 

3 



CURRENT 



0.28706858E-02! 

G* 27758 129£^G2|o , 2 3 4 5 6 7 8 9 

0.42042472E-02 



QiQl6 7 



1 2 3 4 5 6 7 



910 1 23456789 



1234 5 6739,0 



4 0.42042457E-02 

5 0.42991246E-02 

6 7 8 9 jO 1 2 34 5 6 7 8 £j 01 2 3 4 5 Q ^ 9 4 8 7 3 3 4 E* 2 * 04 01 23456789 

! "fl -0.71763993E-02 



POWER 



LOSS IN BRANCHES 



BRANCH NO. 



1 

2 
3 
4 
5 
6 
7 

PARSHLj 



POWER 



0.49445022E-01 
0.77051372E-02 
0.3^351389E~0i 
0.39695946E-01 
0.92412360E-02 
0.43|879351E-05 
0.51j500706E-06 
1! 



01 23456789 



PART I DERIVATIVES 



NODE VOLTAGE WITHRESPECT T!0 RES. IN BRANCH 



SENSITIVITY 



0.69793257E- 
•0.i23633200E 



0. 
0. 



90051467E- 



NODE N,Q. PART IjALS 

1 0. 116322lioE-07 

2 -0.3938866j7E~03 

3 ! 0.15008578E-02 

4 -0.38466744E-03 
PARSHL 2 l| 



PARTIAL DERIVATIIVES 



NODE VOLTAGE WITHRESPECT TlO RES. IN 



06 
01 
01 



23080046EH01 



NODE NO. PART I ALS 

i ; 

1 -0.39728694E-07 -0. 

2 | 0.22852188E-02 0. 

3 h0.8707523|7E-02 -0. 

4 ! 0.22317314E-02 0. 
PARSHLj 3 1 



BRANCH 



SENSITIVITY 



39728653E- 
22852188E- 
87075235E- 
22317314E- 



06 
01 
01 
01 



PARTIAjL DERIVATIIVES I 

NODE ViGLTAGE WITHRESPECT T|0 RES. IN 



NOOE NO. PART I ALS 

1 '! 0.7555589I0E-1O 

2 -0.18261319E-05 

3 i-0.418852lHlE-02 

4 I-0.2864750I6E-05 
PARSHL! 4 l! 



PARTIAL DERIVATIVES 



SENSITIVITY 



0.151U178E-08 
-0.36522638E-04 
-0.i83770541EH01 
-0.57295010E-04 



BRANCH 



6 7 8 9:012 3 



NQDEs VOLT AGE* WITHRESPECT s T|Q i R6S* 6 I ^BRANCH 



01 23456 7 89 



01 23456789 



1 23456789 !.0 



3 



4|0 1 23 456789 



1 23456789 



1 234 5 6739|0 



NODE NO. PART I ALS SENSITIVITY 

6/890 1 73456789:0-1 234567 8 9 jO 1 73456789 jo 1 2 3 4567890 1' 23456789:0 1 23456789 jOl 23456789 

1 0.75556279E-1O 0.32489199E^07 

2 -0.18261312E-05 -0 . 7852364 1E-03 

3 0. 15720065E-04 0.I67596278E-02 

4 -0.28647495E-05 ~0 . 123 1 8423fc-^02 
PARSHL 5 1 ! 



PARTIAL DERIVATIVES j 
NODE VOLTAGE WlTHRESPECT TO RES. IN iBRANCH 5i 



1 23456789 jO 



NODE NO. PART IALS 

1 0.79028616E-07 

2 0.16807263E-03 

3 0.15749805E-01 

4 0.26365753E-03 
PARSHL 6 1 



SENSITIVITY 

0.39514307E-06 
0.84036316E-03 
0.j78749022E-Ol 
0.13182876E-02 



PARTIAL DERIVATIVES 

NODE VOLTAGE Wl THRESPECT TO RES. IN BRANCH 



NODE NO. PART I ALS 

1 0.17422980E-08 

2 0.37502339E-05 

3 ! 0.34721164E-03 

4 -0.8898999j7E-04 
PARSHL 7 I 



PARTIAL DERIVATIVES 



SENSITIVITY 

; i 

0.:84937024EH08 
0. 18282400E-|04 
0.16926567EH02 
■0.43382623EH03 



NODE VOLTAGE W I THRESPECT TlO RES. IN 'BRANCH 



NODE NO. PART I ALS 

1 h0.71763699E-02 

2 j-0.98779051E-03 

3 -0.33975864E-02 

4 h0.9665i630E-03 
PARSHLi 1 2 



| SENSITIVITY 

C71763697E-06 
•0.98779049E-07 
•0.33975864E-06 
0.96651628EH07 



PARTIAL DERIVATIVES 
NODE VOLTAGE WITHRESPECT TO GAIN 



NODE NO 



PART IiALS 



BRANCH 6 TOl BRANCH 



1 -0.3574113I4E-04 

2 | 0.86383292E 00 

3 -0.74362179E 01 

4 'j 0. 13551408E 01 
PARSHL I 3 



SENSITIVITY 

I 

•0.I30792361E-07 
0.74422528E-I03 
■0.64065876E-02 
0.11675059EH02 



6 7 8 9 jo 1 23456789 ; 1 2 3 456739.0 ! 23 4 5 6 7 8 9 j 1 23456789jo 1 23456789|0 1 2 3 4 5 6 7 8 9|0 1 23456789 



PAR SHL 2 3 

4 . ft • ,. • ..: - . 4 •, 6 ' 8 '-' • 2 4 S » * 8 V > 2 "> 4 •> 6 * 8 '■' : 2 3 4 5 A / 8 V i 2 3 4 5 6'7 8 9 I 2 3 4 5 6 / .8 9 i 2 i 4 3 6 / 8 •' 

PARSHL 3 3 

PARSHL 4 3 
PAR SHL 5 3 
PARSHL 6 3 
PARTIAL DERIVATIVES 

NODE VOLTAGE WITHRESPECT TO VOLTAGE SOURCE IN BRANCH 6| 

NODE NO. PART I ALS SENSITIVITY 

1 0.18364523E-04 -0. 15821743E-05 

2 0. 39528981E-01 -0 . 34055737E-02 

3 0.36597508E 01 -0 • 3 1 530160E— 00 

4 -0.93799041E 00 . 808 i 148QE-0 1 
PARSHL 7 3 

PARTIAL DERIVATIVES 

NODE VOLTAGE WITHRESPECT TO VOLTAGE SOURCE IN BRANCH 7 

NODE NO. PARTIALS SENSITIVITY 

1 0.99999591E 00 -0.99999592E 02 \ 

2 0. 13764431E-00 -0.13764431E 02 

3 0.47343888E-00 -0.47343888E 02 

4 0. 13467984E-00 -0. 13467984E 02 ! 
PARSHL 1 4 

PARSHL 2 4 

PARSHL 3 4 

PARSHL 4 4 

PARSHL 5 4 

PARSHL 6 4 

PARSHL 7 4 




I 

: / . a « 6 7 8 9 . 2 3 4 5 6 7 8 9 I i 3 4 6 7 8 9 . t 2 3 4 5 678901 2 3 4 5 6 / 8 9 I :> < 4 S 6 / 8 C 




1 7 







ACNODE 


Bl 




N (1,0), G= .25, C = .1 


B2 




N (2,0), R= 5. 


B3 




N (0,1), L = 2.0, E = 1/0 


B4 




N (1,2), L = 3.0 


Ml 




B (3,4), M= 1.0 






FREQUENCY = 1. 59 






EXECUTE 









T 



jO 1 23 4 5 6 7 8 9 jo I 2 3 4 5 6 7 8 9 jO 1 23456789 jo ) 23456789 jO 1 23456789 jo 1 23456789 |0 1 23456789 O ! 234567 8 




Bl 


NMl,0),6s. 


25,C=.l 


'\2 


NM2,OJ,R.= 5 


• 


13 


NN t Of 1 ) , L =2 


. ,E=1 ./0. 


'34 


Nb( 1,2) ,L=3 


• 


* 1 


BM 3,4) ,M=1 


• 




FJ*EGUENCY = 


1.59 




EXECUTE 





UXECUTfe PROGRAM. 
AC CALLED. 



NODE VpLTAGE SOLUTION VECTpR 



MAGNITUDE 
8j. 146E-03 
li.044E-Cl 



PHASE. ' 
1. 740E+02 
9.049E+D1 



: BRANCH VOLTAGE VICTOR 



MAGN I TCGE 
8.146E-03 
1.044E-G1 
8!. 146E-Q3 
1.055E-0L 



l. 740E + f)2 
9.049E+D1 
5.936E-00 
9. 39GE+01 



PHASE 



0123456789, OI2345678ViO 12345678901 234567890 I 23456/8901 2 34 5 67890 1234567890! 2 3 4 J 6 * H v 



TRANSIENT ANALYSIS EXAMPLE 
6. LCR Voltage Regulator 

The objective of the circuit in this example is to regulate the voltage V3 at 
12,0 volts. This is done by utilizing the cut-off and saturation characteristics of a 
transistor and a diode. If V3, is less than 12. Ov, a comparator senses this and turns 
the transistor on. This cuts off the diode and builds up the voltage on the capacitor. 
When V3S12.0 volts the transistor is turned off The inductor tends to resist, a change 
in current. The return path for the current is through the diode which goes from the 
blocked condition to full conduction. The voltage drops on the capacitor, which is 
sensed by the comparator. There is a hysteresis voltage drop associated with the 
comparator of 0.3 volts. Hence, when V3 - 11. 7< 0, the comparator turns the 
transistor back on and the process is repeated. 



* ♦ 



j Q 

T 15 




O e r =12 u: 



Fig„ 6a LCR Voitage Regulator - Schematic Diagram 



.2 mhy 

( 1, 1M) w ^ 

rA/WSrr Vv^A,V^ r Dl^l; 

B3 B4 



V3 



Bl 



15v. 



t 



2Mf 



B5 



V4 



B2 
(1M,05) 

1 



I *> B6 

-0_ 



B7 
6Q 



Fig. 6b LCR Voltage Regulator - Equivalent Circuit 



< 



B8 



(12.0,11 



4 5 6 7 8. 9 > 2 y 4 s 6 7 H 9 I 7 3 4" *> 6 7 8 t> ■ 2 ! 4 r i 6 "' S '> ! 2 '4 5 6 7 8 v I 2 i 4 .i 6 '■ « 9 1 2 3 4 5 6 7 8 9 1 2 :i 4 :> 6 7 8 9 I 2 I 4 



TRANSIENT 

C 
C 

C LCR VCLTAGE REGULATOR 

C 

C 

Bl N*(0 f l),R*(.l,l.E6) f .E=(i5.»i5.) 

B2 N»(0tl) f R=(i.E6i.5) 

B3 N=(l t 2)tR=.l 

BA N= { 2, 3 ) ,L=.2E-3 

B5 N*(3,A),C=.2E-5 

B6 N=(0,A),R*i. 

B7 N=(0,3),R=6. 

B8 N=(3,0) ,R=M 1.E6.1.E6) ,E={-12. ,~11.7) 
SI B-8 f ( L Y 2t 8) ,OFF 

DELTA*1.0E-6 

KAJOR*20 

lERROR=1.0E-2 

FINISH=200.0E-6 

EXECUTE 
EXECUTE PROGRAM. 
TR CALLED. 



T = 0. 

V 1 = 0.1500E 02 

V 2 « 0.150CE 02 

V 3 = 0.1167E-0A 

V A = 0.1155E-0A 



1 = 0.1788E-0A 

2 = -0.1500E-0A 
3*0. 

A = 0.1500E-05 

5 * 0.U5 5E-0A 

6 » -0.1155E-0A 

7 * -0.19A5E-05 

8 * -0.1200E-OA 



T * 0.2C00E-04 

V 1 * 0.1A87E 02 

V 2 - 0.1A7AE 02 

V 3 * 0.AA36E 01 

V A = 0.5578E 00 



i 




0.1297E 


Oi 


2 




-0.1A87E- 


-OA 


3 




0.1297E 


01 


A 




0.1297E 


01 


5 




0.5578E 


00 


6 




-0.5578E 


00 


7 




-0.739AE 


00 



ssSs ^o-rG« 756AEtQ5 2 



* 5 6 ." 89:0! 2 3 456789! 01 



4 , 6 f 8 teO 0*4CC0E' U -04 2 ' 4 ■.' 1 f' ? ■■ I' ■■ ft « ; » ' " ' ^ i 6 7 8 9 ■ 2 . 4 6 ■ 3 * 2 



V i = 0.1480E 02 

V 2 = 0.1459E 02 O 

V 3 = 0.9586E 01 

V 4 = 0.4516E-00 



1 * 0.2049E 01 

2 = -0.1480E-04 

3 = 0.2049E 01 

4 = 0.2049E 01 

5 = 0.4516E-00 

6 = -0.4516E-0G 

7 = -0.1598E 01 

8 = -0.2414E-05 



T = 0.5378E-04 

V 1 = 0.1477E 02 

V 2 = 0.1454E 02 

V 3 = 0.12CCE 02 

V 4 = 0.3031E-00 



I 1 = 0.2303E 01 

I 2 « -0.1477E-04 

I 3 = 0.2303E 01 

I 4 - 0.2303E 01 

I 5 = 0.3031E-00 

I 6 = -0.3031E-00 

I 7 = -0.2C00E 01 

I 8 = -0.4399E-10 



SWITCH 1 IS ON 



T = 0.5378E-04 

V 1 = -0.1152E 01 

V 2 = -0.1382E 01 

V 3 = 0.1200E 02 

V 4 = 0.3026E-00 



1 » 0.1615E-04 

2 » 0.2303E 01 

3 = 0.2303E 01 

4 « 0.2303E 01 

5 =* 0.3026E-00 ^ 

6 = -0.3026E-00 

7 = -0.2000E 01 

8 = 0.3025E-06 



0.6C00E-04 



T * 



AH 



V - 6 1 s ^ 1*0 » 94 5 56 " 00 2 :< ■'« o 7 S 9 : ' 2 3 4 3 6 8 « 2 ;■ 4 r > 6 ' 8 ■ ? ! 4 6 ' 8 v ■ 1 2 i 4 •> 6 ' 8 * 

V 2 = -0.1135E 01 



V 
V 



3 a 
A = 



C.1183E 02 
-0.8035E-01 



1 « 0-1595E-0A 

2 = C.1891E 01 

3 a 0.1891E 01 
A a 0.1891E 01 

5 = -0.8035E-01 

6 * 0.8035E-01 

7 = -0.1971E 01 

8 = 0.1279E-06 



T = 0.6116E-04 

V 1 a -0.9080E 00 

V 2 = -0.1090E 01 

V 3 * 0.1170E 02 

V A = -0.13A0E-00 



1 a 0.1591E-0A 

2 * 0.1816E 01 

3 = 0.1816E 01 
A * 0.1816E 01 

5 a -0.13A0E-00 

6 * 0.13A0E-00 

7 ». -0.1950E 01 

8 a 0.3791E-10 



SWITCH 1 IS UFF 



T a 0.6116E-04 

V 1^ 0.1A82E 02 

V 2 a G.1A6AE 02 

V 3 ^ 0.1170E 02 

V A = -0.1338E-00 



1 a 0.1816E 01 

2 = -0.1A82E-0A 

3 = 0.1816E 01 
A = 0.1816E 01 

5 a -0.1338E-00 

6 a 0.1338E-00 

7 = -0.1950E 01 

8 * -0.3011E-06 



T * 0.7991E-04 

V 1 = 0.1A79E 02 

V 2 ai 0.1A58E 02 

V 3 = 0.1200E 02 

V A = 0.87A5E-01 



1 , 6 1 a a „ 0.2087E 01 , 
I 2 = -0.1A79E-04 



123456789,0 123 



3 = 0.2087E 01 

4 * 0.2087E 01 

5 6 5 8 *! 1 i 8745E*"!0 1 2 345678 9 |o 1234567890 123456789 

6 H -0.8745E-01 | 



7 4 -0.2000E 



8 * 0.1895EniO 



SWITCH 



01234567890123 



1 = 

2 - 

3 = 

4 *j 



1 4 

2 H 

3 =j 

4 ^ 



IS ON 



T M 0.7991E^04 



-0.1044E 
-0.1252E 
0.1200E 
0.8730E- 



01 



01 
01 
02 
01 



= 0.16046^04 



0.2087E 
0.2087E 
0.2087E 



01 
01 
01 



4 0.8730EH01 
-0.873CE-i01 
-0.2CC0E |01 
0.3007EH06 



T 4 0.8000E- 



04 



0.1041E 01 
0.1249E 01 
0.1200E |02 
0.8142E-01 



0.1604EH04 
0.2081E 101 
0.2081E ;01 
0.2081E 101 
0.8142E-01 
0.8142E-01 
0.2CCCE 01 
0.2987E-06 



T * C.8386E-04 



J _ 



0.9145E 00 
0.1097E |01 
0.117CE 102 
0. 1209E->00 



1 4 0.1591E-04 

2 4 0.1829E 01 

3 4 0.1829E 01 

4 = 0.1829E 01 

5 =4 - 

6 = 

?8 ro"T 

8 * -0.3254E-10 



•0.12C9E-00 
0. 1209E-00 
Q# 19§0E 9 01 



01 23456789 



01 23456789 



r< — "^i 



2 3 4 5 6 



8 9 ; 1 234567 



5/3 

8 9 1 



23456789 



01 23456789 



12 3 4 5 6 7 



1 2 3 4 5 6 7 B 9i0 1 2 3 S^fTCH 1 2 3 J, 5 6 7 f Si° OFF 



4 5 6 7 8 9 



T 4 0.8386E^04 

i ; 

1 * 0.1482E 02 

2 4 0.1463E 02 

3 4 0.1170E 02 

4 H -0.1207EH00 



1 0.1829E 01 

2 * -0.1482E-04 

3 ^ 0.1829E 01 

4 ^ 0.1829E 01 

5 * -0.1207E-00 

6 a 0.1207E-00 

7 = -0.1950E 01 

8 - -0.3011E-06 



0123456789 123456789 jO 123456789 10 123456789 jO 1234 



T * 0.100GE-03 

1 4 0.1479E 02 

2 = 0.1459E 02 

3 4 0.1192E 02 

4 4 0.7623E-01 



1 aj 0.2063E 01 

2 =| -0.1479E-|04 

3 4 0.2063E 01 
A aj 0.2063E 01 

5 * 0.7623E-01 

6 4 -0.7623E-I01 

7 ^ -0.1987E jOl 

8 4 -0.7665E-07 



T =! 0.1016E-I03 



1 4 

2 4 

3 ^ 
A 4 



SKITCN 



0.1479E |02 
0.1458E 02 
0.1200E 02 
0.8494E^01 



0.2085E !0i 
0.1479E-04 
0.2085E 101 
0.2085E 01 
0.8494EH01 
0.8494E-01 
0.2C00E 01 
C2325E-10 



IS ON 



0123456789:01 234567891012345678901 234567 ° 9 j 



A 14 

3 4 5 6 7 8 9]0 1 23456789|0 1 23456789 



1 23456789 1 2 3 4 



T = 0.1016E-03 



,0 1 2 3 4 5 6 7 8 90 1 2 3 y 5 6 ^8 S^O^OilOAZE 

'! V 2 ^ -0.1251E 
I V 3 = 0.12G0E 

i V 4 



1 * 

2 = 

3 = 
A = 

5 * 

6 = 

7 = 

8 = 

I 



Q^23456789 

Oi 
02 

0.8479EH01 



0.1604E^04 
0.2085E 01 
0.2085E 01 
0.208 56 Oi 
0.8479E^01 
-0.8479EH01 
-0.2000E |01 
0.3007EH06 



i 



0.1056E^ 



03 



1 a| -0.9136E 00 

2 J -0.1096E jOl 

3 *i 0.1170E !02 
-0.1228EH00 



1 = 

2 = 

3 = 

4 = 

5 = 

6 = 

7 = 

8 = 



S W I ICh 



0.1591E-H04 
0.1827E 01 
0.1827E 01 
0.1827E 01 
0.1228EH00 
0.1228E-H00 
0.1950E 01 
0.5162E-10 



T J 0.1056EH 



1 = 

2 = 

3 = 

4 = 



IS 



0.1482E 
0.1463E 
0.1170E 



012345678 9j0 1 2 3 4 5 6 7 8 9 iO 1 23456789 



OFF 



03 



02 
02 
02 



-0.1226E-00 



0.1827E 01 
-0.1482E^04 



0.1827E 
0.1827E 

-0.1226E- 
0.1226E- 

-0.1950E 



8 M -0.3010E- 



01 
01 
00 
00 
01 
06 



012345678901234 



I T = 0.1200E-03 i 

I ! I 

i : ! 

V 1 * 0.1480E 02 I ! I 

V 2 = 0.1459E 02 | j 

jO 1 2345678910 1 2 3 Vs 6 3s?^0 l0*l 1 846 9 : Q2 2 3 4 5 6 7 8 9|o 1 2345678 9T0 1 23456789j0 1 23456789 

V 4 = 0.6492E-01 



0123456789 



1 2345678910 t 2 



5 6 7j_8 fejO 1 0«2G39E 9 Ol 2 3456789 ; 1 2 3456789 io 1234567890123456789 

2 =1 -0.1480E^04 I I 

3 A 0.2039E 01 ; | 

4 *| 0.2039E 101 

5 A 0.6A92EH0i 

6 *j -C.6492E^01 

7 H -0.1974E |0i 

8 4 -0.1576EH06 



T * 0.1235E-03 



V 


1 




V 


2 


= 


V 


3 




V 


4 





SWITCH 



T = 0.1235E- 



1 * 

2 = 

3 = 

4 = 



0.1479E 
0.1458E 
0.1200E 
0.8534E- 



0.2085E 
-0.1479£^ 



02 
02 
02 
01 



01 
04 



«i 0.2085E jOl 
0.2085E 01 
0.8534E^0l 
-0.8534E-01 
-0.2000E 01 
4 -0.2384EH12 



IS ON 



■0„1043E 
•0.1251E 
O.12C0E 
0.8518EH 



0.1604E- 

0.2085E 

0.2085E 



04 
01 
01 

0.2085E 101 



03 

01 
01 
02 
01 



0.8519E- 
-0.8518E- 
-0.2000E 

0.3007E- 



0.1275E- 

0.9136E 
0.1096E 
0.1170E 
0.1228E- 



01 
01 
01 
06 



03 

00 
01 
02 
00 



1 2 3 4 5 6 7 



I I * 0.1591E-I04 
I 2 =| 0.1827E 01 

9 1 1 2 3l5 638^j0 lO«1827E9i01 

I 4 = 0.1827E 01 



23456789 



lib 



0123456789 



1 23456789 



01 23456789 



012345678901234 



012345678901234 
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ABSTRACT 



The purpose of this paper is to review the group of roadway design 
programs being utilized and developed by the U c S. Army Engineer Division, 
North Pacific, in cooperation with State Highway Departments of Washington, 
Oregon, and California, U. S. Bureau of Public Roads, U S. Forest Service 
and IBM. 

Programs discussed include horizontal and vertical control, data check, 
template design, alignment offset, quantity computation and several related 
utility programs, all of which are coded in SPS and make up an integrated 
system which is convenient as well as fast. The specific use of these 
programs by the Corps of Engineers is in conjunction with relocation of 
highways, railroads, county roads and municipalities involved in the 
design and construction of large multi-purpose hydroelectric projects. 

The "System" of programs has been developed through a combination 
of initial development, modification of others' programs and adoption of 
others' programs with no modification. Coordination of programming 
efforts is accomplished largely through facilities of HEEP. 

The abstracts of programs at the end of this paper are brief and 
general in nature. Writeups are complete but not necessarily in compliance 
with USERS GROUP Standards. Interested parties should address inquiries 
to: 

Division Engineer 

U e S. Army Engineer Division, North Pacific 
210 Custom House 
Portland 9, Oregon 
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PART I 



GENERAL PROBLEM 



Terrain * - The nature of the earthwork problems being encountered 
by the Walla Walla District may not be entirely unique, but they are 
certainly extreme at times. The Columbia and Snake Rivers have both 
cut deeply into Miocene basalt flows that form the regional bedrock. 
For U5 miles above the John Day Lock and Dam project, the Columbia River 
flows in a relatively deep valley with precipitous walls rising to a 
height of 3,000 feet above the general land surface on the north side, 
and up to 1,000 feet on the south side. The canyon walls consist of a 
series of stepped basalt cliffs interposed with steep talus slopes and 
terminating at the top in relatively flat plateau areas. Because of 
the constricted nature of the canyon, there is very little choice for 
the alignments; however, where the alignment is benched into the basalt 
edges, a minor shift in alignment can mean a major difference in cost 
because of the rock excavation and high fills which may occur. 

Because of the type and shape of the terrain, alignment is quite 
critical and many trial alignments may be required in order to arrive 
at what is considered an optimum one. 

Variety of Criteria. - Another feature that requires considerable 
attention is the number of agencies to be dealt with. In one project 
area there may be up to three State Highway Commissions, four railroads, 
and a dozen counties involved; all of which have from minor to major 
differences in criteria. This condition is emphasized when two rail- 
roads and a highway are aligned adjacently and there is hardly room 
for one of them in the canyon. Something has to move up the hill. 

The John Day Lock and Dam presents a good example of the magnitude 
of the work being done in the North Pacific Division. This will be 
a concrete gravity structure with an ultimate installation of 2,700,000 
KW produced by twenty 1359000KWI units driven by Kaplan turbines 
measuring 305 inches throat diameter. The total excavation for the 
project will be approximately 80,000,000 cubic yards, half of which 
will be hard rock. This total excavation figure is exceeded only by 
Owyhee, Fort Randall, Fort Peck, and Garrison project; all of which 
are earth or rockfill dams. As far as we know, the Uo, 000, 000 cubic 
yards of hard rock excavation will comprise a world record. 

The relocation at this project involves 80 miles of SP&S Railroad, 
57 miles of UP Railroad, kO miles of Washington State Highway, 32 miles 
of Oregon Interstate Highway, and the town of Arlington with a population 
of 633, I960 census. 

When the alignments involved are required to be revised up to 8 
or 10 times in order to reach a compromise with economics and the 
agencies concerned, the workload approaches staggering proportions. 

1 



PART II 



APPROACH TO PROBLEM 



Acquiring Data . - In order to plan and design such a project and 
its accompanying relocation work, surface data must be translated 
into punched cards representing ground cross sections along the align- 
ment centerline. 

First, of course, an approximate alignment must be determined on 
paper. The horizontal alignment and geometric computation programs 
are utilized at this stage for horizontal control. Survey parties then 
either "flag" that alignment in the field, or establish control points. 
Helicopters are utilized to transport the party members to inaccessible 
locations. Aerial photographers are then hired to "fly" the line and 
the resulting film is converted to glass plates. Sufficient area is 
covered to enable the designer to shift the alignment or position several 
roadway alignments such as railroad mainline and shoofly, and highway 
and detour. 

Data Translation and Preparation . - With the pictures taken and plates 
made, the next problem is to translate these to punched cards which are 
input media to the computer. This is accomplished by one of two ways, 
with emphasis on minimum manual work. For planning and design stages, 
if contour maps are available, alignments may be laid out on them and 
cross section cards punched by means of the Benson-Lehner Digital Scale, 
If contour maps are unavailable, or if for other reasons, it is desirable 
to do so, the cross section cards are prepared on a WILD A-T Autograph 
Stereo Plotter connected to an IBM 026 Key Punch through a WILD EK-5. 
A hard copy is prepared concurrently on an electric typewriter. An 
optional feature of this system permits printing out on typewriter the 
X, Y, and Z coordinates at points selected. 

For initial phases, cross sections are taken covering sufficient 
area to allow for alignment shifts and multiple alignments. 

Once the cross section cards are punched, they are checked for 
"detectable bulls"by means of a computer pass using the "Data Check" 
program which senses such errors as overhangs, double or no centerline, 
excessive horizontal or vertical difference between rod readings, etc. 
as may be logically determined. 

At this stage, cross sections can be plotted by the Benson-Lehner 
Electroplotter. To many designers, cross sections are considered 
obsolete but the terrain and materials with which we have to contend 
make cross sections desirable and often necessary. 

Another means by which segments of cross sections are sometimes 
prepared is the Benson-Lehner Oscar Trace Record Reader. This machine 
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uses fathometer graphs as input and produces punched cards in the 
same format produced by the Digital Scale or WILD A-7. This method is 
useful where the toe of proposed fill slopes extend below present water 
surfaces and in navigation channel and powerhouse excavation applications. 
A program is currently being developed by the Walla Walla District which 
will update cross section data cards by extending on either side, by 
inserting additional readings or replacing corrected readings. 
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PART III 



INTEGRATED PROCESSING 



Alignment and Realignment . - With the cross section cards prepared 
and checked, the question arises, "Has the alignment "been revised 
during the process of preparing the data cards?" If so, the new 
alignment computation data is combined with the original or previous 
alignment data and processed through the Earthwork Offset Alignment 
Reference program which computes the offset and skew angle "between 
the two alignments and new stationing* If the skew angle is not too 
great and the original cross section cards cover sufficient distance, 
the output cards of this computer pass serve as header cards which 
are combined with the cross section cards and processed through the 
Earthwork Alignment Shift program to compute new cross section cards 
referenced to the new alignment centerline. 

Template Generation . - Roadway templates to be used for computing 
cut and fill quantities can be prepared on the computer at this stage 
through what is at present a two-pass process. The first pass 
establishes profile grade and the coordinates of the basic roadway- 
points reflecting superelevation for horizontal curves. This basic 
roadway consists of two to four planes defined by from three to five 
points. 

The second pass completes the templates by adding ditches where 
required, and slope readings. This pass also serves as an additional 
data check by indicating slopes that will not catch, etc. The complexity 
of road design involving berms on fill slopes for embankment protection, 
and benches on cut slopes for rock fall protection or due to material 
classification change, has prompted a recent modification to this program 
which enables the designer to specify these features on the terminal 
slopes of the template. Templates can be produced with either slope 
readings or catch point coordinates at extremities. 

Quantity Computation . - For computing quantities, the cross section 
and template cards are merged and processed through the Earthwork 
Quantity Computation program. The output of the pass provides cut and 
fill areas, accumulated cut and fill quantities and mass ordinate. 
This program also has optional features for dredging and levee appli- 
cations, in which case, directed slopes and outside catch points are 
used. It will also make quantity adjustments for curvature correction. 

Design for Contract . - Listings of this output are provided the 
designer who then determines if and where alignment needs to be 
revised vertically or horizontally. This is repeated usually for 
various segments of the line, then when the alignment is satisfactory, 
the complete line is recomputed for quantities. 
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Special Options * - If two alignments are adjacent and the embank- 
ments sometime overlap, the Stage Development program is used to 
combine the cross section and template cards to produce new cross 
section cards which can then be line-shifted and used to compute 
berms for riprap which is to be computed separately from basic 
quantities. 

The answer cards from the quantity run can be processed through 
the Station Interval Quantity summary to produce listings of quantity 
summations at selected intervals. 

Cross section and template cards are processed through the 
Earthwork Data Plot Reduction program to replace template slope 
readings with catch point coordinates and reference all vertical 
coordinates to a common elevation so the cross sections and templates 
can be plotted. 

Pay Quantities Computation . - During and after construction when 
pay quantities are to be computed, the line is re-flown and "as- 
built" templates prepared. These cards differ from the previously 
mentioned templates in that they do not have slope readings but 
extend from catch point to catch point. These templates are merged 
with cross section cards of the final alignment and pay quantities 
are computed through the Earthwork Quantity Computation program. 

The "as-built" template cards are point-plotted over the "as- 
designed" cross sections and used to check for over- or under-built 
conditions. 
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PART IV 



RESULTS 



Opinions of the users of these programs are varied, so in order 
to avoid giving an overly optimistic view, this example was taken 
from the most "conservative" user. Sixty miles of railroad alignment 
on the Little Goose Lock and Dam project was processed through eight 
alignments before arriving at the final choice. By thoroughly 
analyzing the first seven alignments, five million yards of rock 
excavation were saved between the seventh and eighth alignments. If 
hand methods had been used, there would not have been time for more 
than four "shotgun" alignments and the savings resulting from computer 
utilization would not have been possible. 

The cost involved in design stages may be equal to or even 
greater than that incurred by "shotgun" manual methods, but this 
includes increased analyzation of alignments, making possible savings 
such as cited in the example. Even greater savings are being realized 
with a computer on-site in the Walla Walla District so that even more 
alignment trials can be made in certain areas where the alignment is 
most problematical. 

One installation has combined template design and quantity compu- 
tation programs enabling them to obtain quantities in one pass without 
producing intermediate templates. This version also incorporates 
line shifting as an integral part of the single pass. It requires 
a different format of data cards than that being used by our District 
and 60 K. The type of processing usually encountered in the Walla Walla 
District require some insertion of hand-prepared templates before the 
final pass, so we have not seriously considered their degree of 
combination. The two passes of template design might be combined to 
advantage, providing header card storage did not become too restricted. 
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PART V 

ABSTRACTS OF PROGRAMS 

1. HORIZONTAL ALIGNMENT - This program computes horizontal 
alignment from the following "basic input alignment data: 

a. Beginning station and coordinates (RP) 

b. Coordinates of PI*s 

c. Degree or radius of curves 

d. Spiral lengths (highways) or chord lengths and number 

(UPRR) 

Output information includes: 

a. Bearing and length of course "between Pi's 

b. Stations and coordinates of all curve points 

c. Deflection angle (^) and semi-tangents 

d. Degree, length and radius of circular curves 

e. Deflection angle (DE), x and y, u and v of spirals 

Curve definition can be by arc or chord. Distances are rounded 
to nearest .01 foot, coordinates are carried to nearest .001 foot. 
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2. GEOMETRIC COMPUTATION - The present geometric program 
combines the features of several related programs which were used on 
the IBM 650. The following types of computations can be performed 
singly or in combination in one machine pass. 

a. Survey traverse using bearings with azimuth or deflection 

angles « 



b. Compass or transit rule traverse adjustment. 



c. Areas bounded by traverses which may contain circular 

segments. 

d. In a traverse between known points, a combination of two 
unknowns (bearings, 2 distances or a distance and a bearing) may be 
solved for and stored for recall in subsequent problems in same pass 
with interdependency of geometric figures. 

e. Location of points on tangents, circular curves, spirals 
or offset spirals and intersections of any two of these. 

f • The intersection of tangents or circular curves with 
spirals or offset spirals. 

g. Convert selected coordinates from any one plane system 
to another. 

The ability to store data pertaining to any course for use in 
problems which follow makes it possible for the engineer to begin 
with a minimum of data and to solve complex office and field layout 
problems in one pass on the computer. 

Problems are arranged in logical sequence to develop intermediate 
answers as would be done by "people-computers" following traditional 
methods • 



A closed traverse can be adjusted in one problem, for instance, 
and an adjusted loop can be run from it to close on any point for 
which coordinates are known or have been computed. Two alignments can 
be intersected and the bearing and distance to land corners or other 
references readily determined. The area of gusset plates can be 
determined and a host of other geometric problems can be solved. As 
a matter of fact, some users say"it can do anythingl" 
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3. EARTHWORK DATA CHECK 

This program, checks earthwork input data in format used by 
U. S. Army Engineer Division, North Pacific, and its various Districts, 
Errors which would result in erroneous results or machine halts are 
detected and error messages typed or punched identifying the station 
and specifying the type of error detected. 

U. PROFILE GRADE AND TEMPLATE GENERATION 

This program converted from Washington State Highway Department^ 
Five-Point Profile Grade Program, produces the following results: 

a. Generates basic roadbed templates for input to Slope 
Selection - Template Generation Program. 

b. Produces pavement elevation cards supplying elevation of 
five points at each station to .01 foot accuracy. 

c. Computes profile grade at selected or incremental stations. 

d. Reproduces templates inserting new profile grade elevation 
in the reference elevation field. 

e. Punches a summary of vertical control data at end of each 

run. 

5. SLOPE SELECTION - TEMPLATE GENERATION 

This program requires ground line cross section cards plus 
the template cards produced by the Profile Grade Program as basic data. 
Other input consists of slope specification cards and bench-berm data 
cards. Three fill slopes and four cut slopes may be specified for 
appropriate depths of fill or cut. Ditches can be specified with V 
or flat bottom with depth and slope optional. Terminal slopes may 
include bench (cut) or berm (fill). 

The templates produced may include catch points, or the 
terminal rod readings may be replaced with slope readings in the form 
xxx .xx horizontal to 1 foot vertical. 

These templates when collated with ground line cross section 
cards are the input to the Earthwork Quantity Computation Program. 
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6. EARTHWORK QUANTITY COMPUTATIONS 

a. The purpose of this program is to compute quantities of 
earthwork to "be excavated and placed on a given job. The program is 
patterned after the IBM 650 Cut and Fill Program (H-8U1) as modified 
by Bureau of Public Roads, Vancouver, Washington. The computations 
are performed in the conventional manner by taking the original ground 
topography and the design of the completed section, as specified by 
the engineer, and computing the cut and fill end areas. These are 
then used to compute the volumes between sections by the Average-End- 
Area method. This program will apply a volume correction due to 
curvature if desired. Swelling or compaction can be applied to either 
the cut or the fill and these quantities are accumulated throughout 

a project to produce a mass-ordinate. The total accumulated quantities 
of unadjusted cut and fill are also available for each section. 

b. As a by-product of these computations, some design infor- 
mation is available at each section to aid the engineer in his appraisal 
of the design that has just been computed. These include the cut or 
fill at the pivot point, at the toe of both the right and left slope, 
and at centerline. Also, the rod and distance to the catch point of 
the template slope and the original ground line, if these were not 

part of the input data. 

7. EARTHWORK LINE SHIFT 

The purpose of this program is to shift topog card data 
(type "0") to a new or offset centerline. This enables the engineer 
the choice of running multiple trial alignments while using the original 
topog data but shifting the topog data to coincide with the alignment. 
The program can also be used for station adjustment and for the reverse 
line shifting of topog data. It is desirable in some phases of earth- 
work to shift the topog data from the original line to an offset line, 
generate templates at the offset line, stage the ground line, and then 
shift the staged ground back to the original line. By proper program 
control in the Job Control Card, the original line shift headers (used 
to shift original topog to an offset line) and the staged ground at 
the offset line can be used to shift the staged ground back to the 
original line. 

8. EARTHWORK TEMPLATE SHIFT 

This program is designed to take the Type "l" templates that 
were prepared or generated at an offset centerline and punch new templates 
having original centerline stationing and shot distances referenced to 
the original centerline. The new Type M l" templates can then be collated 
with the original terrain cards (Type "0") and a plotter deck is obtained 
which has the original centerline as a base for plotting of both terrain 
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data and Type "l" templates. Use of this program will facilitate 
plotting operations when more than one roadway section is planned and 
is to be plotted on the same cross section of terrain. 

9. EARTHWORK ALIGNMENT OFFSET 

This program is designed to compute offset distances, skew 
angles, and centerline station coordinates of a base line and an 
offset line. Output also includes line shift headers which enables 
the Earthwork Line Shift Program to produce a new set of ground line 
cards referenced to the offset line. 

10. EARTHWORK DATA PLOT REDUCTION 

This program reproduces ground line and template cards adding 
plotter signs where necessary and references both type of cards to a 
common reference elevation. It will also provide demand origin offset 
codes and insert blank cards where necessary. 

11. CENTERLINE DATA PLOT REDUCTION 

This program compiles station numbers and centerline elevations 
of ground line cross sections or template cards in a format for plotting 
of centerline profile. 

12. STATION INTERVAL QUANTITY SUMMARY 

The purpose of this program is to produce a summary of quantities 
of earthwork to be excavated and/or placed on a given job. Two types 
of summary options are available: Station or Interval. The listings 
are used by the engineer for design analyses and as tabulated quantity 
listings for plan and profile drawings for design memorandums and/or 
contract drawings. 

13. PLANE COORDINATE CONVERSION - PHOTOGRAMMETRY LEAST SQUARES METHOD 

The purpose of this program is primarily to convert the machine 
coordinates of a universal stereo-plotter to a local coordinate system, 
applying the least squares method of conversion (Helmert). 

This conversion is mathematically expressed as follows: 

x = p + k X cos a - k Y sin a 

y = q + k Y cos a + k X sin a 

The quantities p and q are ground coordinates of the origin of the 
machine coordinate system, a is the angle of rotation between the two 
systems, and k is the distance scale factor. 
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Charles R. Hebble 1/ 
INTRODUCTION 

The computer program described in this paper, along with certain 
fixed input data, constitutes a detailed mathematical model of a system 
of hydroelectric projects. It is a general program applicable to any 
river system and scheme of development. Natural lake and channel storage 
may be synthesized in addition to reservoir storage and backwater effects. 
The program is capable of accurately simulating the hour-by-hour power 
loading and water regulation of a hydro system. 

Through use of the program it is possible to determine the operating 
characteristics of planned future projects for various- sized power 
installations. Backwater encroachment on upstream reservoirs, pondage 
requirements, and the effects of peaking discharges on downstream river 
stages and reservoirs may be evaluated. 

The program may also be used to compare alternative distributions 
of system load among an existing group of hydraulically and electrically 
integrated projects. It is thus possible, through a trial and error 
approach, to determine optimum load distribution. 

The program has been used in investigating existing system develop- 
ment, along with future planned developments of the lower Columbia and 
lower Snake Rivers. Here a series of run-of-river projects, below 
Grand Coulee Dam, will develop almost all of the available head (see 
Figures 1 and 2). It is anticipated the program will be used to analyse 
future operating problems in regards to power loading and water regulation. 

DESCRIPTION OF EQUIPMENT 

Two versions of the program exist: one for an IBM 1620 having a 
40,000 digit memory with an IBM 1622 Card READ- PUNCH for input- output ; 
the other for an IBM 1920 System (combined IBM 1620/1401) having a 
60,000 digit memory in the IBM 1620 and 4,000 digits in the IBM 1401. 
This latter system has both an IBM 1402 Card READ-PUNCH and an IBM 1403 
on-line printer for input-output. Both systems have the following 
optional features required for program execution: Transfer Numeric 
Strip, Transfer Numeric Fill, Move Flag, and Indirect Addressing. 
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Peripheral equipment required is as follows: card keypunch, card 
verifier, and card sorter. In addition, for the system with the IBM 
1622, a card tabulating machine with a special wired panel is required. 

ANALYSIS AND METH OD OF SOLUTION 

General . - The program is intended for the detailed hour-by-hour 
analysis of hydro system operations; basic input to the program and 
output from the program is hourly. However, periods greater or less 
than one hour may be used. There is no limit to the number of real- 
time operating intervals which can be simulated by the computer. The 
program continues to run as long as there are input data in the READ 
hopper. Ordinarily, however, real-time periods of more than a week 
are not analyzed because of the excessive length of the computer runs 
for all but the smallest systems. As a rough estimate, the program 
requires two minutes to compute one day's output for each hydro project 
in the system. The number of reaches and projects that can be analyzed 
in a given system is dependent on the computer storage capacity. The 
program instructions require approximately 25,000 digits of memory. A 
rule-of- thumb estimate for a given hydro system, including open river 
reaches, is 1,000 digits of memory for each hydro station. Variation 
of memory requirement is due to optional lengths of tables used to 
define system parameters. 

Program Operating Modes . - Discharges from the various hydro projects 
are specified, either directly or indirectly as input to the program for 
each real-time period. Four operating modes exist for this purpose: 
(1) station power loading given; (2) system power load plus generation 
allocation (breakpoint) settings given; (3) total project discharge 
given; and (4) project forebay elevations specified. Any of the fore- 
going, in combination with other fixed and variable data, determines the 
individual project discharges. Different projects may be operated 
under different operating modes at the same time, and the operating mode 
of any project may be changed at any time. 

Variable Input . - Each hydro station obtains operating data for 
each real-time period from coded variable input cards. Variable input 
data include project power loads (or discharges), system loads with 
project breakpoint settings for load- frequency control, local inflow, 
optional spill, and number of generating units synchronized on the line 
and available. These data may be varied as desired during the course 
of a run. 

Operation Sequence . - The program begins with the upstream project 
of the series and proceeds downstream, routing flows through each open 
channel reach or reservoir. Local inflow between projects is added to 
the routed flows. The outflow from the projects is computed as the sum 
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of power discharge, spill, and average fixed release. (The average 
fixed release, which must remain constant throughout the run, includes 
losses due to lockage, leakage, and useage.) The entire process is 
repeated for successive 15-minute computation (routing) periods. 
Project and desired reach answers are obtained upon completion of the 
fourth routing period. The variable data for the next hour is then read 
in and the entire process repeated. Routing and output intervals may 
be varied by a special job definition routine. The routing interval 
may be different for different projects or river reaches but must be 
a multiple of the basic routing interval. The computer program has 
two fundamental parts. One is the method of simulating a hydro-power 
station, the other is the flow routing procedure. Each part will be 
considered separately and in some detail. 

Hydro-Power Station Simulation . - Each generating station synthe- 
sized in the program is represented by individual unit characteristics 
as illustrated in Figure 3. This is in contrast to the method used in 
other power programs for coarser increments of time which consider only 
generating station characteristics as a whole. Such programs assume 
some operating efficiency between best efficiency and, say, full-gate 
efficiency. This assumption, while adequate for studies having a basic 
time period of day, week, or month, during which there are many swings 
in station power loadings, are inadequate for the detailed hourly com- 
putations required in the study of peaking operations. The number of 
generating units synchronized on the line may be specified as variable 
input or may be automatically selected by the program. The program has 
the ability to automatically add units as required to meet the load up 
to the maximum number available. When a power load in excess of the 
maximum capability of the installation is specified, generation will 
automatically be adjusted to equal plant capability as limited by head. 
The output is coded to indicate this alteration by the computer program. 

Under automatic selection of units the program computes the number 
of units for best efficiency operation. The total station load (or 
discharge) is divided by the unit best-efficiency loading (or discharge) 
for the particular head existing on the station to determine the desired 
number of units. Fractional numbers of units are truncated; thus unit 
loadings are at or higher than the best efficiency point. The number of 
units selected is limited to the maximum available. 

The program automatically causes all excess water above maximum 
pool elevation to be spilled. This is termed "mandatory" spill. Power 
discharge, fixed release, optional spill and mandatory spill are summed 
to arrive at the total project discharge. The alarm section of the out- 
put is coded calling attention to such mandatory spill. 

Should a project attempt to draft below its minimum pool elevation 
as a result of releases exceeding inflow; the power discharge and hence, 



3 




4 



V 



generation is automatically reduced to prevent the overdraft, 

irrespective of the power demand or discharge called for by the input , 
data. Here again the alarm section of the output is coded. 1 

The computation of power, discharge, head, etc. are accomplished 
by means of table look-up and interpolation. For existing projects, 1 
unit performance tables may be prepared from observed data. In the 

case of planned projects, the basic data must first be computed. A ' 
separate computer program has been developed which can quickly compute 

performance characteristics. The program is based on the turbine < 

performance characteristics of a unit machine, i.e., a machine runner 

diameter of one foot operating under a net head of one foot. These 

data are summarized in a unit performance hill in which power output 

is plotted against the peripheral speed coefficient. It is thus 

possible to quickly evaluate the performance of different sized units 

and units having different characteristics. 

Centralized Load- Frequency Control . - A feature of the program is 
a provision for simulating centralized control of the station power 
loadings in a manner similar to that of centralized load- frequency 
control equipment. Such equipment is presently installed at the major 
Federal hydro stations on the Columbia River. The power loadings of 
these projects are centrally controlled from the system load dispatcher's 
office in Portland, Oregon. 

When the load- frequency control feature of the program is used, 
participating settings for each of the hydro stations are then included 
as input data and a single system load rather than individual station 
loadings are specified. The system load is apportioned among the hydro 
stations in accordance with station participation values. Other than 
this, the program operates in the same manner as when individual loadings 
are prescribed for the several hydro stations. 

This centralized, automatic dispatching of power loadings affords 
the opportunity for efficient coordinated system Operation of the Federal 
projects. Use of the load- frequency control feature of the computer 
program will allow planning- studies and scheduling of power operations 
in conjunction with the equipment, in addition to its other uses. It is 
envisioned that ultimately the load dispatching of these hydro stations 
will be directly and automatically controlled by computer. The present 
program could well serve as a basis for such a future program. 

Flow Routing . - The total discharge from a project (power discharge, 
spill, and fixed release) is routed downstream by a method of flow routing 
known as incremental storage routing. This approach consists of sub- 
dividing reservoirs or river channels into incremental reaches. Each 
incremental reach is represented in the program by two tables: one giving 
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the relationship between storage and elevati on for the particular reach; 
the other giving the relationship between discharge and elevation . By 
proper selection of routing interval and choice of the number of reaches, 
actual hydraulic conditions may be accurately approximated. The discharge 
from any reach is computed as a function of the reach elevation, or where 
backwater effects must be considered, as a function of the next down- 
stream reach as well. The change in storage content of the incremental 
reservoirs is related to the difference between inflow and outflow. 

The combined effect of the storage-elevation and discharge-elevation 
relationships mentioned in the preceding paragraph causes a time constant 
to be introduced into the routing. This is known as "time of storage" 
and is defined as the change in storage per unit change in discharge . 

Ts = A S/ A Z = A S 
A Q/ A Z A Q 

where , 

Ts = time of storage 
S = storage volume 
Q = discharge rate 
Z = elevation 

Dimensionally , using the units adopted for program use, 

Ts (min) = A S(kcfs-m i n)/ A Z(feet) 
A Q(kcfs)/ A Z(feet) 

The routing procedure used in the program does not explicitly consider 
time of storage; however, the routing interval used in the program must 
be shorter than the minimum time of storage for any given reach. Failure 
to observe this criterion will result in oscillations of the reach eleva- 
tion and discharge. These oscillations tend to increase in amplitude 
until the range of tables is exceeded. A reservoir which is encroached 
upon by backwater from a downstream reservoir may have an extremely short 
time of storage. 

The time- of- storage concept may be visualized by assuming a constant 
rate of increase in inflow to an incremental reservoir. Outflow from the 
reservoir will eventually reach an equilibrium condition where it is 
increasing at the same rate as the inflow; at this time, the outflow 
hydrograph will be displaced from the inflow hydrograph by the time of 
storage. Figure 4 illustrates a discharge hydrograph routed through six 
identical incremental reservoirs. Storage-elevation and discharge- 
elevation curves are shown on this same sheet. 
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Overlapping reservoirs are normally subdivided into two reaches: 
a forebay reach and a tailwater reach for the next upstream project. 
The sum of the storages of the individual reaches must equal the actual 
reservoir storage. 

In subdividing open channel reaches into incremental reaches, no 
definite rules can be given for the number of increments. In general, 
the greater the number of reaches, the greater the translation of the 
discharge hydrograph with minimum attenuation. Discharge hydrographs 
can be modified, to a lesser degree, by choice of routing interval. 
Here again, selection of routing interval beyond that required to pre- 
vent oscillation, is subject to rules established for particular problem 
definition. 

For river systems where all tailwater reaches have extremely short 
times of storages, an alternate routing method is available. This 
method permits use of relatively long routing intervals, up to the answer 
output interval if desired. In effect it considers the time of storage 
to be zero. Tailwater elevation is first determined as a function of 
inflow (and next downstream reach elevation if necessary); the storage 
corresponding to the new tailwater elevation is then computed as is the 
change in storage during the routing interval; this change in storage 
is converted to an average rate of flow during the routing interval and 
added (or subtracted) to tailwater inflow to arrive at the outflow. 
The algebraic sum added to the inflow is dependent on change of tailwater 
elevation during the routing period, i.e., a rising tailwater will cause 
outflow to be less than inflow with the opposite occurring on a declining 
tailwater condition. 

The flow and storage- routing procedure used in the program has a 
sound theoretical basis which permits analysis of future reservoirs as 
well as existing reservoirs and river reaches. The computation of 
discharge, elevation, storage content, etc. are accomplished by means 
of table look-up and interpolation. Volume- elevation relationships for 
open river reaches and reservoirs are generally obtained from survey 
data. Water surface profiles are computed for various forebay elevations 
and rates of flow by independent means. (Here again separate computer 
programs are available for determining water surface profiles, and 
volume relationship, in reservoirs or open channel reaches.) From these 
profiles, the discharge- elevation- volume relationships are derived. 



Output Data . - Program output consists of the following data for 
each project: 



PROGRAM OUTPUT 



Date and time 
Abbreviated project name 
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End- of- period forebay elevation 

End-of- period tailwater elevation 

Average head 

Forebay storage content 

Storage change during period 

Local inflow 

Routed inflow 

Power discharge 

Spill 

Total discharge 
Generation 

Number of units operating 
Alarm codes 
Operating mode number 
Project number 

Output may be obtained for selected reaches, in addition to project 
data, at users option. Reach output consists of the following data: 

Date and time 
Abbreviated reach name 
End- of- period elevation 
End-of-period total discharge 
Reach number 

The 1920 System program output consists of both on-line printing 
and punched cards. The same data are given by both forms of output. 
Choice of output mode is optional with user. The 1620 System output is 
punched cards. The 1920 System uses an auxiliary listing program to list 
the punched cards for each project or optional reach and system Summary 
Cards . 

Output data consists of Project Cards, selected Reach Cards and 
Summary Cards containing end-of-period system generation totals and 
system potential energy remaining in storage. The sequence of the 
punched card output and on-line printout is as follows: An initial 
Summary Card gives potential energy in the system based on starting 
profile data. This is followed by data for each hydro project and for 
optional selected reaches for the initial time period. A Summary Card 
concludes the data for the initial time period. Project, Reach and 
Summary data are then punched out for the next subsequent output interval. 
Output intervals may be varied by special control. In the absence of 
special control, the output interval is one hour. Sample output formats 
of a typical study are shown on Figure 5 for the IBM 1620 System and on 
Figure 6 for the IBM 1920 System. 
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36*650 






500 






2 


1 . 506 




36*650 






500 






2 


1 4 * *) 

1.268 




36.620 


5*690 




6* 190 


133 


2.0 


2 


731 




36.650 


16*600 




19* 100 


439 


4*0 


2 


594- 




36.620 


50*370 




50*870 


1.197 


13*0 


2 


1 .602- 




36.650 


74*590 




75.090 


1,749 


20.0 


2 


1.991- 




36.650 


63*930 




84*430 


1,944 


22.7 


2 


2.015- 




36.620 


64*490 




84,990 


1.944 


23.0 


2 


2.026- 




36*650 


84*780 




85.280 


1.944 


23.0 


2 


2.035- 




36*620 


84*970 




85.470 


1 .944 


23.0 


2 


1.677- 




36.650 


76.380 




76.880 


1,749 


21.0 


2 


1.364- 




36.650 


68*890 




69.390 


1.580 


18.5 


2 


1.262- 




36.620 


66*420 




66.920 


1.527 


18.0 


2 


1.460- 




36*650 


71*190 




71.690 


1.638 


19.0 


2 


2.022- 




36*620 


84*650 




85.150 


1.944 


23.0 


2 


2.034- 




36*650 


84*970 




85*470 


1,944 


23.0 


2 


2.042- 




36*650 


85*160 




85*660 


1,944 


23.0 


2 


1*880- 




36.620 


81*240 




81.740 


1,853 


22.0 


2 


1*374- 




36*650 


69*130 




69,630 


1.580 


19.0 


2 


562- 




36.620 


49,610 




50,110 


1.139 


13*0 


2 


242 




36.650 


30*340 




30.840 


704 


8*0 


2 


728 




36*650 


18.690 




19, 190 


439 


5*0 


2 


7*223- 




880*510 






29,513 







22866*42 



760.474 



1.281.880 



1,293,880 



346*2 



•n 
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ELfcV. tNO PER 100 


AVG. 


FORE BAY 


DELTA 


••••••••• INFLOW 


•••••••• 


••••••• 


DISCHARGE 




AVG. 


AVERAGE 


OP 










FOKEBAY 


TAILNTR 


HEAD 


STORAGE 


STORAGE 


LOCAL 


ROUTEO 


TOTAL 


POWER 


SPILL 


TOTAL 


GEN. 


NUMBER 


CD 


DATE AND 


TIME 


FEtT 


FEET 


FEET 


CFS-DYS 


CFS-DYS 


CFS 


CFS 


CFS 


CFS 


CFS 


CFS 


NW 


UNITS 


• 


15 


01 


85 


0100 


1262.08 


945.41 


315.30 


204 » 062 


1,319 





37,750 


37,750 


5,740 





6,240 


133 


2.0 


2 


lb 


01 


85 


0200 


1262.16 


944.47 


317.19 


205,485 


1,423 





36,700 


36,700 


2,050 





2,550 


45 


1.0 


2 


IS 


01 


85 


0300 


1262.24 


944.05 


317.94 


206,991 


1,506 





36,650 


36,650 








500 





0.0 


2 


Id 


01 


85 


0*00 


1262.33 


943.89 


318.27 


208,497 


1,506 





36,650 


36,650 








500 





0.0 


2 


li> 


01 


85 


0*00 


1262.40 


944.48 


318.02 


209,765 


1,268 





36,620 


36,620 


5,690 





6,190 


133 


2.0 


2 


15 


01 


85 


0600 


1262.44 


945.99 


317.05 


210,496 


731 





36,650 


36,650 


18,600 





19,100 


439 


4.0 


2 


15 


01 


85 


0700 


1262.41 


950.20 


314.43 


209,903 


594- 





36,620 


36,620 


50,370 





50,870 


1,197 


13.0 


i 


15 


01 


85 


0800 


1262.32 


953.79 


310.61 


208,301 


1,602- 





36,650 


36,650 


74,590 





75,090 


1,749 


20.0 


2 


IS 


01 


85 


0900 


1262.20 


955.99 


307.45 


206,310 


1,991- 





36,650 


36,650 


83,930 





84,430 


1,944 


22.7 


4. 


15 


01 


85 


1000 


1262.09 


957.01 


305.73 


204,295 


2,015- 





36,620 


36,620 


84,490 





84,990 


1,944 


23.0 


2 


1* 


01 


85 


1 100 


1261.97 


957.57 


304.79 


202,268 


2,026- 


o 


4o,o>U 


36,650 


84,780 





85,280 


1 ,944 


23.0 


2 


15 


01 


85 


1200 


1261.86 


957.89 






C |UJJ 





36,620 


J W , M&V 


84, 970 


o 


85,470 


1 • 944 


23.0 


2 


lb 


01 


85 


1300 


1261.76 


957.38 


304.13 


198,556 


1.677- 







36,650 


76,380 





76,880 


1,749 


21.0 


2 


15 


01 


85 


1400 


1261.68 


956.50 




IOT |Q) 


i . 

l , JOt* 





36,650 


Jvf 


WO , 9 w V 


o 


69, 390 


1 , 580 


18.5 




15 


01 


85 


1500 


1261.61 


955.83 


305.43 


195,930 


1,262- 





36,620 


36,620 


66,420 





66,920 


1,527 


18.0 


2 


15 


01 


85 


1600 


1261.53 


955.84 


305.74 


194,470 


1,460- 





36,650 


36,650 


71,190 





71,690 


1,638 


19.0 


2 


15 


01 


85 


1700 


1261.41 


956.91 


305.18 


192,448 


2,022- 





36,620 


36,620 


84,650 





85,150 


1,944 


23.0 


2 


15 


01 


85 


1800 


1261.29 


957.52 


304.19 


190,413 


2,034- 





36,650 


36,650 


84,970 





85,470 


1,944 


23.0 


2 


15 


01 


85 


1900 


1261.18 


957.85 


303.58 


188,371 


2,042- 





36,650 


36,650 


85,160 





85,660 


1,944 


23.0 


2 


15 


01 


85 


2000 


1261.07 


957.72 


303.33 


186,491 


1,880- 





36,620 


36,620 


81,240 





81,740 


1,853 


22.0 


2 


15 


01 


85 


2100 


1260.99 


956.69 


303.74 


185,117 


1,374- 





36,650 


36,650 


69,130 





69,630 


1,580 


19.0 


2 


15 


01 


85 


2200 


1260.96 


954.43 


305.17 


184,555 


562- 





36,620 


36,620 


49,610 





50,110 


1,139 


13.0 


2 


15 


01 


85 


2300 


1260.97 


951.15 


307.96 


184,797 


242 





36,650 


36,650 


30, 340 





30,840 


704 


8.0 


2 


15 


01 


85 


2400 


1261.02 


947.86 


311.19 


185,524 


728 





36,650 


36,650 


18,690 





19, 190 


439 


5.0 


2 




TOTALS 










17,223- 














29,513 








AVERAGES 


1261.75 


952.77 


308.97 


198,353 







36,690 


36,690 


53,410 





53,910 


1,230 


14.4 






\ J 



